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Abstract 

This review concerns the origin and the possible effects of magnetic fields in 
the early Universe. We start by providing to the reader with a short overview 
of the current state of art of observations of cosmic magnetic fields. We then 
illustrate the arguments in favour of a primordial origin of magnetic fields in 
the galaxies and in the clusters of galaxies. We argue that the most promising 
way to test this hypothesis is to look for possible imprints of magnetic fields 
on the temperature and polarization anisotropies of the cosmic microwave 
background radiation (CMBR). With this purpose in mind, we provide a 
review of the most relevant effects of magnetic fields on the CMBR. A long 
chapter of this review is dedicated to particle physics inspired models which 
predict the generation of magnetic fields during the early Universe evolution. 
Although it is still unclear if any of these models can really explain the 
origin of galactic and intergalactic magnetic fields, we show that interesting 
effects may arise anyhow. Among these effects, we discuss the consequences 
of strong magnetic fields on the big-bang nucleosynthesis, on the masses and 
couplings of the matter constituents, on the electroweak phase transition, 
and on the baryon and lepton number violating sphaleron processes. Several 
intriguing common aspects, and possible interplay, of magnetogenesis and 
baryogenesis are also dicussed. 
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Introduction 

Magnetic fields are pervading. Planets, stars, galaxies and clusters of galaxies 
have been observed that carry fields that are large and extensive. Though 
strong homogeneous fields are ruled out by the uniformity of the cosmic 
background radiation, large domains with uniform fields are possible. 

A crucial ingredient for the survival of magnetic fields on astrophysical 
scales is for them to live in a medium with a high electrical conductivity. As 
we shall see in Chap.|^, this condition is comfortably fulfilled for the cosmic 
medium during most of the evolution of the Universe. As a consequence, it 
is possible for magnetic fields generated during the big-bang or later to have 
survived until today as a relic. 

To establish the existence and properties of primeval magnetic fields 
would be of extreme importance for cosmology. Magnetic fields may have 
affected a number of relevant processes which took place in the early Uni- 
verse as well as the Universe geometry itself. Because of the Universe high 
conductivity, two important quantities are almost conserved during Universe 



evolution: the magnetic flux and the magnetic helicity (see Sec. p.. 41) . As we 
will see, valuable information about fundamental physics which took place 
before the recombination time may be encoded in these quantities. 

In the past years a considerable amount of work has been done about 
cosmic magnetic fields both on the astrophysical and the particle physics 
point of view. The main motivations of such wide interest are the following. 

The origin of the magnetic fields observed in the galaxies and in the clusters 
of galaxies is unknown. This is an outstanding problem in modern cosmology 
and, historically, it was the first motivation to look for a primordial origin 
of magnetic fields. Some elaborated magnetohydrodynamical (MHD) mech- 
anisms have been proposed to amplify very weak magnetic fields into the fiG 
fields generally observed in galaxies (see Sec. |l.l| ). These mechanisms, known 



as galactic dynamo, are based on the conversion of the kinetic energy of the 
turbulent motion of the conductive interstellar medium into magnetic energy. 
Today, the efficiency of such a kind of MHD engines has been put in question 
both by improved theoretical work and new observations of magnetic fields 



in high redshift galaxies (see Sec. |1.2| ). As a consequence, the mechanism re- 
sponsible for the origin of galactic magnetic fields has probably to be looked 
back in the remote past, at least at a time comparable to that of galaxy for- 
mation. Furthermore, even if the galactic dynamo was effective, the origin 



of the seed fields which initiated the processes has still to be identified. 

Even more mysterious is the origin of magnetic fields in galaxy clusters. 
These fields have been observed to have strength and coherence size compa- 
rable to, and in some cases larger than, galactic fields. In the standard cold 
dark matter (CDM) scenario of structure formation clusters form by aggrega- 
tion of galaxies. It is now understood that magnetic fields in the inter-cluster 
medium (ICM) cannot form from ejection of the galactic fields (see Sec. [1.2|) . 



Therefore, a common astrophysical origin of both types of fields seems to 
be excluded. Although, independent astrophysical mechanisms have been 
proposed for the generation of magnetic fields in galaxies and clusters, a 
more economical, and conceptually satisfying solution would be to look for 
a common cosmological origin. 

Magnetic fields could have played a significant role in structure formation. It 
may be not a coincidence that primordial magnetic fields as those required 
to explain galactic fields, without having to appeal to a MHD amplifica- 
tion, would also produce pre-recombination density perturbations on proto- 
galactic scales. These effects go in the right direction to solve one of the 
major problems of the CDM scenario of structure formation (see Sec. [1. 31 ). 
Furthermore, if primordial magnetic fields affected structure formation they 
also probably left detectable imprints in the temperature and polarization 
anisotropies, or the thermal spectrum, of the cosmic microwave background 
radiation (CMBR) (see Chap.l). 

Field theory provides several clues about the physical mechanisms which 
may have produced magnetic fields in the early Universe. Typically, magne- 
togenesis requires an out-of-thermal equilibrium condition and a macroscopic 
parity violation. These conditions could be naturally provided by those phase 
transitions which presumably took place during the big-bang. Well known 
examples are the QCD (or quark confinement) phase transition, the elec- 
troweak (EW) phase transition, the GUT phase transition. During these 
transitions magnetic fields can be either generated by the turbulent motion 
induced in the ambient plasma by the rapid variation of some thermodynamic 
quantities (if the transition is first order) or by the dynamics of the Higgs 
and gauge fields. In the latter case the mechanism leading to magnetogene- 
sis shares some interesting common aspects with the mechanism which have 
been proposed for the formation of topological defects. On the other hand, if 
cosmic strings were produced in the early Universe they could also generate 



cosmic magnetic fields in several ways. Inflation, which provides a consistent 
solution to many cosmological puzzles, has also several features which make 



it interesting in the present context (see Sec . [4 .51 ). Although to implement an 
inflationary scenario of magnetogenesis requires some nontrivial extensions 
of the particle physics standard model, recent independent developments in 
field theory may provide the required ingredients. Magnetic fields might also 
be produced by a preexisting lepton asymmetry by means of the Abelian 
anomaly (see Sec. [4.4|) . Since the predictions about the strength and the spa- 
tial distribution of the magnetic fields are different for different models, the 
possible detection of primeval magnetic fields may shed light on fundamental 
physical processes which could, otherwise, be unaccessible. 

Even if primordial magnetic fields did not produce any relevant effect after 
recombination, they may still have played a significant role in several funda- 
mental processes which occurred in the first 100,000 years. For example, we 
shall show that magnetic fields may have affected the big-bang nucleosynthe- 
sis, the dynamics of some phase transitions, and baryogenesis. Since big-bang 
nucleosynthesis (BBN) has been often used to derive constraints on cosmo- 
logical and particle physics parameters, the reader may be not surprised to 
learn here that BBN also provides interesting limits on the strength of pri- 
mordial magnetic fields (see Chap.|^). Even more interesting is the interplay 
which may exist between baryogenesis and magnetogenesis. Magnetic fields 
might have influenced baryogenesis either by affecting the dynamics of the 
electroweak phase transition or by changing the rate of baryon number vi- 
olating sphaleron processes (see Chap.|^). Another intriguing possibility is 
that the hypercharge component of primeval magnetic fields possessed a net 
helicity (Chern-Simon number) which may have been converted into baryons 
and leptons by the Abelian anomaly (see Chap.^. In other words, primor- 
dial magnetic fields may provide a novel scenario for the production of the 
observed matter- antimatter asymmetry of the Universe. 

An interesting aspect of magnetic fields is their effect on the constituents 
of matter. This in turn is of importance on many aspects of the processes 
that took place in the early times. Masses of hadrons get changed so that 
protons are heavier than neutrons. The very nature of chirality could get 
changed see (Chap.|^). However the characteristic field for this to happen 
is H = rra^ which is about 10^^ G. These fields cannot exist at times when 
hadrons are already existing and therefore are probably not relevant. Near 



cosmic superconductive strings the story may be different. 

Clearly, this is a quite rich and interdisciplinary subject and we will not be 
able to cover with the same accuracy all its different aspects. Our review 
is manly focused on the production mechanism and the effects of magnetic 
fields before, or during, the photon decoupling from matter. 

In Chap.[^ we shortly review the current status of the observations. In 
order to establish some relation between recent time and primeval magnetic 
fields we also provide a short description of some of the mechanisms which 
are supposed to control the evolution of magnetic fields in the galaxies and 
in the intergalactic medium. We only give a very short and incomplete 
description of the effect of magnetic fields on structure formation. Some basic 
aspects of this subject are, anyhow, presented in Chap.^ where we discuss 
the effect of magnetic fields on the anisotropics of the cosmic microwave 
background radiation. From a phenomenological point of view Chap.^ is 
certainly the most interesting of our review. The rapid determination of 
the CMBR acoustic peaks at the level of a few percent will constrain these 
fields significantly. We briefly touch upon the recent determination of the 
second acoustic peak. In Chap .^ we describe several effects of strong magnetic 
fields on the BBN and present some constraints which can be derived by 
comparing the theoretical predictions of the light elements relic abundances 
with observations. Since it can be of some relevance for BBN, propagation 
of neutrinos in magnetized media is also shortly discussed at the end of that 
chapter. In Chap.^ we review several models which predict the generation 
of magnetic fields in the early Universe. In the same chapter some possible 
mutual effects of magnetogenesis and baryogenesis are also discussed. Some 
aspects of the effects which are described in Chapts.|^ and ^, which concern 
the stability of strong magnetic fields and the effect that they may produce 
on matter and gauge fields, are discussed in more details in Chap.|^. At the 
end we report our conclusions. 



Chapter 1 

The recent history of cosmic 
magnetic fields 

1.1 Observations 

The main observational tracers of galactic and extra-galactic magnetic fields 
are (comprehensive reviews of the subject can be found in Refs.|l|, 0]): the 
Zeeman splitting of spectral lines; the intensity and the polarization of syn- 
chrotron emission from free relativistic electrons; the Faraday rotation mea- 
surements (RMs) of polarized electromagnetic radiation passing through a 
ionized medium. 

Typically the Zeeman splitting, though direct, is too small to be useful 
outside our galaxy. Unfortunately, although the synchrotron emission and 
RMs allows to trace magnetic fields in very distant objects, both kind of 
measurements requires an independent determination of the local electron 
density ng. This is sometimes possible, e.g. by studying the X-ray emission 
from the electron gas when this is very hot, typically when this is confined 
in a galaxy cluster. Otherwise Ue may be not always easy to determine, 
especially for very rarefied media like the intergalactic medium (IGM). In 
the case of synchrotron emission, which intensity is proportional to n^B'^, an 
estimation of B is sometimes derived by assuming equipartition between the 
magnetic and the plasma energy densities. 

If the magnetic field to be measured is far away one relies on Faraday 
rotation. The agreement generally found between the strength of the field 
determined by RMs and that inferred from the analysis of the synchrotron 



emission in relatively close objects gives reasonable confidence on the reli- 
ability of the first method also for far away systems. It should be noted, 
however, that observations of synchrotron emission and RMs are sensitive 
to different spatial components of the magnetic field 0. The RM of the 
radiation emitted by a source with redshift Zg is given by 

RM{zs) ^ ^^ = 8.1 X 10^ r n,B\\{z){l + z)-Hl{z) ^ (1.1) 
A(A"^) JO m"^ 

where B\\ is the field strength along the line of sight and 

dl{z)= lQ-^HQ\l + z){l + Vlz)-'^dz Mpc . (1.2) 

Hq is the Hubble constant. The previous expression holds for a vanishing 
cosmological constant and modification for finite A is straightforward. This 
method requires knowledge of the electron column and possibility of field 
reversals. For nearby measurements in our own galaxy pulsar frequency and 
their decays can pin down these effects. Otherwise these stars are too far to 
help. For this reason to determine the magnetic field of the IGM by Faraday 
RMs is quite hard and only model dependent upper limits are available. 
We now briefly summarize the observational situation. 

Magnetic fields in galaxies. The interstellar magnetic field in the Milky 
Way has been determined using several methods which allowed to obtain 
valuable information about the amplitude and spatial structure of the field. 
The average field strength is 3 — 4/iG. Such a strength corresponds to an 
approximate energy equipartition between the magnetic field, the cosmic 
rays confined in the Galaxy, and the small-scale turbulent motion |[l[] 

B^ 

Pm = ^ ^ Pt^ PCR ■ (1-3) 

OTT 

Remarkably, the magnetic energy density almost coincides with energy den- 
sity of the cosmic microwave background radiation (CMBR). The field keeps 
its orientation on scales of the order of few kiloparsecs (kpc), comparable 
with the galactic size, and two reversals have been observed between the 
galactic arms, suggesting that the Galaxy field morphology may be symmet- 
rical. Magnetic fields of similar intensity have been observed in a number 
of other spiral galaxies. Although equipartition fields were observed in some 
galaxies, e.g. M33, in some others, like the Magellanic Clouds and M82, the 
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field seems to be stronger than the equipartition threshold. Concerning the 
spatial structure of the galactic fields, the observational situation is, again, 
quite confused with some galaxies presenting an axially symmetrical geom- 
etry, some other a symmetrical one, and others with no recognizable field 
structure f^. 

Magnetic fields in galaxy clusters. Observations on a large number of 
Abel clusters [^, some of which have a measured X-ray emission, give valu- 
able information on fields in clusters of galaxies. The magnetic field strength 
in the inter cluster medium (ICM) is well described by the phenomenological 
equation 

where L is the reversal field length and h^Q is the reduced Hubble constant. 
Typical values of L are 10 — 100 kpc which correspond to field amplitudes 
of 1 — 10 /iG. The concrete case of the Coma cluster [Q can be fitted with 

a core magnetic field B ~ S.S/ifoo G tangled at scales of about 1 kpc. A 
particular example of clusters with a strong field is the Hydra A cluster for 
which the RMs imply a 6 /iG field coherent over 100 kpc superimposed with a 
tangled field of strength ~ 30 /iG [H . A rich set of high resolution images of 
radio sources embedded in galaxy clusters shows evidence of strong magnetic 
fields in the cluster central regions 0. The typical central field strength 
~ 10 — 30 /iG with peak values as large as ~ 70 /iG. It is noticeable that 
for such large fields the magnetic pressure exceeds the gas pressure derived 
from X-ray data suggesting that magnetic fields may play a significant role in 
the cluster dynamics.lt is interesting, as it has been shown by Loeb and Mao 
1^ that a discrepancy exists between the estimate of the mass of the Abell 
cluster 2218 derived from gravitational lensing and that inferred from X-ray 
observations which can be well explained by the pressure support produced 
by a magnetic fields with strength ~ 50 /iG. It is still not clear if the apparent 
decrease of the magnetic field strength in the external region of clusters is 
due to the intrinsic field structure or it is a spurious effect due to the decrease 
of the gas density. Observations show also evidence for a filamentary spatial 
structure of the field. According to Eilek [^ the filaments are presumably 
structured as a flux rope, that is a twisted field structure in which the field 
lies along the axis in the center of the tube, and becomes helical going away 
from the axis. 



It seems quite plausible that all galaxy clusters are magnetized. As we 
will discuss in the next section, these observations are a serious challenge 
to most of the models proposed to explain the origin of galactic and cluster 
magnetic fields. 

Magnetic fields in high redshift objects. High resolution RMs of very 
far quasars have allowed to probe magnetic fields in the distant past. The 
most significative measurements are due to Kronberg and collaborators (see 
Ref. Ul and Refs. therein). RMs of the radio emission of the quasar 3C191, at 
z = 1.945, presumably due a magnetized shell of gas at the same redshift, are 
consistent with a field strength in the range 0.4 — 4 fiG. The field was found to 
maintain its prevailing direction over at least ~ 15 kpc, which is comparable 
with a typical galaxy size. The magnetic field of a relatively young spiral 
galaxy at z = 0.395 was determined by RMs of the radio emission of the 
quasar PKS 1229-021 lying behind the galaxy at z = 1.038. The magnetic 
field amplitude was firmly estimated to be in the range 1 — 4 /iG. Even more 
interesting was the observation of field reversals with distance roughly equal 
to the spiral arm separation, in a way quite similar to that observed in the 
Milky Way. 

Intergalactic magnetic fields. The radio emission of distant quasars is 
also used to constrain the intensity of magnetic fields in the IGM which we 
may suppose to pervade the entire Universe. As we discussed, to translate 
RMs into an estimation of the field strength is quite difficult for rarefied 
media in which ionized gas density and field coherence length are poorly 
known. Nevertheless, some interesting limits can be derived on the basis of 
well known estimates of the Universe ionization fraction and adopting some 
reasonable values of the magnetic coherence length. For example, assuming a 
cosmologically aligned magnetic field, as well as f2 = 1, A = 0, and h = 0.75, 
the RMs of distant quasar imply Bjgm ^ 10^^^ G [|l|]. A field which is aligned 
on cosmological scales is, however, unlikely. As we have seen in the above, 
in galaxy clusters the largest reversal scale is at most 1 Mpc. Adopting this 
scale as the typical cosmic magnetic field coherence length and applying the 
RM(zs) up to Za ~ 2.5, Kronberg found the less stringent limit Bjgm ^ 10~^ 
G for the magnetic field strength at present time. 

A method to determine the power spectrum of cosmic magnetic fields 
from RMs of a large number of extragalactic sources has been proposed by 
Kolatt [Q. The result of this kind of analysis would be of great help to 
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determine the origin and the time evolution of these fields. 

Another interesting idea proposed by Plaga is unfortunately not cor- 
rect. The idea here is to look to photons from an instantaneous cosmo logical 
source, like a gamma burst or a supernova, and check for the existence of a 
delayed component of the signal. This new component would be due to an 
original photon creating an electron positron pair and in turn the charged 
particle sending a photon in the original direction by inverse Compton scat- 
tering. For sources at cosmological distances the delay would be sensitive 
to a small B field, say 10^^^ G that would affect the motion of the charged 
intermediate particle. Unfortunately, the uncontrollable opening of the pair 
will produce a similar delay that cannot be disentangled from the time delay 
produced by the magnetic field. 

1.2 The alternative: dynamo or primordial ? 

For a long time the preferred mechanism to explain the aforementioned ob- 
servations was the dynamo mechanism |T^. Today, however, new observa- 



tional and theoretical results seems to point to a different scenario. Before 
trying to summarize the present status of art, a short, though incomplete, 
synthesis of what is a dynamo mechanism may be useful to some of our read- 
ers. More complete treatments of this subject can be found e.g. in Refs. 

§0,0,0 0. 

A dynamo is a mechanism responsible for the conversion of kinetic energy 
of an electrically conducting fluid into magnetic energy. It takes place when 
in the time evolution equation of the magnetic field (see e.g. Ref . |T5| ) 



fjT3 1 

^ = Vx(vxB) + — V^B. (1.5) 



where a is the electric conductivity, the first term on the RHS of Eq.( |1.5| 
(frozen-in term) dominates the second one which account for magnetic diffu- 



sion. As we will see in Sec. |1.4| this statement can be reformulated in terms of 
the magnetic Reynolds number which has to be much larger than unity. As 
it is clear from Eq. (|1.5|) , a novanishing seed field is needed to initiate the dy- 
namo process. Other three key ingredients are generally required. They are 
hydrodynamic turbulence, differential rotation and fast reconnection of mag- 
netic lines. In the frozen-in limit magnetic lines are distorted and stretched 



by turbulent motion. It can be shown |jT3| that in the same limit the ratio 
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B/ p of the magnetic field strength with the fluid density behaves like the 
distance between two fluid elements. As a consequence, a stretching of the 
field lines result in an increase of B. However, this effect alone would not 
be sufficient to explain the exponential amplification of the field generally 
predicted by the dynamo advocates. In fact, turbulence and global rota- 
tion of the fluid (e.g. by Coriolis force) may produce twisting of closed flux 
tubes and put both part of the twisted loop together, restoring the initial 
single-loop configuration but with a double flux (see Fig. 2 in Ref. ||12||). The 
process can be iterated leading to a 2"-amplification of the magnetic field 
after the n-th cycle. The merging of magnetic loops, which produce a change 



in the topology (quantified by the so called magnetic helicity, see Sec. [1.4|) 
of the magnetic field lines, requires a finite, though small, resistivity of the 
medium. This process occurs in regions of small extension where the field 



is more tangled and the diffusion time is smaller (see Sec.p..4[). As a conse 



quence, the entire magnetic configuration evolves from a small-scale tangled 
structure towards a mean ordered one. 

The most common approach to magnetic dynamo is the so called mean 
field dynamo. It is based on the assumption that fluctuations in the magnetic 
and velocity fields are much smaller of the mean slowly varying components 
of the corresponding quantities. Clearly, mean field dynamo is suitable to 
explore the amplification of large scale magnetic structures starting from 
small scale seed fields in the presence of a turbulent fluid. The temporal 
evolution of the mean component of the magnetic field is obtained by a 
suitable averaging of Eq.(|1.5|) (below, mean quantities are labelled by a 
and random quantities by a 1) 

rfVK 

^ = Vx(aBo + voxBo)-Vx [(r/ + /5) V x Bq] , (1.6) 

where 

« = -\rc{^i • V X vi) /5 = lr,(v?) , (1.7) 

7] = l/Ana is the magnetic diffusivity, and Tc is the correlation time for the 
ensemble of random velocities. The coefficient a is proportional to the helicity 
h = (vi ■ V X vi) of the flow; h measures the degree to which streamlines are 
twisted. A macroscopic parity violation is required to have a oc /i 7^ 0. One of 
the possible sources of this violation can be the Coriolis force produced by the 
rotation of the galaxy |Tl[]. The term V x (/3V x Bq) describes the additional 
field dissipation due to turbulent motion. Turbulence plays another crucial 
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role in the generation of a toroidal component of the large scale magnetic 



fields which is essential for the stability of the entire field configuration |T3 
Indeed the helicity, through the a-term, is responsible for the generation of 
an electric field parallel to Bq Q. This field provides a mode for conversion 
of toroidal into poloidal magnetic field components. This is the so called 
a-effect. To complete the "dynamo cycle" B^ ^ -Bp, another mechanism 
is required to convert the poloidal component into a toroidal one. This 
mechanism is provided by the differential rotation of the galactic disk which 
will wrap up the field line producing a toroidal field starting form a poloidal 
component; this is the cj-effect. The combination of the a and u effects gives 
rise to the, so called, a — u galactic dynamo. As a result of the coexistence of 
the poloidal and toroidal magnetic components, one of the main prediction 
of the oi a — uj dynamo is the generation of an axially symmetric mean fields. 
In the case the (3 term can be neglected, the solution of the mean field 



dynamo equation ( p..6|) can be written in the form |T^ 



Bo = {±smkz,coskz,0)e'\ (1.8) 

where z is the coordinate along the galaxy rotation axis, and 7 = —rjk'^^ak, 
fc ~ 1/L being the wavenumber. The field grows exponentially with time for 
non-zero helicity and if the scale L is sufficiently large. A general prediction 
of a dynamo mechanism is that amplification ends when equipartition is 
reached between the kinetic energy density of the small-scale turbulent fiuid 
motion and the magnetic energy density. This correspond to a magnetic 
field strength in the range of 2 — 8 fiG. Depending on the details of the 
model and of the local properties of the medium, the time required to reach 
saturation, starting from a seed magnetic fields with intensity as low as 10~^° 
G, may be of 10^ — 10^ years. It should be noted that such an estimation 
holds under the assumption that the Universe is dominated by CDM with no 
cosmological constant. If, however, as recent observations of distant type-IA 



supernovae [0 and CMB anisotropy measurments |jT^ suggest, the Universe 
posses a sizeable cosmological constant, the available time for the dynamo 
amplification increases and a smaller initial seed field may be required. This 



point has been recently rised by Davis, Lilley and Tornkvist [Q who showed 
as the required seed field might be as low as 10^^° G. 



^Our readers with some experience in field theory may recognize that by producing 
parallel electric and magnetic fields the a term is responsible of a sort of macroscopic CP 
violation. 
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In the last decade the effectiveness of the mean field dynamo has been 
questioned by several experts of the field (for a recent review see Ref. [0]). 



One of the main arguments rised by these authors against this kind of dynamo 
is that it neglects the strong amplification of small-scale magnetic fields which 
reach equipartition, stopping the process, before a coherent field may develop 
on galactic scales. 

The main, though not the unique, alternative to the galactic dynamo is 
to assume that the galactic field results directly from a primordial field which 
gets adiabatically compressed when the protogalactic cloud collapse. Indeed, 
due to the large conductivity of the intergalactic medium (see Sec.[1.4|), mag- 
netic flux is conserved in the intergalactic medium which implies that the 
magnetic field has to increase like the square of the size of the system /. It 
follows that 

/ . \2/3 

-Bprim, = -Bgal ''°'''""' • (1-9) 

V Pgal / 

Since the present time ratio between the interstellar medium density in the 
galaxies and the density of the IGM is picivi/Pgai — 10~^, and -Bgai ~ 10^^ G, 
we see that the required strength of the cosmic magnetic field at the galaxy 
formation time (z ~ 5), adiabatically rescaled to present time, is 

Spnm, 0-10-'° G. (1.10) 

This value is compatible with the observational limit on the field in the IGM 
derived by RMs, with the big-bang nucleosynthesis constraints (see Chap.|]), 
and may produce observable effects on the anisotropics of the cosmic mi- 
crowave background radiation (see Chap.^. Concerning the spatial structure 
of the galactic field produced by this mechanism, differential rotation should 
wrap the field into a symmetric spiral with field reversal along the galactic 
disk diameter and no reversal across the galactic plane 0. 

To decide between the dynamo and the primordial options astrophysicists 
have to their disposal three kind of information. They are: 

• the observations of intensity and spatial distribution of the galactic 
magnetic fields; 

• the observations of intensity and spatial distribution of the intergalactic 
magnetic fields; 

• the observations of magnetic fields in objects at high redshift. 
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Observations of the magnetic field intensity in a some galaxies, including the 
Milky Way, show evidence of approximate equipartition between turbulent 
motion and magnetic energies, which is in agreement with the prediction of 
linear dynamo. There are however some exceptions, like the M82 galaxy and 
the Magellanic Clouds, where the field strength exceed the equipartition field. 
An important test concerns the parity properties of the field with respect to 
the rotations by tt about the galactic center. As we have discussed above, 
the primordial theory predicts odd parity and the presence of reversals with 
radius (a symmetric spiral field), whereas most dynamo models predict even 
parity (axially symmetric spiral) with no reversal. Although most galaxies 
exhibit no recognizable large-scale pattern, reversals are observed between 
the arms in the Milky Way, M81 and the high redshift galaxy discussed in 
the previous section, though not in M31 and IC342. Given the low statis- 
tical significance of the sample any conclusions are, at the moment, quite 
premature [0. 

As we reported in the previous section only upper limits are available for 
the intensity of magnetic fields in the intergalactic medium. Much richer is 
the information that astrophysicists collected in the recent years about the 
magnetic fields in the inter-cluster medium (ICM). As we have seen, mag- 
netic fields of the order of 1 — 10 /iG seems to be a common features of galaxy 
clusters. The strength of these fields is comparable to that of galactic fields. 
This occurs in spite of the lower matter density of the ICM with respect to 
the density of interstellar medium in the galaxies. It seems quite difficult to 
explain the origin of the inter-cluster magnetic fields by simple ejection of 
the galactic fields. Some kind of dynamo process produced by the turbulent 
wakes behind galaxies moving in the ICM has been proposed by some au- 
thors but criticized by some others (for a review see Ref.[|T|). This problem 
has become even more critical in light of recent high-precision Faraday RMs 
which showed evidence of magnetic fields with strength exceeding 10 nG in 
the cluster central regions. According to Kronberg [|I| , the observed indepen- 
dence of the field strength from the local matter density seems to suggest that 
galactic systems have evolved in a magnetic environment where B >1 /iG. 
This hypothesis seems to be corroborated by the measurements of the Fara- 
day rotations produced by high redshift protogalactic clouds. As we wrote in 
the previous section such measurements show evidence for magnetic fields of 
the order of 1 /xG in clouds with redshift larger than 1. Since, at that time 
galaxies should have rotated few times, these observations pose a challenge 
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to the galactic dynamo advocates. We should keep in mind, however, that 
galaxy formation in the presence of magnetic fields with strength > 10~^ G 
may be problematic due to the magnetic pressure which inhibits the collapse 
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It is worthwhile to observe that primordial (or pre-galactic) magnetic 
fields are not necessarily produced in the early Universe, i.e. before recom- 
bination time. Several alternative astrophysical mechanisms have been pro- 
posed like the generation of the fields by a Biermann battery effect pO| 



see 



also Ref. ||1[). It has been suggested that the Biermann battery may produce 
seed fields which are successively amplified on galactic scale by a dynamo 



powered by the turbulence in the protogalactic cloud [jT4|, ^. This mech- 
anism, however, can hardly account for the magnetic fields observed in the 
galaxy clusters. Therefore, such a scenario would lead us to face an unnatu- 
ral situation where two different mechanisms are invoked for the generation 
of magnetic fields in galaxies and clusters, which have quite similar charac- 
teristics and presumably merge continuously at the border of the galactic 
halos. 

Another possibility is that magnetic fields may have been generated by 
batteries powered by starbursts or jet-lobe radio sources (AGNs). In a sce- 



nario recently proposed by Colgate and Li |]22[ strong cluster magnetic fields 
are produced by a dynamo operating in the accretion disk of massive black 
holes powering AGNs. We note, however, that the dynamics of the pro- 
cess leading to the formation of massive black holes is still unclear and that 
preexisting magnetic fields may be required to carry away the huge angular 
moment of the in- falling matter (see e.g. Ref. [Il9| ) . For the same reason, 



preexisting magnetic fields may also be required to trigger starbursts (see 
the end of next section). This suggests that seed fields produced before 
recombination time may anyway be required. 

In conclusion, although the data available today do not allow to answer 
yet to the question raised in this section, it seems that recent observations 
and improved theoretical work are putting in question the old wisdom in 
favour of a dynamo origin of galactic magnetic fields. Especially the recent 
observations of strong magnetic fields in galaxy clusters suggest that the 
origin of these fields may indeed be primordial. 

Furthermore, magnetic fields with strength as large as that required for 
the primordial origin of the galactic fields through gravitational compression 
of the magnetized fluid, should give rise to interesting, and perhaps necessary, 
effects for structure formation. This will be the subject of the next section. 
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1.3 Magnetic fields and structure formation 

The idea that cosmic magnetic fields may have played some role in the forma- 
tion of galaxies is not new. Some early work has been done on this subject, 



e.g. by Peblees ||2^, Rees and Rheinhardt ||2^ and especially by Wasserman 
| 25| . A considerable amount of recent papers testify the growing interest 



around this issue. A detailed review of this particular aspect of cosmology is, 
however, beyond the purposes of this report. We only summarize here few 
main points with the hope to convince the reader of relevance of this subject. 



Large scale magnetic fields modify standard equations of linear density 
perturbations in a gas of charged particles by adding the effect of the Lorentz 
force. In the presence of the field the set of Euler, continuity and Poisson 
equations become respectively ||2^ 



/5v a V ■ Vv\ 

+ _v + 

\ot a a J 

86 „ a V ■ (pv) 

— + 3-P+ ^-^ 
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Here a is the scale factor and other symbols are obvious. This set of equations 
is completed by the Faraday equation 



dja^B) 
dt 



and 



V X fv X a^B) 



V-B = 



;i.i4) 
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The term due to the Lorentz force is clearly visible in the right hand side 
of the Euler equation. It is clear that, due to this term, an inhomogeneous 
magnetic field becomes itself a source of density, velocity and gravitational 
perturbations in the electrically conducting fiuid. It has been estimated 
p5| that the magnetic field needed to produce a density contrast (5 ~ 1 



as 



required to induce structure formation on a scale /, is 



Bo{l) ~ 10 



-9 



1 Mpc 



h'Q G . 



:i.l6) 
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In his recent book Peebles Ref.[PB| pointed-out a significabt coincidence: the 
primordial magnetic field required to explain galactic fields without invoking 
dynamo amplification (see Eq. |1.10| ) would also play a relevant dynamical role 
in the galaxy formation process. 

The reader may wonder if such a dynamical role of magnetic fields is re- 
ally required. To assume that magnetic fields were the dominant dynamical 
factor at the time of galaxy formation and that they were the main source of 
initial density perturbations is perhaps too extreme and probably incompati- 
ble with recent measurements of the CMBR anisotropics. A more reasonable 
possibility is that magnetic fields are an important missing ingredient in the 
current theories on large scale structure formation (for a recent review on 
this subject see Ref.[B^). It has been argued by Coles EH] that an inho- 



mogeneous magnetic field could modulate galaxy formation in the cold dark 
matter picture (CDM) by giving the baryons a streaming velocity relative to 
the dark matter. In this way, in some places the baryons may be prevented 
from falling into the potential wells and the formation of luminous galax- 
ies on small scales may be inhibited. Such an effect could help to reconcile 
the well know discrepancy of the CDM model with clustering observations 
without invoking more exotic scenarios. 

Such a scenario received some support from a paper by Kim, Olinto and 
Rosner |^ which extended Wasserman's [^ pioneering work. Kim et al. 



determined the power spectrum of density perturbation due to a primordial 
inhomogeneous magnetic field. They showed that a present time rms mag- 
netic field of 10~^° G may have produced perturbations on galactic scale 
which should have entered the non-linear grow stage at z ~ 6, which is com- 
patible with observations. Although, Kim et al. showed that magnetic fields 
alone cannot be responsible of the observed galaxy power spectrum on large 
scales, according to the authors it seems quite plausible that in a CDM sce- 
nario magnetic fields played a not minor role by introducing a bias for the 
formation of galaxy sized objects. 

A systematic study of the effects of magnetic fields on structure forma- 
tion was recently undertaken by Battaner, Florido and Jimenez- Vicente PD 



Florido and Battaner |31|, and Battaner, Florido and Garcia- Ruiz |32]. Their 
results show that primordial magnetic fields with strength Bq < 10^^ in the 
pre-recombination era are able to produce significant anisotropic density in- 
homogeneities in the baryon-photon plasma and in the metric. In particular, 
Battaner at al. showed that magnetic fields tend to organize themselves and 
the ambient plasma into filamentary structures. This prediction seems to 



be confirmed by recent observations of magnetic fields in galaxy clusters . 
Battaner et al. suggest that such a behavior may be common to the entire 
Universe and be responsible for the very regular spider-like structure ob- 
served in the local supercluster [^] as for the filaments frequently observed 



in the large scale structure of the Universe |^. Araujo and Opher [Q] have 
considered the formation of voids by the magnetic pressure. 



An interesting hypothesis has been recently rised by Totani [^. He 
suggested that spheroidal galaxy formation occurs has a consequence of star- 
bursts triggered by strong magnetic fields. Totani argument is based on two 
main observational facts. The first is that magnetic fields strengths observed 



in spiral galaxies sharply concentrate at few microgauss (see Sec. LI), quite 
independently on the galaxy luminosity and morphology. The second point 
on which Totani based his argument, is that star formation activity has been 



observed to be correlated to the strength of local magnetic field [^. A clear 
example is given by the spiral galaxy M82, which has an abnormally large 
magnetic field of ~ 10 fiG and is known as an archetypal starburst galaxy. 
Such a correlation is theoretical motivated by the so-called magnetic braking 



||19|| : in order for a protostellar gas cloud to collapse into a star a significant 
amount of angular moment must be transported outwards. Magnetic fields 
provide a way to fulfill this requirement by allowing the presence of Alfven 
waves (see Sec. ^.2D which carry away the excess of angular moment. Whereas 
it is generally agreed that galaxy bulges and elliptical galaxies have formed 
by intense starburst activity at high redshift, the trigger mechanism lead- 
ing to this phenomenon is poorly known. According to Totani, starbursts, 
hence massive galaxy formation, take place only where the magnetic field is 
stronger of a given threshold, which would explain the apparent uniformity 
in the magnetic field amplitude in most of the observed galaxies. The value 
of the threshold field depends on the generation mechanism of the galactic 
magnetic field. Totani assumed that a seed field may have been produced 
by a battery mechanism followed by a dynamo amplification period. Such 
an assumption, however, looks not necessary and a primordial field may very 
well have produced the same final effect. 

1.4 The evolution of primordial magnetic fields 

A crucial issue for the investigation of a possible primordial origin of present 
time galactic and intergalactic magnetic fields is that concerning the time 
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evolution of the magnetic fields in the cosmic medium. Three conditions are 
needed for the persistence of large static fields: 

a) intrinsic stability of the field; 

b) the absence of free charges which could screen the field; 

c) to have a small diffusion time of the field with respect to the age of the 
Universe. 

Condition a) does not hold for strong electric fields. It is a firm prediction 
of QED that an electric field decays by converting its energy in electron- 
positron pairs if e\E\ > ml [^, ^. This, however, is a purely electric 
phenomenon. Although, at the end of the sixties, there was a claim that 
strong magnetic fields may decay through a similar phenomenon [^ the 
argument was proved to be incorrect. Only very strong fields may produce 
nontrivial instabilities in the QCD (if B > 10^^ G) and the electroweak 
vacuum (if B > 10^^ G) which may give rise to a partial screening of the 
field. These effects (see Chap.^ may have some relevance for processes which 
occurred at very early times and, perhaps, for the physics of very peculiar 
collapsed objects like magnetars ||3^. They are, however, irrelevant for the 
evolution of cosmic magnetic fields after BBN time. The same conclusion 
holds for finite temperature and densities effects which may induce screening 
of static magnetic fields (see e.g. Ref.|]4T[|). 

Condition b) is probably trivially fulfilled for magnetic fields due to the 
apparent absence of magnetic monopoles in nature. It is interesting to ob- 
serve that even a small abundance of magnetic monopoles at present time 
would have dramatic consequences for the survival of galactic and intergalac- 
tic magnetic fields which would lose energy by accelerating the monopoles. 
This argument was first used by Parker ^^ to put a severe constraint on 



the present time monopole flux, which is Fm < 10 ^^ cm ^s ^sr ^. It was 



recently proposed by Kephart and Weiler ||3^ that magnetic monopoles ac- 
celerated by galactic magnetic fields could give rise to the highest energy 
cosmic rays {E < 10^^ eV) and explain the violation of the famous Greisen- 
Zatsepin-Kuzmin cut-off [Q. 

Also the condition c) does not represent a too serious problem for the 
survival of primordial magnetic fields. The time evolution law of a magnetic 



field in a conducting medium has been already written in Eg. (1.5) 
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Neglecting fluid velocity this equation reduces to the diffusion equation 
which implies that an initial magnetic configuration will decay away in a time 

T^,s{L) = AnaL\ (1.17) 

where L is the characteristic length scale of the spatial variation of B. In 
a cosmological framework, this means that a magnetic configuration with 
coherence length Lq will survive until the present time to (t = corresponds 
to the big-bang time) only if Tdifr (i^o) > ^o- In our convention, Lq corresponds 
to the present time length scale determined by the Hubble law 

Lo = L{U) ^ , (1.18) 

a[ti) 

where a{t) is the Universe scale factor and L{ti) is the length scale at the 
time at which the magnetic configuration was formed. Note that Lq may 
not coincide with the actual size of the magnetic configuration since other 
effects (see below) may come-in to change the comoving coherence length. 
As we see from Eq.( |1.17D the relevant quantity controlling the decay time 



of a magnetic configuration is the electric conductivity of the medium. This 
quantity changes in time depending on the varying population of the available 
charge carriers and on their kinetics energies. However, since most of the 
Universe evolution takes place in a matter dominated regime, during which 
all charge carriers are non-relativistic, our estimate of the magnetic diffusion 
length is simpler. In general, electric conductivity can be determined by 
comparing Ohm law J = crE with the electric current density definition 
J = nev, where for simplicity we considered a single charge carrier type with 
charge e, number density n and velocity v. The mean drift velocity in the 
presence of the electric field E is v ~ eEr/rri where m is the charge carrier 
mass and r is the average time between particle collisions. Therefore the 
general expression for the electron conductivity is 

a = . 1.19 

m 

After recombination of electron and ions into stable atoms Universe conduc- 
tivity is dominated by residual free electrons. Their relative abundance is 



2 



In the case the average coUision time of the charge carrier is larger than the Universe 



age Tu , the latter has to be used in place of r in Eq. |1.19| [|45| 
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roughly determined by the value that this quantity took at the time when 
the rate of the reaction p + e <-^ i7 + 7 became smaller than the Universe 
expansion rate. In agreement with the results reported in Ref. ^^, we use 

ne{z) ~ 3 X 10"^° cm"^ Qoh (1 + zf , (1.20) 

where Qq is the present time density parameter and h is Hubble parameter. 
Electron resistivity is dominated by Thomson scattering off cosmic back- 
ed 
ground photons. Therefore r ~ I/ti^cxt, where ax = ^ is the Thomson 

cross section, and n^ = 4.2 x 10^(1 + z)^. Substituting these expressions in 
Eq. dLlOD we get 

2 
HP 

a = — ~10^ifio/i s-\ (1.21) 

It is noticeable that after recombination time Universe conductivity is a con- 
stant. Finally, the cosmic diffusion length, i.e. the minimal size of a magnetic 
configuration which can survive diffusion during the Universe time-life to, is 
found by substituting to = 2 x (I^q/i^)"^/^ s~^ into Eq.( |1.17D which, adopting 
r2o = 1 and h = 0.6, gives 

Ldiff ^ 2 X 10^^ cm ~ 1 A.U. (1.22) 

It follows from this result that magnetic diffusion is negligible on galactic 
and cosmological scales. 

The high conductivity of the cosmic medium has other relevant conse- 
quences for the evolution of magnetic fields. Indeed, as we already mentioned 
in the Introduction, the magnetic flux through any loop moving with fluid is 
a conserved quantity in the limit a ^ 00. More precisely, it follows from the 
diffusion equation (|1.5| ) and few vector algebra operations (see Ref. |15|)that 



d<^s{B) 



— / V X (V X B) ■ dS , (1.23) 

Ka Js 



dt Ana 

where S is any surface delimited by the loop. On scale where diffusion can be 
neglected the field is said to be frozen-in, in the sense that lines of force move 
together with the fluid. Assuming that the Universe expands isotropically Q, 



3t 



In the Sec. 2.1 we shall discuss under which hypothesis such an assumption is consistent 



with the presence of a cosmic magnetic field. 
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and no other effects come-in, magnetic flux conservation implies 
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B{t) = B{U) 4^ . (1.24) 



Tliis will be one of the most relevant equations in our review. It should be 
noted by the reader that B{t) represents the local strength of the magnetic 
field obtained by disregarding any effect that may be produced by spatial 
variations in its intensity and direction. Eq. (|1.24|) is only slightly modified 
in the case the a uniform magnetic field produce a significative anisotropic 
component in the Universe expansion (see Sec. |2.1| ). 



Another quantity which is almost conserved due to the high conductivity 
of the cosmic medium is the, so-called, magnetic helicity, defined by 

n= f d^xB-A , (1.25) 

Jv 

where A is the vector potential. Helicity is closely analogous to vorticity 
in fluid dynamics. In a field theory language, TC can be identified with the 
Chern-Simon number which is known to be related to the topological prop- 
erties of the field. Indeed, it is known from Magnetohydrodynamics (MHD) 
that 7i is proportional to the sum of the number of links and twists of the 
magnetic field lines |^^. As it follows from Eq. (|1.5|) , the time evolution of 



the magnetic helicity is determined by 

— / d^xB ■ (V X B) . (1.26) 

■Tca Jv 



dU _ 

dt Ana 



As we shall show in Chap.^, several models proposed to explain the origin of 
primordial magnetic fields predict these fields to have some relevant amount 
of helicity. 

In the previous section we have already mentioned the important role 
played by magnetic helicity in some MHD dynamo mechanisms driving an 
instability of small scale magnetic fields into large scale fields. A similar ef- 
fect may take place at a cosmological level leading to significative corrections 
to the simple scaling laws expressed by the Eqs.( LT8| , 1.24|) . Generally, these 



kind of MHD effects occur in the presence of some, turbulent motion of the 
conductive medium (note that Eqs. (|1.18|J1.24|) has been derived under the 
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assumption of vanishing velocity of the fluid v = 0). Hydrodynamic turbu- 
lence is generally parameterized in terms of the Reynolds number, defined 
by 

Re = — , (1.27) 

where u is the kinematic viscosity. Fluid motion is said to be turbulent if 
Re ^ 1. In the presence of a magnetic field another parameter turns-out to 
be quite useful. This is the magnetic Reynolds number defined by 

vL 
RcM = — , 1.28 

where t] = . When Rcm ^ 1 transport of the magnetic lines with the 

Ana 
fluid dominates over diffusion. In this case hydrodynamic turbulence in a 

conducting medium give rise to magnetic turbulence. It is often assumed 
in MHD that a fully developed magnetic turbulence give rise to equipartion 
between the kinetic and the magnetic energy of the fluid. Whether equipar- 
tition hypothesis is valid or not is a controversial issue. 

Both the hydrodynamic and magnetic Reynolds numbers can be very 
large in the early Universe. This is a consequence of the high electric con- 
ductivity and low viscosity of the medium and, especially, of the large scales 
which are involved. The electric conductivity of the early Universe has been 
computed by several authors. A first simple estimation of a in the radiation 



dominated era was performed by Turner and Widrow |^5[. In terms of the 

a 
resistivity rj their result is r/ ~ — . A more exact series of calculations can 



be found in Ref . [^ which include logarithmic corrections due to Debye and 
dynamical screening. As a result a more correct expression for rj is 

r7~|ln(l/a). (1.29) 

Other detailed computations of the Universe conductivity close to the QCD 
and the electroweak phase transitions are available in Refs. ||49|| . The kine- 



matic viscosity follows the behavior ||50 



" ^ aTHl/a) ■ ^'-^'^ 

In the early Universe u ':^ rj, i.e. ReM ^ -Re. Concerning the absolute 
value of these parameters, using the previous expressions it is easy to verify 
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that for a reasonable choice of the velocity field that may be produced by a 
phase transition, v < 10~^, both Re and Rcm are much larger than unity 
by several orders of magnitude, even for very small scales (for more details 
see Chap.^). It seems that the early Universe was a quite "turbulent child" 
! Turbulence is expected to cease after e+e~ annihilation since this process 
reduces the plasma electron population and therefore increases the photon 
diffusion length hence also the kinematic viscosity. This should happen at a 
temperature around 1 MeV. 

Turbulence is expected to produce substantial modification in the scaling 
laws of cosmological magnetic fields. This issue has been considered by sev- 



eral authors. Brandenburg, Enqvist and Olesen ||5T[ first consider MHD in 
an expanding Universe in the presence of hydro-magnetic turbulence. MHD 
equations were written in a covariant form and solved numerically under 
some simplifying assumptions. The magnetic field was assumed to be dis- 
tributed randomly either in two or three spatial dimensions. In the latter 
case a cascade (shell) model was used to reduce the number of degree of 
freedom. In both cases a transfer of magnetic energy from small to large 
magnetic configurations was observed in the simulations. In hydrodynamics 
this phenomenon is know as an inverse cascade. Cascade processes are know 
to be related to certain conservation properties that the basic equations obey 
[ |52| |. In the two-dimensional inverse cascade, the relevant conserved quantity 



is the volume integral of the vector potential squared, J d'^xA'^, whereas in 
the three-dimensional cases it is the magnetic helicity. It was recently showed 
by Son 0] that no inverse cascade can develop in 3d if the mean value of TC 
vanishes. If this is the case, i.e. in the presence of non-helical MHD turbu- 
lence, there is still an anomalous grow of the magnetic correlation length with 
respect to the scaling given in Eq.( |1.18| ) but this is just an effect of a selective 



decay mechanism: modes with larger wavenumbers decay faster than those 

whose wavenumbers are smaller. Assuming Universe expansion is negligible 

with respect to the decay time, which is given by the eddy turnover time 

tl ~ L/vl, and the decay of the large wavenumber modes does not affect 

those with smaller wavenumbers. Son found that the correlation length scale 

with time as 

/t \2/5 
L(t)~ - L, (1.31) 



where Tj = Li/vi is the eddy turnover time at t = 0. Assuming equipartition 
of the kinetic and magnetic energies, that is vl ~ B^, it follows that the 
energy decay with time like t~^^^. When the Universe expansion becomes not 
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negligible, i.e. when t > to, one has to take into account that the correlation 
length grows as r^/^, where r is the conformal time. Since r ~ T"^, it follows 
the T~^/^ law. In the real situation, the final correlation length at the present 
epoch is. 

In the above, the first factor comes from the growth of the correlation length 
in the time interval < t < to when eddy decay is faster than Universe 
expansion; the second factor comes from the growth of L in the t > to 
period; the last factor comes from trivial redshift due to the expansion of 
the Universe. Tnt is the temperature of the Universe when the fluid becomes 
non-turbulent. As we discussed Tnt ~ 1 MeV. If, for example, we assume 
that turbulence was produced at the electroweak phase transition, so that 
Ti = Tew ~ 100 GeV, that t;, ~ 0.1 and L^ ~ IQ-^rniTEw) ~ 10"^ cm, one 
find Lq ~ 100 AU. This result has to be compared with the scale one would 
have if the only mechanism of dissipation of magnetic energy is resistive 
diffusion which, as we got in Eg. ( 1.221) is ~ 1 AU. 



A larger coherence length can be obtained by accounting for the magnetic 
helicity which is probably produced during a primordial phase transition. 
The conservation of Ti has an important consequence for the evolution of 
the magnetic field. When 7i is non-vanishing, the short-scale modes are 
not simply washed out during the decay: their magnetic helicity must be 
transferred to the long-scale ones. Along with the magnetic helicity, some 
magnetic energy is also saved from turbulent decay. In other words, an inverse 
cascade is taking place. Assuming maximal helicity, i.e. that B ■ (V x B) ~ 
LB'^, the conservation of this quantity during the decay of turbulence implies 
the scaling law 



B, ~ B. y- 

This corresponds to "line averaging" , which gives a much larger amplitude of 
the magnetic field than the usual "volume averaging" . Equipartition between 
magnetic and kinetic energy implies 



L 



This relation together with the expression for the eddy decay time, tl = 
L/vl, leads to the following scaling law for the correlation length of helical 
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magnetic structures 



L 



LA I 



2/3 



:i.33) 



Comparing this result with Eq.(1.31), we see that in the hehcal case the 



correlation length grow faster than it does in the turbulent non-helical case. 
The complete expression for the scaling of L is finally obtained by including 
trivial redshift into Eq.( |1.33| ). Since in the radiation dominated era T^^ ~ 
a ~ t^/^, we have [^ 
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(1.35) 



According to Son, [^ helical hydromagnetic turbulence survives longer than 
non-helical turbulence allowing Tnt to be as low as 100 eV. If again we as- 
sume that helical magnetic turbulence is generated at the electroweak phase 
transition (which will be justified in Chap.^ we find 



L, 



Tew\'^' 



T 



nt 



To 



100 pc 



;i.36) 



which is much larger than the result obtained in the non-helical case. It is 
worthwhile to observe that, as the scale derived in the previous expression 
is also considerably larger than the cosmological magnetic diffusion length 
scale given in Eq.( p..22|) , that magnetic field fields produced by the EW phase 
transition may indeed survive until present. 



27 



Chapter 2 

Effects on the Cosmic 
Microwave Background 



2.1 The effect of a homogeneous magnetic 
field 

It is well know from General Relativity that electromagnetic fields can affect 
the geometry of the Universe. The energy momentum tensor 

T:^ = ^ {-F-'F^ + ^g-^F.^F^'^^ , (2.1) 

where F^^" is the electromagnetic field tensor, acts as a source term in the 
Einstein equations. In the case of a homogeneous magnetic field directed 
along the 2;-axis 

yOO ^ yll ^ ^22 ^ _y33 =p^ = ^ T°' = . (2.2) 

Clearly, the energy-momentum tensor becomes anisotropic due to the pres- 
ence of the magnetic field. There is a positive pressure term along the x 
and y-axes but a "negative pressure" along the field direction. It is known 
that an isotropic positive pressure increases the deceleration of the universe 
expansion while a negative pressure tends to produce an acceleration. As a 
consequence, an anisotropic pressure must give rise to an anisotropy expan- 
sion law |5^ . 



Cosmological models with a homogeneous magnetic field have been con- 



sidered by several authors (see e.g. |5J]). To discuss such models it is beyond 
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the purposes of this review. Rather, we are more interested here in the pos- 
sible signature that the pecuhar properties of the space-time in the presence 
of a cosmic magnetic field may leave on the Cosmic Microwave Background 
Radiation (CMBR). 



Following Zeldovich and Novikov [g^ we shall consider the most general 



axially sjTiimetric model with the metric 

ds^ = dt" - a\t){dx^ + dy^) - b^{t)dz^ . (2.3) 

It is convenient to define a = -; P = j', and 

a b 

r = ^ cy = a-(3. (2.4) 

Prad 

Then, assuming r, cr < 1, the Einstein equations are well approximated by 

(2.5) 
(2.6) 



d f a\ 

Jt \h) 


( a\-i -2 ^ Ar 
\h) 7t ' 7t 


dr 
dt 


2^ /.or ^ 



where H = {2a+P) and 7 are defined by the equation of state p = (7 — l)p. It 
is easy to infer from the first of the previous equations that the magnetic field 
acts so as to conserve the anisotropy that would otherwise decay with time 
in the case r = 0. By substituting the asymptotic value of the anisotropy, 
i.e. a — > 6r, into the evolution equation for r in the RD era one finds 

where g is a constant. Therefore, in the case the cosmic magnetic field is 
homogeneous, the ratio of the magnetic and blackbody radiation densities is 
not a constant, but falls logarithmically during the radiation era. 

In order to determine the temperature anisotropy of the CMBR we as- 
sume that at the recombination time tree the temperature is everywhere Tree- 
Then, at the present time, to, the temperature of relic photons coming from 
the X (or y) and z directions will be respectively 

Tx,y = Tree- = T^e CXp ( - / adt ] T^ = T^^^T = ^rec CXp ( - / (3dt 

(2.8) 
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Consequently, the expected temperature anisotropy is 

—- = -— = l-exp(/ {a-(3)dt] 

-^ -^ rec V-^trcc / 

rk) 1 /-to 

^ / {P~a)dt = — / adlnt . (2.9) 

" tree ^ " t-rcc 

By using this expression, Zeldovich and Novikov estimated that a cosmolog- 
ical magnetic field having today the strength of 10~^ -i- 10"^'^ Gauss would 
produce a temperature anisotropy 6T/T < 10~^. 

The previous analysis has been recently updated by Barrow, Ferreira and 



Silk [£5|. In that work the authors derived an upper limit on the strength 
of a homogeneous magnetic field at the recombination time on the basis 
of the 4-year Cosmic Background Explorer (COBE) microwave background 
isotropy measurements [Q. As it is well know COBE detected quadrupole 
anisotropics at a level ST/T ~ 10~^ at an angular scale of few degrees. By 
performing a suitable statistical average of the data and assuming that the 
field keeps frozen-in since the recombination till today, Barrow at al. obtained 
the limit 

Bito) < 3.5 X 10~' f/\nohl,)'/' G . (2.10) 

In the above / is a 0(1) shape factor accounting for possible non-Gaussian 
characteristics of the COBE data set. 

From this results we see that COBE data are not incompatible with a 
primordial origin of the galactic magnetic field even without invoking a dy- 
namo amplification. 



2.2 The effect on the acoustic peaks 

We will now focus our attention on possible effects of primordial magnetic 
fields on small angular scales. That is, temperature, as well polarization, 
anisotropics of the CMBR. By small angular scale (< 1°) we mean angles 
which correspond to a distance smaller than the Hubble horizon radius at 
the last scattering surface. Therefore, what we are concerning about here 
are anisotropics that are produced by causal physical mechanisms which are 
not related to the large scale structure of the space-time. 

Primordial density fiuctuations, which are necessary to explain the ob- 
served structures in the Universe, give rise to acoustic oscillations of the 
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primordial plasma when they enter the horizon some time before the last 
scattering. The oscillations distort the primordial spectrum of anisotropics 



by the following primary effects ||5^: a) they produce temperature fluctua- 
tions in the plasma, b) they induce a velocity Doppler shift of photons, c) 
they give rise to a gravitational Doppler shift of photons when they climb-out 
or fall-in the gravitational potential well produced by the density fluctuations 
(Sachs- Wolfe effect). 

In the linear regime, acoustic plasma oscillations are well described by 
standard fluid-dynamics (continuity + Euler equations) and Newtonian grav- 
ity (Poisson's equation). In the presence of a magnetic field the nature of 
plasma oscillations can be radically modified as Magneto-Hydro-Dynamics 
(MHD) has to be taken into account. 

To be pedagogical, we will first consider a single component plasma and 
neglect any dissipative effect, due for example to a finite viscosity and heat 
conductivity. We will also assume that the magnetic field is homogeneous on 
scales larger than the plasma oscillations wavelength. This choice allows us to 
treat the background magnetic field Bq as a uniform field in our equations (in 
the following symbols with the subscript stand for background quantities 
whereas the subscript 1 is used for perturbations). Within these assumptions 
the linearized equations of MHD in comoving coordinates are |J5^ Q: 



V ■ vi 
S + - - = 0, (2.11) 



where a is the scale factor. 



d 4 , V01 Bo X vi X Bo Bo X V X Bi 

a a a 47ra* ATTpoa° 



V X ( vi X Bo ) 
dtBi = ^ ^ , (2.13) 




V^0i = AirGpo U + -j^—^ ) (2.14) 

and 

V-Bi = 0, (2.15) 



^Similar equations were derived by Wasserman [g5| to the purpose to study the possible 
effect of primordial magnetic fields on galaxy formation. 
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where B = Ba^ and 6 = — , 0i and Vi are small perturbations on the 

Po 
background density, gravitational potential and velocity respectively, cs is 

the sound velocity. Neglecting its direct gravitational influence, the magnetic 

field couples to fiuid dynamics only through the last two terms in Eg. ( p. 121) . 



The first of these terms is due to the displacement current contribution to 
V X i?, whereas the latter account for the magnetic force of the current 
density. The displacement current term can be neglected provided that 

VA = , ^° < Cs , (2.16) 

^ATcip + p) 

where va is the, so called, Alfven velocity. 

Let us now discuss the basic properties of the solutions of these equations, 
ignoring for the moment the expansion of the Universe. In the absence of the 
magnetic field there are only ordinary sound waves involving density fiuctu- 
ations and longitudinal velocity fiuctuations (i.e. along the wave vector). By 
breaking the rotational invariance, the presence of a magnetic field allows 
new kind of solutions that we list below (useful references on this subject are 



59|, 13). 



1. Fast magnetosonic waves 

In the limit of small magnetic fields these waves become the ordinary 
sound waves. Their velocity, c+, is given by 

cl r^ cl + v\sm^ e , (2.17) 

where 6 is the angle between k and Bq. Fast magnetosonic waves 
involve fiuctuations in the velocity, density, magnetic field and gravi- 
tational field. The velocity and density fiuctuations are out-of-phase 
by 7r/2. Eq. (|2.17|) is valid for va << cs- For such fields the wave is 
approximatively longitudinal. 

2. Slow magnetosonic waves 

Like the fast waves, the slow waves involve both density and velocity 
fiuctuations. The velocity is however fiuctuating both longitudinally 
and transversely even for small fields. The velocity of the slow waves 
is approximatively 

ci^vlcos^e. (2.18) 
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3. Alfven waves 

For this kind of waves Bi and vi lie in a plane perpendicular to the 
plane through k and Bq. In contrast to the magnetosonic waves, the 
Alfven waves are purely rotational, thus they involve no density fluc- 
tuations. Alfven waves are linearly polarized. Their velocity of propa- 
gation is 

c'^^ = vices'^ e . (2.19) 

Detailed treatments of the evolution of MHD modes in the matter dominated 
and radiation dominated eras of the Universe can be found in Refs. pT| , |B2[ . 



5^ + -V^-^ = (2.22) 



The possible effects of MHD waves on the temperature anisotropies of the 
CMBR has been first investigated by Adams et al. [^ In the simplest case of 
magnetosonic waves, they found that the linearized equations of fluctuations 
in the Fourier space are 

6^ + Vb-3^ = 0, (2.20) 

H + V. - cie6, + ei. + ""-""^^^ - ^^^ - -i-k . (Bo X (k X B,)) = , 

(2.21) 
for the baryon component of the plasma and 

4 
3 

V^ - k^{-b^ - o^) - k^i) - aneariVb - V"^) = , (2.23) 

for the photon component. In the above \^ = ik ■ v, i? = = and 

Pj + Pi ^Pi 
Cb is the baryon sound velocity in the absence of interactions with the photon 

gas. As it is evident from the previous equations, the coupling between the 

baryon and the photons fluids is supplied by Thomson scattering with cross 

section ax- 

In the tight coupling limit (H ~ V-y) the photons provide the baryon fluid 

with a pressure term and a non-zero sound velocity. The magnetic fleld, 

through the last term in Eq.( |2.2lD , gives rise to an additional contribution 

to the effective baryon sound velocity. In the case of longitudinal waves this 

amounts to the change 

cl ^ cl + vlsin^ e . (2.24) 

In other words, the effect of the fleld can be somewhat mimicked by a vari- 
ation of the baryon density. A complication arises due to the fact that the 
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Figure 2.1: The effect of a cosmic magnetic field on the multipole moments. 
The sohd hne shows the prediction of a standard CDM cosmology (fi = l,/i = 
0.5, f^B = 0.05) with an n = 1 primordial spectrum of adiabatic fluctuations. 
The dashed line shows the effect of adding a magnetic field equivalent to 
2 X 10"'' Gauss today. From Ref.jEH 



velocity of the fast waves depends on the angle between the wave- vector and 
the magnetic field. As we mentioned previously, we are assuming that the 
magnetic field direction changes on scales larger than the scale of the fluctu- 
ation. Different patches of the sky might therefore show different fluctuation 
spectra depending on this angle. 



The authors of Ref . [^ performed an all-sky average summing also over 
the angle between the field and the line-of-sight. The effect on the CMBR 
temperature power spectrum was determined by a straightforward modifi- 
cation of the CMBFAST |Q numerical code. From the Fig. |2.1| the reader 
can see the effect of a field i?o = 2 x 10~^ G on the first acoustic peak. The 
amplitude of the peak is reduced with respect to the free field case. This 
is a consequence of the magnetic pressure which opposes the in-fall of the 
photon-baryon fluid in the potential well of the fluctuation. Although this is 
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not clearly visible from the figure, the variation of the sound velocity, hence 
of the sound horizon, should also produce a displacement of the acoustic 
peaks. The combination of these two effects may help to disentangle the 
signature of the magnetic field from other cosmological effects (for a compre- 
hensive review see [^) once more precise observations of the CMBR power 



spectrum will be available. Adams at al. derived an estimate of the sensi- 



tivity to B which MAP [^ and PLANCK [^ satellites observations should 
allow to reach by translating the predicted sensitivity of these observations 
to Qb- They found that a magnetic field with strength today i?o > 5 x 10~^ 
G should be detectable. 

It is interesting to observe that a magnetic field cannot lower the ratio of 
the first to second acoustic peak as showed by recent observations pS . 



Alfven waves may also leave a signature on the CMBR anisotropics. There 
are at least three main reasons which make this kind of wave of considerable 
interest. The first is that Alfven waves should leave a quite peculiar im- 
print on the CMBR power spectrum. In fact, as we discussed in the above, 
these waves do not involve fluctuations in the density of the photon-baryon 
fluid. Rather, they consist only of oscillations of the fluid velocity and of 
the magnetic field. Indeed, by assuming that the wavelength is smaller than 
the Hubble radius and that relativistic effects are negligible, the equations 
describing Alfven waves are 



Sb =6^ = , (2.25) 

Vfe +-v,+ ^^' --^ \ ^ ^ Bi = 0, 2.26 

a n AirpbO. 

v^ —aUeCTTi'Vb — v^) = (2.27) 

(p = . (2.28) 

Since the gravitational Doppler shift (Sachs- Wolfe effect) is absent in this 

case, the cancellation against the velocity Doppler shift which occurs for the 



acoustic modes ||5^ does not take place for the Alfven waves. This could 
provide a more clear signature of the presence of magnetic fields at the last 
scattering surface [Q. 

The second reason why Alfven waves are so interesting in this contest is 
that they are vector (or rotational) perturbations. As a consequence they 
are well suited to probe peculiar initial condition such as those that might be 
generated from primordial phase-transitions. It is remarkable that whereas 
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vector perturbations are suppressed by universe expansion and cannot arise 



from small deviations from the isotropic Friedmann Universe for t — » |l53 
this is not true in the presence of a cosmic magnetic field 0. 

The third reason of our interest for Alfven waves is that for this kind of 
waves the effect of dissipation is less serious than what it is for sound and fast 
magnetosonic waves. This issue will be touched upon in the next section. 

A detailed study of the possible effects of Alfven waves on the CMBR 
anisotropics has been independently performed by Subramanian and Barrow 
|69| and Durrer at al.[^ who reached similar results. We summarize here 



the main points of the derivation as given in Ref. |70 



In general, vector perturbations of the metric have the form 
where B and H are divergence-free, 3d vector fields supposed to vanish at 



infinity. Two gauge-invariant quantities |^ are conveniently introduced by 
the authors of Ref. 

a = H - B and n = ^r-B . (2.30) 

which represents the vector contribution to the perturbation of the extrin- 
sic curvature and the vorticity. In the absence of the magnetic field, and 
assuming a perfect fluid equation of state, the vorticity equation of motion 
is 

n + (i~3ct)-n = o . (2.31) 

a 

In the radiation dominated era the solution of this equation is fi = const. 

which clearly does not describe waves and, as we mentioned, is incompatible 

with an isotropic universe when t ^ 0. In the presence of the magnetic field, 

Durrer et al. found 

fB k")^ 
n = \°" ' , n and (2.32) 

47r(p, +p,) 

O = , '^''^ , Bi . (2.33) 



^ Collisionless matter, like e.g. gravitons after the Planck era, may however support 
nonzero vorticity even with initial conditions compatible with an isotropic universe p8]. 
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These equations describe Alfven waves propagating at the velocity va{g ■ k), 
where va is the Alfven velocity and e is the unit vector in the direction of 
the magnetic field 0. In this case some amount of initial vorticity is allowed 
which is connected to the amplitude of the magnetic field perturbation Bi 

|Oo| = {va/Bo)\Bi\ . (2.34) 

The general form of the CMBR temperature anisotropy produced by vec- 
tor perturbations is 



-1 =-V-n\-l+ a-i^dX, (2.35) 

\ J- / •'tdec 



where V = f2 — a is a gauge-invariant generalization of the velocity field. We 
see from the previous equation that besides the Doppler effect Alfven waves 
gives rise to an integrated Sachs- Wolfe term. However, since the geometric 
perturbation a is decaying with time, the integrated term is dominated by 
its lower boundary and just cancels a in V. Neglecting a possible dipole 
contribution from vector perturbations today, Durrer at al. obtained 

— {n,'k)c:in-n(k,tdec) = n-nosm(vAktdec{e-k)^ . (2.36) 



As predicted in Ref . [pg] , Alfven waves produce Doppler peaks with a period- 
icity which is determined by the Alfven velocity. Since, for reasonable values 
of the magnetic field strength, f a ^ 1 this peaks will be quite difficult to 
detect. 

Durrer et al. argued that Alfven waves may leave a phenomenologi- 
cally more interesting signature on the statistical properties of the CMBR 
anisotropies. In the absence of the magnetic field all the relevant information 
is encoded in the C/s coefficients defined by 



'ST ST , 
— (n) — (n 



{n.n'=M) ^ 



where /i = n-n'. By introducing the usual spherical harmonics decomposition 

srrji 

— (n) = ^a£„Y£^(n), (2.38) 



•^Differently form the authors of Ref. pq] , Durrer at aL assumed a homogeneous back- 
ground magnetic field . This however is not a necessary condition for the validity of the 
present considerations. 
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the C^'s are just 

Ci = {aem(^*iJ ■ (2-39) 

Because of its spin-l nature, the vorticity vector field induces transitions 
i —>■ i ±1 hence a correlation between the multipole amplitudes a^+i,m and 
ae-i,m- This new kind of correlation is encoded in the coefficients 

De{m) = (a^_i,„a|+i,„) = (a^+i,„a;_i,„) . (2.40) 

Durrer at al. |7^ determined the form of the Ci and De coefficients for the 
case of a homogeneous background magnetic field in the range —7 < n < —1, 
where n determine the vorticity power spectrum according to 

(f]o.(k)fioi(k)) = i5,,-kk)M\M) (2.41) 

Mk)= Ao -^, k<ko. (2.42) 

On the basis of these considerations they found that 4-year COBE data allow 
to obtain a limit on the magnetic field amplitude in the range —7<n< —3 
on the order of (2 — 7) x lO^^Gauss. 



2.3 Dissipative effects on the MHD modes 

In the previous section we neglected any dissipative effect which may possibly 
affect the evolution of the MHD modes. However, similar to the damping of 
baryon-photon sound waves by photon shear viscosity and heat conductivity, 
damping of MHD perturbations may also occur. This issue was studied in 



detail by Jedamzik, Katalinic and Olinto |J72| who first determined the damp- 



ing rates of fast and slow magnetosonic waves as well as of Alfven waves. 



Furthermore, it was shown in Refs. [[f^ , [73| that dissipation of MHD modes 



produce an effective damping of inhomogeneous magnetic fields. The dissi- 
pation process goes as follows. A spatially tangled magnetic field produces 
Lorentz forces which accelerate the plasma and set up oscillations. Since the 
radiation-baryon pressure is much larger than the magnetic pressure, as long 
as the photon mean-free-path is smaller than the scale of the magnetic tan- 
gle, the motions can be considered as being largely incompressible. In this 
situation mainly Alfven waves, which do not involve density fiuctuations, 
are excited. In the absence of dissipation, this process will continue until, 
for all scales A with magnetic field relaxation time r ~ X/va shorter than 
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the Hubble time t^, an approximate equipartition between magnetic and ki- 
netic energies is produced. If the fluid is non-ideal, however, shear viscosity 
will induce dissipation of kinetic energy, hence also of magnetic energy, into 
heat. In this case dissipation will end only when the magnetic field reaches 
a force-free state. 

In the absence of magnetic fields it is known that in the diffusive regime 
(i.e. when the perturbation wavelength is much larger than the mean free 
path of photon or neutrinos) acoustic density fluctuations are effectively 
damped because of the finite viscosity and heat conductivity (Silk damping 



||74|| ) . At recombination time, dissipation occurs for modes smaller than the 
approximate photon diffusion length, d^ ~ (Z^tj^)^/^, where l^ is photon mean 
free path. The dissipation of fast magnetosonic waves proceeds in a quite sim- 
ilar way. Indeed, it was showed in Ref. ||72| that the dissipation length scale 



of these kind of waves coincide with the Silk damping scale. More interesting 
is the result found in Refs. [^, ^ which shows that damping of Alfven and 
slow magnetosonic waves is significantly different from damping of sound and 
fast magnetosonic waves. The reason for such a different behavior is that, for 
a small background magnetic field f a <S 1 so that the oscillation frequency 
of an Alfven mode (vAk/a) is much smaller than the oscillation frequency 
of a fast magnetosonic mode with the same wavelength (fsound^/o)- While 
all magnetosonic modes of interest satisfy the condition for damping in the 
oscillatory regime (f sound "^ l-yk/a), an Alfven mode can become overdamped 
when the photon (or neutrino) mean-free-path becomes large enough for dis- 
sipative effects to overcome the oscillations {v a cos 9 ~ l^{T)k/a, where 9 is 
the angle between the background magnetic field and the wave vector). Be- 
cause of the strong viscosity, that prevent fluid acceleration by the magnetic 
forces, damping is quite inefficient for non-oscillating overdamped Alfven 

modes with 

2ttUT) 
A < Ao, ~ ^^ ^ . 2.43 

va cos 9 

As a result, the damping scale of overdamped Alfven modes at the end of 
the diffusion regime is smaller than the damping scale of sound and fast 
magnetosonic modes (Silk damping scale) by a factor which depends on the 
strength of the background magnetic field and the 9 angle. La ~ v a cos 9d^. 
From the previous considerations it follows that the results discussed in 
the previous section hold only under the assumption that the magnetic field 
coherence length is not much smaller than the comoving Silk damping scale 
[Ls ~ 10 Mpc), in the case of fast magnetosonic waves, and not smaller than 
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La for Alfven waves. 

Some other interesting work has been recently done by Jedamzik, Katal- 
inic and Ohnto [^ concerning the effects of dissipation of small-scale mag- 



netic fields on the CMBR. The main idea developed in the paper by Jedamzik 
et al. is that the dissipation of tangled magnetic fields before the recombi- 
nation epoch should give rise to a nonthermal injection of energy into the 
heat-bath which may distort the thermal spectrum of CMBR. It was showed 



by the authors of Ref. |^ that once photon equilibration has occurred, mainly 
via photon-electron scattering and double-Compton scattering, the resultant 
distribution should be of Bose-Einstein type with a non-vanishing chemi- 
cal potential. The evolution of the chemical potential distortions at large 
frequencies may be well approximated by |7^ 



dt tDc[z) p^ 

where, in our case, Qb = dps/dt is the dissipation rate of the magnetic 
field and tnc = 2.06 x 10^'^ s (fi{,/i^)~^z~^'^ is a characteristic time scale for 
double-Compton scattering. Jedamzik et al. assumed a statistically isotropic 
magnetic field configuration with the following power spectrum 



l^\2 = Bn^) i^if^ for k<kr, (2.45) 




and zero otherwise, normalized such that (6^) = Bq. The energy dissipation 
rate was determined by substituting this spectrum in the following Fourier 
integral 



Qb = TTTT /^'^^^ = ^TTT I d'k\h\^ (2Imcu) expf-2 j \mujdt 
87rk%- J dt 8Trk%- J \ J ^ 

(2.46) 

together with the mode frequencies for Alfven and slow magnetosonic waves 

determined in Ref.([]72|) 



-""^^—"©^^'oTiyl!)^- P'^^' 



where 3(p^ +P-y)v' = V^ ^^^ V is the shear viscosity. For k^ ^ k'jyZ^^'^, where 
k% = (15n°aT/i/2.39 x 10^9 s)^/^, an analytic solution of Eq.( p.44| ) was then 
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found to be 

r>2 / lO X (n+3) 

/^ = ^Aff^-f) • (2-48) 



In the above i^ is a numerical factor of order 1, the precise value depending 
on the spectral index n and z^ is the characteristic redshift for "freze-out" 
from double-Compton scattering. This redshift equals z^ = 2.5 x 10^ for 
typical values fi^/i^ = 0.0125, and Yp = 0.24. The scale k^zf/"^ has a simple 
interpretation. It is the scale which at redshift 2;^ is damped by one e-fold. 
For the above values of fibh"^ and Yp the corresponding comoving wavelength 
is Xn = (27r)/(A;0,z3/2) = ^g^^^ 

The present upper limit on chemical potential distortion of the CMBR 
come from the COBE/FIRAS data: |/i| < 9 x 10^^ at 95% confidence level 
|]77| . Comparing this limit with the prediction of Eq.( p.48| ) it follows that 
primordial magnetic fields of strength ^3 x 10~^ G, and comoving coherence 
length ^ 400 pc are probably excluded. On slightly larger scales, dissipation 
of spatially tangled magnetic fields may give to a different kind of CMBR 
distortion which may be described by a superposition of blackbodies of differ- 
ent temperature, i.e. a Compton y distortion [^. The absence of this kind 
of distortions in the observed CMBR thermal spectrum disallow magnetic 
fields of > 3 X 10~®Gauss on scales ~ 0.6 Mpc. 

2.4 Effects on the CMBR polarization 

Thomson scattering is a natural polarizing mechanism for the CMBR. It 
is enough if the photon distribution function seen by the electrons has a 
quadrupole anisotropy to obtain polarization. At early times, the tight cou- 
pling between the photons and the electron-baryon fluid prevents the de- 
velopment of any photon anisotropy in the baryon's rest frame, hence the 
polarization vanishes. As decoupling proceeds, the photons begin to free- 
stream and temperature quadrupole anisotropies can source a space depen- 
dent polarization. For this reason temperature and polarization anisotropies 
are expected to be correlated (for a comprehensive review on the subject see 
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The expected polarization anisotropy is not large, perhaps about 10~^. 
Currently the best polarization limit comes from the Saskatoon experiment 
PDI , with a 95% confidence level upper limit of 25 /^K at angular scales of 



about a degree, corresponding to 9 x 10 ^ of the mean temperature. Future 
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balloons and satellites observations, like e.g. the PLANCK [^ mission to be 
launched in 2007, are expected to have a good chance to measure the CMBR 
polarization power spectrum. 

Kosowsky and Loeb [^ first observed that the possible presence of mag- 



netic fields at the decoupling time may induce a sizeable Faraday rotation in 
the CMBR. Since the rotation angle depend on the wavelength, it is possi- 
ble to estimate this effect by comparing the polarization vector on a given 



direction at two different frequencies. The basic formula is [|T5 



iir'^m'^c^ 



(2.49) 



where is the amount by which the plane of polarization of linearly polarized 
radiation has been rotated, after traversing a distance L in a homogeneous 
magnetic field B in a direction q. Xg is the ionized fraction of the total 
electron density ng and m is the electron mass. Finally A is the radiation 
wavelength. 

Although the magnetic field strength is expected to be larger at early 
times, the induced Faraday rotation depends also on the free electron den- 
sity (see Eq.( |2.49| )) which drops to negligible values as recombination ends. 
Therefore, rotation is generated during the brief period of time when the 
free electron density has dropped enough to end the tight coupling but not 
so much that Faraday rotation ceases. A detailed computation requires the 



solution of the radiative transport equations in comoving coordinates [pi 



At + ifc/i(Ar + ^) = -$ - ztfAr - Ar(0) - /xH + -P2(;u)5p](2.50) 

AQtkfiAg = -k[AQ-^il-P2{fi)Sp + 2cuBAu] (2.51) 
Au + ikfiAu = -kAu - 2uobAq . (2.52) 

In the above Ay, Ag and Au respectively represent the fluctuations of tem- 
perature and of the of the Stokes parameters Q and U |jl5|. The linear 



polarization is Ap = JAq + Afj. The numerical subscripts on the radiation 
brightnesses A^ indicate moments defined by an expansion of the directional 
dependence in Legendre polynomials -P^(yu): 



Axi{k) = -J_^di2Pe{i2)Ax{k,fi). (2.53) 
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Therefore, the subscripts 0, 1, 2 label respectively monopole, dipole and 

quadrupole moments. k = XeTieCrT- is the differential optical depth and 

a 
the quantities H and R have been defined in the previous section. Derivatives 

are respect to conformal time. Finally, ub is the Faraday rotation rate 

^^ = / = o 2 2 2 - • 2.54 

From Eqs.( p.5l] , |2.52| ) wee see that Faraday rotation mixes Q and U Stokes 
parameters. The polarization brightness Ag is induced by the function Sp = 
—At (2) — Aq(2) + Aq(0) and Au is generated as Aq and Au are rotated 
into each other. In the absence of magnetic fields, Au retains its tight- 
coupling value of zero. The set of Eqs. (|2.50| - |2.52|) is not easily solved. A 



convenient approximation is the tight coupling approximation which is an 
expansion in powers of krc with tc = k^^ . This parameter measures the 
average conformal time between collisions. At decoupling the photon mean 
free path grows rapidly and the approximation breaksdown, except for long 
wavelength as measured with respect to the thickness of the last scattering 
surface. For these frequencies the approximation is still accurate. 

Kosowsky and Loeb assumed a uniform magnetic field on the scale of 
the width of the last scattering surface, a comoving scale of about 5 Mpc. 
This assumption is natural if the coherent magnetic field observed in galaxies 
comes from a primordial origin, since galaxies were assembled from a comov- 
ing scale of a few Mpc. The mean results was obtained by averaging over the 
entire sky. Therefore the equations still depend only on k and /i = cos(k ■ q) 
and not on the line-of-sight vector q and the perturbation wave-vector k 
separately. The evolution of the polarization brightnesses, for given values 
of k and /i is represented in Fig. p.2| as a function of the redshift. 

By following the approach described in the above Kosowsky and Loeb es- 
timated the polarization angle produced by a magnetic field on CMB photons 
with frequency z/i < z/2 to be 

,1/2 _, .J. "W ( Bo ^/30GHz•^2 



(,y-^i.r i--|u^^j^ J . (2.55) 



A 10 % correction may apply to this expression to account for the effects 
Ofe/i^ and ^oh^ (in the range fife/i^ > 0.007 and ^oh^ < 0.3). 

For a primordial field of Bq ~ 10^^ G which could result in the ob- 
served galactic field without dynamo amplification, one can therefore expect 
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Figure 2.2: The evolution of the polarization brightnesses, for k = 0.16 Mpc~^ 
and /i = 0.5 (in arbitrary units). Also plotted as a dotted line is the differ- 
ential visibility function fe~'^ in units of Mpc~^. From Ref. [plj . 



a rotation measure of order 1.6 deg cm ^ = 280 rad m ^. This rotation is 
considerable by astrophysical standards and could in principle be measured. 

We noticed at the beginning of this section that temperature and polar- 
ization anisotropics of the CMBR are generally expected to be correlated. 
The statistical properties of such correlation may be affected by the presence 
of a magnetic field at the decoupling time in a peculiar way. In fact, it was 



showed by Scannapieco and Ferreira |Q that such a field may induce an 
observable parity odd cross correlation between polarization and tempera- 
ture anisotropics. Any polarization pattern on the sky can be separated into 
"electric" {E) and "magnetic" {B) components. The nomenclature reflects 
the global parity property. Like multipole radiation, the harmonics of an 
ii^-mode have (—1)^ parity on the sphere, whereas those of a 5-mode have 
(—1)^+^ parity. Indeed, given a measurement of the Stokes parameters Q and 
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U, this data can be decomposed into a sum over spin ±2 spherical harmonics 

{Q±zU){n)=J2af^ ±2r,„(n) . (2.56) 

Im 

Under parity inversion, gY^m, -^ (—1)^ s^^m so that 2^^™ ± -2^fm are parity 
eigenstates. It is then convenient to define the coefficients 

E /' 2 I -2\ B '^ / 2 -2\ (c) r^\ 

SO that the E'-mode remains unchanged under parity inversion for even £, 
whereas the i?-mode changes sign. 

In an isotropic Universe, cross correlation between the B and E polar- 
izations is forbidden as this would imply parity violation. Magnetic fields, 
however, are maximally parity violating and therefore they may reveal their 
presence by producing such a cross correlation, Faraday rotation being the 



physical process which is responsible for this effect. The authors of Ref. ||8^ 
determined the expected cross correlation between temperature and E and 
B polarization modes. On the basis of such result they concluded that mag- 
netic fields strengths as low as 10~^ (present time value obtained assuming 
adiabatic scaling) could be detectable by the PLANCK satellite mission. It 
is worthwhile to note that Scannapieco and Ferreira only considered homo- 
geneous magnetic fields. We note, however, that most of their considerations 
should apply also to the case of magnetic fields with a finite coherence length. 
In this case measurements taken in different patches of the sky should present 
different temperature-polarization cross correlation depending on the mag- 
netic field and the line-of-sight direction angle. 

The consequences of Faraday rotation may go beyond the effect they pro- 
duce on the CMBR polarization. Indeed, Harari, Hayward and Zaldarriaga 



[|83| observed that Faraday rotation may also perturb the temperature power 
spectrum of CMBR. The effect mainly comes as a back-reaction of the radi- 
ation depolarization which induces a larger photon diffusion length reducing 
the viscous damping of temperature anisotropies. 

In the absence of the magnetic field {ujb = 0), to the first order in the 
tight- coupling approximation one finds 

Ac/ = 0, Aq = '^Spsin\e) (2.58) 
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and 

Sp = -^At(2) = ^tkrcAril) = -^tcAq (2.59) 

Aq = Ar(0) + $. Obviously, all multipoles with / > 3 vanish to this order. 



Replacing all quantities in terms of Aq one obtains p4 



■: ( R 16 k^Tc\ : k^ ^ k^ ,^ ,, „,,, 

(2.60) 
that can be interpreted as the equation of a forced oscillator in the presence 
of damping. 

In the presence of the magnetic field ub 7^ 0. The depolarization de- 
pends upon two angles: a) the angle between the magnetic field and wave 
propagation and b) the angle of the field with the wave vector k. Since we 
assume that the vector k is determined by stochastic Gaussian fluctuations, 
its spectrum will have no preferred direction. Therefore this dependence will 
average out when integrated. It is also assumed that for evolution purposes, 
the magnetic field has no component perpendicular to k. This imposed axial 
symmetry is compatible with the derivation of the above written Boltzmann 



equations. Under these assumptions Harari et al. found ^3 



^ ^ A A 3 Sp sin^ 9 

A„ = -Fcos9A« A« = jj^^^j^j^ (2.61) 

where the coefficient F was defined by 

Fcose = 2ujBTc (2.62) 

which gives 

Physically, F represents the average Faraday rotation between two photon- 
electron scattering. Note that assuming perfect conductivity 
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and therefore F is a time independent quantity. Faraday rotation between 
collisions becomes considerably large at frequencies around and below z/^. 
This quantity is implicitly defined by 



F . (^) (2.65) 



which gives 



.,.8.4Gh.9|^)^~27Gh. (^1 . (2.66) 



From Eqs. ( |2.61|) and, the definition of Sp given in the first part of this section, 



one can extract 

Aqo = ldo{F)Sp, Aq, = -^d2iF)Sp, (2.67) 

At, = -Sp{l - h) (2.68) 

5 

and, from the equation for At in the tight coupling, 

4 4 

Sp = — -ikrcAr, = z ^r^Ao . (2.69) 

3(3 -2d) ^ ' 3(3 -2rf) ^ " ^ ^ 

In the above the coefficients are defined so that di ^ 1 + 0{F^) for small F, 
i.e. small Faraday rotation, while di —>■ 0{1/F) as F ^ oo (for the exact 
definition see Ref.|^). Equations (|2.61| , |2.68D and (p.69|) condense the main 



effects of a magnetic field upon polarization. When there is no magnetic field 
{F = 0,d = 1) A[7 = and Ag = — ^At, sin^ 9. A magnetic field generates 
A(7, through Faraday rotation, and reduces Ag. In the limit of very large 
F (large Faraday rotation between collisions) the polarization vanishes. The 
quadrupole anisotropy At, is also reduced by the depolarizing effect of the 
magnetic field, by a factor 5/6 in the large F limit, because of the feedback 
of Aq upon the anisotropy or, in other words, because of the polarization 
dependence of Thomson scattering. The dipole At^ and monopole Atq are 
affected by the magnetic field only through its incidence upon the damp- 
ing mechanism due to photon diffusion for small wavelengths. Indeed, the 
equation for Aq = At, + $, neglecting 0{R^) contributions, now reads 

V / R 16(5-3rf) k^Tc \a ^^ a k^ r^ r. r.^.^ 

(2.70) 
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Figure 2.3: Numerical integration for the multipoles of the anisotropy cor- 
relation function in a standard CDM model without a primordial magnetic 
field (F = 0), and with -F = 1, 4, 9, which correspond to z/q = z/^, z^d/2, t'd/S 
respectively, with u^ ^ 27 GHz (E./O.OlGauss)^/^^ From Ref.|3|. 



which is the equation of a damped harmonic oscillator. 

The damping of the temperature anisotropics on small angular scales can 
be determined by solving the radiative transfer equation to second order in 
the tight-coupling approximation. By assuming solutions of the form 



Ax(t) = Axe 



(2.71) 



for X = T, Q, and U, and similarly for the baryon velocity Vb, Harari et al. 
p3| found the following solution for Eq.( |2.7CI| ) 



k 



UJ 



+ i'y 



'3(1 + R) 
where the photon-diffusion damping length-scale is 



(2.72) 



7(rf) 



k'^Tc / 8 (5 - 3d) 






6(1 + i?) U5 (3 - 2rf) 
48 



+ 



R^ 



R 



(2.73) 



The damping affects the multipole coefficients of the anisotropy power spec- 
trum which are defined by 

Ci = {47if J edkP{k)\AT,{k,To)\' . (2.74) 

The average damping factor due to photon diffusion upon the C^'s is given 
by an integral of e~^'^ times the visibihty function across the last scattering 
surface |^ ^. It depends upon cosmological parameters, notably R, and 
upon the recombination history. 

In the Fig. |2.3| it is represented the correction to the temperature power 
spectrum expected for several values of the parameter F. We see from that 
figure that on small angular scales the effect of the magnetic field is to increase 
the temperature anisotropics. The magnitude of this effect was estimated to 
be up to 7.5% in a CDM Universe on small angular scales (/ ~ 1000) at a 
level that should be reachable from future CMBR satellite experiments like 



MAP §§ and PLANCK §^. The frequency at which the effect should be 
detectable will, however, depend on the strength and coherence length of 
the magnetic field at the recombination time. Both experiment should be 
sensitive to magnetic fields around -B^=iooo = O.IG or, equivalently, Bq = 
10~^G a level that is comparable to be BEN limit (see Chap.|^). 
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Chapter 3 

Constraints from the Big Bang 
Nucleosynthesis 



The study of the influence of magnetic flelds on the Big Bang Nucleosynthesis 
(BBN) began with the pioneering works of Matese and O'Connell |^, 



and Greenstein [^ . It is remarkable that most of the more relevant effects 



were already pointed-out in those early papers. 

In their flrst paper on the subject Matese and O'Connell |^ showed that 
in the presence of very strong magnetic fields, B > B^ = eB /ml = 4.4 x 10^^ 
G (above this field strength quantized magnetic levels, "cyclotron lines", 
appear), the f3 decay rate of neutrons is significantly increased. This is mainly 
a consequence of the periodicity of the electron wave function in the plane 
normal to the field which turns into an enlarging of the electron available 
phase-space. Since the magnetic fields required to obtain any sizeable effect 
cannot be reached in the laboratory in the foreseeable future, Matese and 
O'Connell addressed their attention to the early Universe. The effects of 
primordial magnetic fields on the production of ^He during BBN was first 
considered in Ref.p^. On the basis of the results obtained in their previous 



work [^, Matese and O'Connell argued that strong magnetic fields should 
suppress the '^He relic abundance with respect to the standard case. Briefiy, 
their argument was the following. Since, after the neutron to proton ratio has 
been frozen, it takes some time for neutrons to be bounded into composite 
nuclei, a faster neutron decay due to the magnetic field implies smaller relic 
abundances of ^He and of the heavier elements. 



In Ref. [ii] two other possible effects of a magnetic field on BBN were 



shortly considered. The first of these effects consists in the variation that a 
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strong magnetic field induces on the energy density of tlie electron-positron 
gas. This effect is a consequence of the grow of the electron and positron 
phase-space in the presence of over-critical {B > Be) magnetic fields. Below 
we shall show how such an effect may have relevant consequences on the BBN 
through its action on the expansion rate of the Universe and the entropy 
transfer from the e~^e~ gas to the photons. The second effect touched by 
Matese and O'Connell concerns the influence of a uniform magnetic field 
on the Universe geometry and its consequences on the BBN Q. Matese and 
O'Connell analysis of these two effects was only qualitative and, as far as we 
know, no further work was published by these authors about these issues. 



In spite of the large number of effects considered in Ref . ||8^ Matese and 
O'Connell did not include in their analysis a simpler and quantitatively more 
relevant effect of magnetic fields on the BBN, namely the direct contribution 
of the magnetic field energy density to the expansion rate of the Universe. 
The relevance of such effect was realized by Greenstein [^ shortly after the 



publication of the Matese and O'Connell paper. Greenstein showed that by 
increasing the Universe expansion rate the presence of the magnetic field also 
increases the temperature at which the neutron-proton equilibrium ratio is 
frozen. Since this ratio is roughly given by |^| 



{n/pU=i^] exp(-g/T), (3.1) 

\mpj 

where Q = uin — rrip, a small change in the freezing temperature gives rise 
to a large variation in the neutron relative abundance hence in the relic 
abundance of the light elements. In his paper Greenstein also noted that 
if the magnetic field is sufficiently tangled over distances small compared 
to the events horizon, it will have no effect on the Universe geometry. An 
explicit calculations of the ^He relic abundance as a function of the mag- 
netic field strength were reported in a previous paper by the same author 



90| . Greenstein concluded that the effect of the magnetic field energy den- 
sity overcomes that of the magnetic field on the neutron decay discussed by 
Matese and O'Connell. Furthermore, from the requirement that the relic ^He 
mass fraction does not exceed the 28%, he inferred the upper limit B < 10^^ 
Gauss at the time when T = 5 x 10^ °K. 



In a following paper by Matese and O'Connell 88 , the authors performed 



a more careful analysis of the effects of a magnetic field on the weak reactions 



^This issue was previously considered by Thorne [Q 
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which keep neutron and protons in thermal equihbrium considering, this 
time, also the direct effect of the magnetic field on the Universe expansion 
rate. Their final conclusions were in agreement with Greenstein's result. 

The recent activity about the origin of magnetic fields during phase tran- 
sitions in the early Universe (see Chap.^ renewed the interest on the BBN 
bounds on primordial magnetic fields and induced several authors to recon- 
sider the work of Matese and O'Connell and Greenstein. It is remarkable 
that after about twenty years and a large number of new astrophysical ob- 
servations Greenstein's and Matese and O'Connell upper limit remains today 
roughly unchanged. Moreover, this is the case in spite of important develop- 
ments of the BBN numerical computations codes. 

We shall now abandon our historical approach to this section and proceed 
to give a more detailed description of the subject. 

3.1 The effect of a magnetic field on the neutron- 
proton conversion rate 

The reactions which are responsible for the chemical equilibrium of neutrons 
and protons in the early Universe are the weak processes 

(3.2) 
(3.3) 
(3.4) 

In the absence of the magnetic field and in the presence of a heat-bath, the 
rate of each of the previous processes takes the generic form 

r(12 - 34) = (ll j{^^^ (27r) V(E.P.)|-Mr /i/2(l - /3)(1 - h), (3-5) 

where Pi is the four momentum, Ei is the energy and fi is the distribution 
function of the i-th particle species involved in the equilibrium processes. 
All processes ( p.2| , |3.3| , |3^) share the same amplitude A^ determined by the 
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standard electroweak theory. 

The total neutrons to protons conversion rate is 
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(3.6) 
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where q and e are respectively the neutron-proton mass difference and the 
electron, or positron, energy, both expressed in units of the electron mass 
me- We neglect here the electron and neutrino electron chemical potentials 



as these are supposed to be small quantities during BBN [^]. The rate 1/r 
is defined by 

r" 27r3 ^^-'^ 

where G is the Fermi constant and a = Qa/ dv — —1-262. For T — > the 
integral in Eq. (p.6|) reduces to 

I = r dtt{t - qf\/t^ - 1 ~ 1.63 (3.8) 

and r„ = t jl is the neutron time-life. 

The total rate for the inverse processes (p — >■ n) can be obtained reversing 
the sign of q in Eq.(p.6|). It is assumed here that the neutrino chemical 



potential is vanishing (at the end of Sec.^]J it will also be discussed the case 
where such an assumption is relaxed). Since, at the BBN time temperature 
is much lower than the nucleon masses, neutrons and protons are assumed 
to be at rest. 



As pointed out by Matese and O'Connell p6|, ^, the main effect of a 
magnetic field stronger than the critical value B^ on the weak processes ( p.2[ - 
|3.4|) comes-in through the effect of the field on the electron, and positron, 
wave function which becomes periodic in the plane orthogonal to the field 
As a consequence, the components of the electron momentum in that 



plane are discretized and the electron energy takes the form 

Er^B) = [pI + \e\B{2n + 1 + s) + m^'] ' (3.9) 

where we assumed B to be directed along the z axis. In the above, n denotes 
the Landau level, and s = ±1 if, respectively, the electron spin is along or 
opposed to the field direction. Besides the effect on the electron dispersion 
relation, the discretization of the electron momentum due to the magnetic 
field has also a crucial effect on the phase-space volume occupied by these 
particles. Indeed, in the presence of a field with strength larger than Be the 
substitution 

J-0J3fMEo) -^ \e\Bf:{2-Sr,o) J 0^jFD{En{B),T) , (3.10) 
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has to be performed |^]. Since we only consider here magnetic fields which 
are much weaker than the proton critical value {eB <^ fn^), we can safely 
disregard any effect related to the periodicity of the proton wave function. 



The squared matrix element for each of the reactions (|3.2| - |3.4| ) is the same 
when the spin of the initial nucleon is averaged and the spins of the remaining 
particles are summed. Neglecting neutron polarization, which is very small 
for B < 10^^ G, we have |8§ 



T.\M{n)\ = l 



spins 



1-6, 



nO 



Pz 

En 



(3.11) 



It is interesting to observe the singular behaviour when a new Landau level 
opens up {En = Pz)- Such an effect is smoothed-out when temperature is 
increased P^ . 

Expressions ( p.9| ) and ( |3.10| ) can be used to determined the rate of the 
processes ( |3.2| - |3.4|) in a heat-bath and in the presence of an over-critical 
magnetic field . We start considering the neutron /?-decay. One finds 
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1 + e T, 



(3.12) 



where 7 = B/B^ and Umax is the maximum Landau level accessible to the 
final state electron determined by the requirement Pzin)"^ = q^ —m?^ — 2neB > 
0. It is noticeable that for 7 > |(5'^ — 1)^ = 2.7 only the n = term survives 
in the sum. As a consequence the /3-decay rate increases linearly with 7 above 



such a value. The computation leading to ( |3.12|) can be readily generalized 
to determine the rate of the reactions ( ^.2[ ) and ( |3.3| ) for 7 7^ 
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(3.14) 



By using the well know expression of the Euler-MacLaurin sum (see e.g. 
Ref. pip it is possible to show that in the hmit 5^-0 Eqs.( p.l^3.14| ) reduces 
to the standard expressions derived in the absence of the magnetic field 0. 

The global neutron to proton conversion rate is obtained by summing the 
last three equations 
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e T, 



(3.15) 



It is noticeable that the contribution of Eq. (|3.12|) to the total rate (|3.15|) 
is canceled by the second term of ( p. 14 ). As a consequence it follows that 
Eq. ( |3.15| ) does not depend on rimax and the n ^ p conversion grows linearly 
with the field strength above B^. From Fig.^]I]the reader can observe that, in 
the range considered for the field strength, the neutron depletion rate drops 
quickly to the free-field when the temperature grows above few MeV's. Such 
a behaviour is due to the suppression of the relative population of the lowest 
Landau level when eB 3> T^. 

In the absence of other effects, the consequence of the amplification of 
Tn^p due to the magnetic field would be to decrease the relic abundance of 
^He. In fact, a larger r„ ^p implies a lower value of the temperature [Tp) 
at which the neutron to proton equilibrium ratio is frozen because of the 
expansion of the Universe. It is evident from (|3.1| ) that the final value of 



{n/p) drops exponentially as Tp is increased. Furthermore, once n/p has 
^For a different approach see Ref.[p3[ 
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Figure 3.1: The neutron-depletion rate T^^p, normalized to the free-field 
rate, is plotted as a function of the temperature for several values of 7. From 

T3„f ITTTI 



Ref-H 



been frozen, occasional neutron /^-decays can still reduce the relic neutron 
abundance [^ . As it follows from Eq. (|3.12| ) , the presence of a strong mag- 



netic field accelerates the process which may give rise to a further suppression 
of the n/p ratio. In practice, however, neutron decay takes place at a times 
when the magnetic field strength has already decreased significantly due to 
the Universe expansion so that the effect is negligible. 

The result of Matese and O'Connell has been confirmed by Cheng et al. 

Among other effects, the authors of 



93| and by Grasso and Rubinstein 
Ref. 01 considered also QED and QCD corrections in the presence of strong 
magnetic fields . In principle these corrections may not be negligible in the 
presence of over-critical magnetic fields and their computation requires a 
proper treatment. In order to give to the reader a feeling of the relevance of 
this issue, we remind him the wrong result which was derived by O'Connell 
by neglecting QED radiative correction to the electron Dirac equation in 



the presence of a strong magnetic field. By assuming the electron anomalous 
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magnetic moment to be independent on the external field O'Connell found 



Er, 



pI + 
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(3.16) 



For B > {An /a) Be this expression give rise to negative values of the ground 
state energy which, according to O'Connell, is the manifestation of the in- 
stability of the vacuum to spontaneous production of electron-positron pairs. 
This conclusion, however, is in contradiction with standard electrodynam- 
ics from which we know that a constant magnetic field cannot transfer en- 
ergy. This problem was solved by several authors (see e.g. Ref.[^) by 
showing that by properly accounting for QED radiative corrections to the 
Dirac equation no negative value of the electron energy appear. The effect 
can be parametrized by a field dependent correction to the electron mass, 
rrie — > rrie + M, where 
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(3.17) 



Such a correction was included in Refs. ||9^, p5| . It is interesting to observe 
that although pair production cannot occur at the expense of the magnetic 
field , this phenomenon can take place in a situation of thermodynamic 
equilibrium where pair production can be viewed as a chemical reaction 
e~^ + e~ <-*> 7, the magnetic field playing the role of a catalysts agent PB| . 
We will return on this issue in Sec. |3.3| . 

Even more interesting are the corrections due to QCD. In fact. Bander 
and Rubinstein showed that in the presence of very strong magnetic fields the 
neutron-proton effective mass difference q becomes [^ (for a more detailed 
discussion of this issue see Chap.^ 



Q{B) = O.UfiNB -mn + mp + f{B) 



(3.18) 



The function f{B) gives the rate of mass change due to colour forces being 
affected by the field, ^n is the nucleon magnetic moment. For nucleons the 
main change is produced by the chiral condensate growth, which because 
of the different quark content of protons and neutrons makes the proton 
mass to grow faster [^. Although, as a matter of principle, the correction 
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to Q should be accounted in the computation of the rates that we reported 
above, in practice however, the effect on the final result is always negligible. 
More subtle it is the effect of the correction to Q on the neutron-to-proton 
equilibrium ratio. In fact, as it is evident from Eq.( p.l| ), in this case the 
correction to Q enters exponentially to determine the final neutron-to-proton 
ratio. However, the actual computation performed by Grasso and Rubinstein 



|9J] showed that the effect on the light element abundances is sub-dominant 



whenever the field strength is smaller than < 10^* Gauss. 

3.2 The effects on the expansion and cooHng 
rates of the Universe 

In the previous section we discussed how the presence of strong magnetic 
fields affects the rates of the weak reactions which are responsible for the 
chemical equilibrium of neutrons and protons before BBN. The knowledge 
of such rates is, however, not sufficient to predict the relic abundances of 
the elements synthesized during BBN. In fact, the temperature Tp at which 
{n/p)eq is frozen is determined by the competition of the weak reaction and 
the Universe expansion according to the condition 



Tn-.piTp) = H{Tp) . (3.19) 

From this expression it is clear that in order to determine Tp the knowledge 
of the Universe expansion rate H(T) is also required. 

In the absence of a cosmological term and assuming the effect of the 
magnetic field on the Universe geometry to be negligible, H is determined 
by 

H = - = i ^-^ . (3.20) 



a V 3 

where, according to the standard notation, a is the scale factor of the Uni- 
verse, G the Newton gravitational constant and p{T) is the energy density 
of the Universe. In the presence of a magnetic field 

p(r,5)=Pe,n(T,5)+p, + pB(5) (3.21) 

where Pem(T, B) is the energy density of the standard electromagnetic com- 
ponent (photons + electrons and positrons) of the heat-bath and p^ is the 
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energy density of all neutrino species (the reason why pem depends on the 
magnetic field strength will be discussed in the next section) . We see that to 
the standard components of p it adds now the contributions of the magnetic 
field energy density 

Pb{T) = ^ . (3.22) 

It is worthwhile to observe that, concerning its direct contribution to p, the 
magnetic field behaves like any relativistic component of the heat-bath. In 
fact, by assuming that the field is not too tangled on scales smaller than the 
magnetic dissipation scale, and that the Universe geometry is not affected by 
the magnetic field (see Sec. |2.1| ), the magnetic fiux conservation during the 
Universe expansion implies 

BocR-^ocT^ — ^ Pb{T)ocT\ (3.23) 

which is the same behaviour of the radiation. 



In the absence of other effects, the relation (|3.23| ) would allow to parametrize 
the effect of the magnetic field in terms of a correction to the effective number 
of massless neutrino species AN^ 19^ . Indeed, by comparing the contribu- 



tion of Nu light {m,y <^ 1 MeV) neutrino species with the energy density of 
the Universe, which is 



with ( |3.22| ) one gets 



AN^ = i^6^ (3.25) 

where b = B/T^. 

Before closing this section we have to mention another possible conse- 
quence of the faster Universe expansion induced by the presence of the mag- 
netic field. The effect is due to shortening of time between weak reactions 
freeze-out and breaking of the deuterium bottleneck. It follows that neu- 
trons have less time to decay before their confinement into nucleons take 
place which turns into a larger abundance of ^He. In Ref. [^] it was showed 
that such an effect is generally sub-dominant. 
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3.3 The effect on the electron thermodynam- 
ics 

In the above we discussed how the phase-space of electrons and positrons is 
modified by the presence of strong magnetic fields and how this effect changes 
the weak processes rates. The consequences of the variation of the electron 
phase-space, however, go well beyond that effect. Electron and positron 
thermodynamics functions will also be affected. In fact, by applying the 
prescription ( p.lO|) , we find that the number density, the energy density and 
the pressure of the electron-positron gas are now given by 

ne{B) = 7Tr^E(2-M/ fFD{T)dp, (3.26) 

eB °° r+°° 

p,{B) = — — ^(2-M/ EJ-MT)dp. (3.27) 

Pe{B) = 75^E(2-'5„„)y_^ ^^^/pD(r)<ift (3.28) 



/™(r) - i+,,E„„.| (3-29) 



where 



is the Fermi-Dirac distribution function, and En(pz) is given by ( p.9|) . As for 
the case of the weak processes rates, it is possible to show that Eqs.( |3.2(| - 
p.28 ) reduce to the their well know standard from in the limit B —>■ (see 



e.g. Ref. 

Numerical computations |^4[ show that, for small T, pe grows roughly 
linearly with B when B > Be- This effect is mainly due to - 1) the reduction, 
for each Landau level, of the area occupied by the cyclotron motion of the 
electron in plane perpendicular to the field; 2) the growth of the energy gap 
among the lowest Landau level and the n > levels, which produces an 
overpopulation of the lowest Landau level. The first effect is the dominant 
one. The net number density and the pressure of the electron-positron gas 
follow a similar behaviour. As we already noted in Sec. |3.1| , the energy cost of 
producing the electron-positron pairs excess cannot be paid by the magnetic 
field which is supposed to be static. Rather, the "power bill" is paid from 
the heat-bath, or better, from its photon component BB, 03. Especially 



in the context of BBN this point is quite relevant since the energy transfer 
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from the photons to the lowest Landau level of the electron-positron gas will 
affect the expansion rate of the Universe, its cooling rate and the effective 



baryon-to-photon ratio rj ||9^, ^, 0. We start discussing the first two effects. 
We observe that the growth of the electron and positron energy density, due 
to the presence of the magnetic field, gives rise to a faster expansion rate 
of the Universe. This point was first qualitatively discussed by Matese and 



O'Connell ^^ and recently analyzed in more detail by Grasso and Rubinstein 
||94|| . The time-temperature relation will also be modified. The relevance of 
latter effect has been first showed by Kernan, Starkman and Vachaspati ||9^ 
by solving numerically the relation 

"-t n.TTpem~r Pern fo on\ 

at dpem/dT 

where pem = Pe+P-y and pem = P-y+Pe are the energy density and the pressure 
of the electromagnetic component of the primordial heat-bath. In agreement 
with our previous considerations, Eq.( ^.3(]| ) has been obtained by imposing 
energy conservation of the electromagnetic component plasma. 

For small values of the ratio eB/T"^, the most relevant effect of the mag- 
netic field enters in the derivative dpem/dT^ that is smaller than the free field 
value. This effect can be interpreted as a delay in the electron-positron an- 
nihilation time induced by the magnetic field. This will give rise to a slower 
entropy transfer from the electron-positron pairs to the photons, then to a 
slower reheating of the heat bath. In fact, due to the enlarged phase-space 
of the lowest Landau level of electrons and positrons, the equilibrium of the 
process e~^e~ <-> 7 is shifted towards its left side. Below we will discuss as 
this effect has a clear signature on the deuterium and ^He predicted abun- 
dances. Another point of interest is that the delay in the e~^e~ annihilation 
causes a slight decrease in the T^/T ratio with respect to the canonical value, 
(4/ll)V3 [0], 

The delay in the entropy transfer from the e+e gas to the heat-bath in- 
duces also an increment in the value of baryon-to-photon ratio r/. In the ab- 
sence of other effects a larger value of 77 would induce smaller relic abundances 
of Deuterium and ^He. This effect, was first predicted in the Refs.[0, 95 



Furthermore, it is interesting to observe that in the case the primordial mag- 
netic field is inhomogeneous and it is confined in finite volume regions where 
its strength exceed the cosmic mean value (e.g. flux tubes), this effect may 
give rise to spatial variation in the relic element abundances. 
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3.4 Derivation of the constraints 

In order to account for all the effects discussed in the previous sections, 
the use of a numerical analysis is required. Usually, this is done by modi- 
fying properly the famous BBN numerical code developed by Wagoner and 
improved by Kawano ||100|| . After some discussion around the relative impor- 



tance of the different effects, the results of different groups have converged 
to a common conclusion: the most relevant effect of a cosmological mag- 
netic field on BBN is that produced by the energy density of the field on 
the Universe expansion rate. This is in qualitative agreement with the early 
result of Greenstein |5^. From a more quantitative point of view, however. 



the effect of the magnetic field on the electron thermodynamics cannot be 
totally neglected. In fact, it was showed in |^ that such an effect produces 



sizable changes in the relic abundance of ^He, Deuterium and ^He (see e.g. 
Fig.^]2| for the ^He relic abundance prediction) As a consequence, we think 



that the effect of the magnetic field on the BBN cannot be simply parame- 
terized in terms of a contribution to the effective number of neutrino species. 
Although, in this respect, a different conclusion was reached in |^8|, |101|| it 



should be noted that, differently from [Q, in those papers only approximate 
expressions for the electron thermodynamic quantities in the presence of a 
strong magnetic field were used. Such an approximation may be not justified 
when eB > T^. 

According to the standard procedure, the upper limit on the strength of 
the cosmological magnetic field was obtained in ||9^, ^ |102|] by comparing the 



numerically predicted relic abundance of '^He with the observational upper 
limit. In [^ however, the information about Deuterium and ^He was also 
used. In fact, since the effective value of rj is also affected by the magnetic 
field , 7] was chosen in the actual numerical simulation so to saturate the 
predicted value of D +^ He/H to the observational upper limit. This choice 
assured the minimal predicted abundance of ^He for each considered value 
of 5. 

By fixing N^ = 3, requiring Yp < 0.245 and D +^ He/H < 1.04 x 10"^ 
Grasso and Rubinstein derived the upper limit 

B{T = 10^ K) < 1 X 10^=^ G . (3.31) 



Similar results has been obtained by the authors of Refs. [|101| , |103 |. 



It is useful to remind to the reader under which assumptions the previous 
limit has been derived. They are the following. 
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Figure 3.2: The ^He predicted abundance is represented in function of the 
parameter 7, considered at T = 10^ °K, in three different cases: only the 
effect of the magnetic field energy density is considered (dashed line); only 
the effect of the field on the electron statistics is considered (dotted-dashed 
line); both effects are considered (continuous line). The dotted line represents 
the observational upper limit. From Ref. p5| 



1. Universe dynamics was assumed to be well described by Friedman- 
Robert son- Walker metric. In other words we are assuming that the 
magnetic field does not lead to a sizeable anisotropic component in the 
Universe expansion 0. 

2. The effective number of neutrino species is three. This means, for 
example, that if neutrinos are of Dirac type their mass and magnetic 
moment are negligible. This in order do not populate right handed 
degree of freedom by the magnetic dipole interaction of the neutrinos 
with the field. 

3. The neutrino chemical potential is negligible. 

4. Fundamental physical constants are equal to their present time values 
(for a discussion on this issue see Ref. ||104||). 



''BBN in the presence of anisotropic Universe, possibly due to a homogeneous cosmic 
magnetic field, has been considered by Thorne |54l 
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Some of these assumptions will be relaxed in the following part of this chap- 
ter. 

In order to translate our limit ( |3.31| ) into a bound on the magnetic field 
at the time of galaxy formation some caution is required. If we just assume 
that the magnetic field re-scales adiabatically with the Universe expansion, 
according to Eg. ( p. 23 ), the BBN limit reads 



Bo<7x 10-^ G . (3.32) 

We should keep in mind, however, that in this case we are neglecting any 
possible nonadiabatic evolution of the magnetic field as that which could be 
induced by a non-trivial topology of the field. Even assuming an adiabatic 
evolution, we note that the limit (|3.32|) cannot be directly interpreted as 



a limit on the progenitor of galactic magnetic fields. The reason for that 
is that BBN probes magnetic fields on scales of the order of the horizon 
radius at BBN time (the Hubble comoving radius at BBN time is ~ 100 
pc) which are much smaller than typical protogalaxy sizes (~ 1 — 10 Mpc). 
Therefore, if cosmic magnetic field are tangled on scales smaller than the 
protogalactic size, the progenitor magnetic field has to be interpreted as a 
proper average of smaller flux elements. Clearly, the result of such an average 
will depend on the statistical properties of the random magnetic field. If the 
field vector were to perform a random walk in 3d volume, the scaling would 
be B{L) = (-B)rms, L ~ A^~^/2 ||105|| , where Lq is the comoving coherence 



length of the magnetic field and A^ = L/Lq is the number of steps. An 
argument based on the statistical independence of conserved flux elements 
gives B{L) ~ A^~^ ||106|| , whereas another argument based on the statistical 



independence of the field in neighboring cells predicts B{L) ~ A^ ^/^ ||107 



Adopting a phenomenological point of view, one may just write that the rms 



field computed on the scale L at the time t is [lOS 



where p is an unknown parameter {p = 3/2, 1, 1/2 respectively in the three 
cases discussed in the above). The meaning of Bq is now understood as 
Bq = \im.L^ooB{L,to) 0. If, for example, we adopt the value p = 1 and 



4 



A detailed discussion about average procedures of tangled magnetic fields can be found 



in Ref.[109| 
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assume Lq = 100 pc, the limit ( |3.32| ) implies 

(S(lMpc,to))rms< 10-^G. (3.34) 

Therefore, although the BBN bound is much more stringent than what is 
usually claimed in the literature, it cannot exclude a primordial origin of 
galactic magnetic fields by the adiabatic compression of the field lines. 

For the same reasons which we have explained in the above, BBN limits on 
primordial magnetic fields cannot be directly compared with bounds derived 
from the analysis of CMBR anisotropies. In fact, unless the magnetic fields 
is uniform through the entire Universe, CMBR offers a probe of magnetic 
fields only on comoving scales which are much larger than the horizon radius 
at BBN time. 

We shall now consider how the previous limits changes by relaxing one of 
the assumptions under which the constraint ( p.31| ) has been derived, namely 
that related to the neutrino chemical potential. The effects of a possible 
neutrino-antineutrino asymmetry in this context has been recently consid- 



ered by Suh and Mathews ||110| . This issue is interesting since several recent 
leptogenesis scenarios predict the formation of such asymmetry during the ra- 
diation era. It is well know that even in the absence of a primordial magnetic 
field a non-vanishing neutrino chemical potential can affect the predictions 



of BBN (see Ref . \\L 1 1|| and references therein). In fact, a degeneracy of the 
electron neutrino changes both the weak reaction rates and the neutron-to- 
proton equilibrium ratio, whereas a degeneracy in any of the neutrino species 
modifies the expansion rate of the Universe. Clearly, the presence of any of 
these effects would affect the BBN limit on the strength of a primordial mag- 
netic field . Suh and Mathews found that if the limit is Bq < 5.8 x 10~^ G 
with ^e = fJ'Ue/Tuc = (in good agreement with the limit (|3.32|) ), it becomes 
Bq < 2.8 X 10^^ Gauss if ^e = jJ'vjTu^ = 0.15. Therefore, we see that in 
the presence of phenomenologically acceptable values of the neutrino chem- 
ical potential the BBN constraint on the magnetic field can be considerably 
relaxed. 
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3.5 Neutrino spin-oscillations in the presence 
of a magnetic field 

It is interesting to consider how the hmit obtained in the previous section 
changes if neutrinos carry non-standard properties which may change the 
effective neutrino number during BBN. We are especially interested here to 
the possibility that neutrinos carry non-vanishing masses and magnetic mo- 
ments. If this is the case, the dipole interaction of the neutrinos with the 
magnetic field may give rise to spin-oscillations of the neutrinos, i.e. periodic 
conversion of a helicity state into another. In the case of Dirac neutrinos, this 
phenomenon may have crucial consequences for BBN. In fact, spin-oscillation 
can populate the right-handed helicity state of the neutrino which, being 
practically sterile (for rrii, <^ T) to weak interactions, would otherwise play 
no effective role. By adding a new degree of freedom to the thermal bath, 
such an effect may produce dangerous consequences for the outcome of BBN. 
This problem was first pointed-out by Shapiro and Wasserman [|112|| and, in- 
dependently, by Lynn ||113|| who used the argument to put a constraint on 



the product of the magnetic field with the neutrino magnetic moment. In 
both works, however, the important role played by neutrino refractive prop- 
erties determined by the neutrino collective interaction with the heat-bath, 
as well as that played by neutrino scattering with leptons, were disregarded. 
A more complete treatment was developed by Fukugita et al. [|114]| . They 
showed that the conditions under which the neutrino wrong-helicity state 
can be effectively populated are the following: 

1. the spin-oscillation frequency Ai^magn = 2/ij,i? has to exceed the Uni- 
verse expansion rate; 

2. since neutrino scattering destroy the phase relationship between the 
left-handed and right-handed helicity states, Ai^magn has to be larger 
than the scattering rate; 

3. since the refractive indices for left-handed and right-handed states, ul 
and riR, are not equal, oscillations can occur only if 

Ai^magn ^ Ai^refr (3.35) 

where Ai^refr = (^l - nR)Ey. 
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The BBN is affected only if such conditions are simultaneously satisfied at 
some temperature T^^c in the range Tdoc < T^^c < Tqcd where T^ec ~ 1 MeV 
is the neutrino decoupling temperature and Tqcd ~ 200 MeV. Note that in 
the case right-handed neutrinos decouples before the QCD phase transition, 
the huge amount of entropy which is expected to be released during this 
transition would dilute their relative abundance so to prevent any effect on 



the BBN. From the previous considerations the authors of Ref. [|114| | derived 
the limit 

/^.<10-^Vb (-^j ' (3-36) 

where fiB is the Bohr magneton. The work of Fukugita et al. has been re- 
considered by several authors. For example, Enqvist, Olesen and Semikoz 
115|| , improved the previous analysis by considering the effect of the neu- 



trino refractive properties on the left-right transition probability. Elmfors, 



Grasso and Raffelt |p.l6|| accounted for the effect of the magnetic field on 
the neutrino refractive properties and used an improved treatment of neu- 
trino collisions. First of all, Elmfors et al. noted that by affecting the ther- 
modynamics properties of the electromagnetic component of the heat-bath 
(see Sec.p.3|) a strong magnetic field changes also the neutrino potentials. 
This may have relevant consequences both for neutrino spin-oscillations and 



flavour oscillations in a magnetized medium ||117|] . The interplay between 
spin-oscillations and collisions was then accounted in Ref.[ 116 | by means of 
the following evolution equation [|118 | 



^ = VxP-DPt , (3.37) 

where P is the neutrino polarization vector and Pt is its transverse compo- 
nent respect to the neutrino direction of motion. V is a vector of effective 
magnetic interaction energies which can be decomposed into its transverse 
and longitudinal components 

(3.38) 
(3.39) 
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is the left-right neutrino energy difference in the magnetized medium [[ It is 
worthwhile to note that in Eq.( |3.40|) the expression (|3.27|) has to be used for 
Pe{B). As we wrote above, collisions destroy the phase coherence between the 
left-handed and right-handed component of a neutrino state, which amounts 
to a damping of the transverse part Pt of the polarization vector. The 
main contribution to the damping rate D comes from neutrino elastic and 
inelastic scattering with leptons and equals half the total collision rate of the 
left-handed component ||121|| . In the early Universe at T ~ 1 MeV one finds: 



(/^) = /D^G^T^ 



(3.41) 



where /d is a order one numerical factor. Inserting the previous expressions 
in Eq.( |3.37| ) it is easy to derive the neutrino depolarization rate Fdepoi- In 
the small mixing angle limit. 



tan26' 






2/i5a 

AE- 



<1 



rcfr 



one finds 



depol 



(2^5t)^ (D) 



fl 



Ttt 



G^TS 



(3.42) 



(3.43) 



where /l = 1 for /i and r neutrinos, while for e neutrinos J'l ~ 3.6. By 
requiring this rate to be smaller than the Universe expansion rate H{T) in 
the temperature interval Tjec < T < Tqcd, 
upper limit Q 

'10-9 G^ 



Elmfors at al. 1161 found the 



17 



;U^ < 7 X 10 " flB 



Bn 



(3.44) 



which is not too different from the limit ( p.36 ) previously found by the au 



thors of Ref. 



Limits on fi^, were also found by the authors of Refs. 



122|| who considered the case of random magnetic fields. 



In principle, right-handed neutrinos could also be populated by direct 
spin-flip interactions mediated by virtual photons produced by scattering 
on charged particles or by annihilation processes [ 123 |, as well as by the 
interaction with small scale magnetic fields produced by thermal fluctuations 
124[. In practice, however, bounds on n,y from a possible large scale magnetic 



""For a computa tion of the neutrino refractive properties in a magnetized medium see 
also Refs.jnl, 



12C 



Note that in Ref. |116| Bq was defined as the magnetic field at BBN time. 
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field are found to be more stringent even for very weak magnetic fields. 
The most stringent upper limit on Dirac type neutrino magnetic moment 
with mass rriiy < 1 MeV, comes from stellar evolution considerations. It is 
/^v ^ 3 X 10~^^/iB [P-25| , |126|| . It is interesting to observe that if one of the 
neutrinos saturate this limit, Eq.( p.44|) implies the following quite stringent 
bound on the present time cosmic magnetic field, Bq < 10^^'^ G. 

In the particle physics standard model, neutrinos have no magnetic dipole 
moment. However, if the neutrino has a Dirac mass rriu, radiative corrections 



automatically give rise to a finite dipole moment ||113 



/x^ = 3.2 X 10"'^ /iB 



/ m^ 



VleV 



(3.45) 



even without invoking any further extension of the standard model beside 
that required to account for the finite neutrino mass. On the basis of this 
consideration Enqvist et al. ||127|| derived the following upper limit for the 



present time local magnetic field 



Bo< 



2 X lO^'^ G 






(3.46) 



eV 



Clearly, this limit cannot compete with the constraint derived in the previous 
section. 

Spin oscillations in the presence of twisted primordial magnetic fields (i.e. 



magnetic field with a nonvanishing helicity, see SecO) have been considered 
by Athar ||128|| . Athar showed that in such a situation the left-right conversion 
probabilities for neutrino and antineutrinos may be different. This result may 
open the interesting possibility that a neutrino-antineutrino asymmetry may 
be generated during the big-bang by a preexisting non trivial topology of a 



primeval magnetic field. As we shall see in Sec.O, the production of a net 
magnetic helicity of the Universe is indeed predicted by some models. 

It is also interesting to speculate on the effects when the number of di- 
mensions change, and these are large ||129|| . In fact BBN is one of the most 
serious objections to this idea, together with the background diffuse gamma 
radiation ||130|| . Detailed studies of effects of magnetic fields in these scenarios 



are not available yet. 
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Chapter 4 

Generation of magnetic fields 



4.1 Magnetic fields from primordial vorticity 

The idea that primordial magnetic fields can be produced by plasma vortical 
motion during the radiation era of the early Universe has been first pro- 
posed by Harrison ||131|| . Since this mechanism has been reviewed in several 



papers (see e.g. |1|, |132|| ) we shall not discuss it in detail here. Harrison's 



mechanism is based on the consideration that the electron and ion rotational 
velocities should decrease differently in the expanding Universe in the pre- 
recombination era. The reason is that Thomson scattering is much more 
effective for electrons than ions. Therefore electrons remains tightly coupled 
for a longer time to the radiation and behave like relativistic matter whereas 
ions are already non-relativistic. It follows that during Universe expansion 
angular velocity decreases like u oc a^^ for electrons and like uj oc a~^ for 
ions, where a is the Universe scale factor. The difference between these two 
velocities causes an electromotive force, hence an electric current which gen- 
erate a magnetic field. Harrison ||133|| showed that if a primordial turbulence 



was present at the recombination this mechanism may lead to present time 
intergalactic magnetic fields as large as 10^^ G on a scale-length of 1 Mpc. 

A problem, however, arises with this scenario. In fact, it was noted by 
Rees [|T9[ that since rotational, or vector, density perturbations decay with 
cosmic expansion, in order to produce sizeable effects at recombination time 
this kind of fluctuations should has been dominant at the radiation-matter 
equality time. This seems to be incompatible with the standard scenario for 
galaxy formation. 
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Another related problem is that, in contrast to scalar or tensor perturba- 
tions, rotational perturbations cannot arise from small deviations from the 
isotropic Friedmann Universe near the initial singularity. This is a conse- 
quence of the Helmholtz-Kelvin circulation theorem which states that the 
circulation around a closed curve following the motion of matter is con- 
served. Such a problem, however, may be partially circumvented if collision- 
less matter was present during the big-bang (e.g. decoupled gravitons after 



the Planck era). Rebhan [^ showed that in this case the Helmholtz-Kelvin 
theorem does not apply and growing modes of vorticity on superhorizon scale 
can be obtained. In fact, a non perfect fluid can support anisotropic pressure 
which may generate nonzero vorticity even if it was zero at the singular- 
ity. It follows that the only constraint to the amount of primordial vorticity 
comes from the requirement that it does not produce too large anisotropies 
in the CMBR. Rebhan showed that this requirement implies the following 
upper limit to the strength of a present time intergalactic magnetic field, 
with coherence length L, produced by vortical plasma motion 

Bo{L) < 3 X 10"^^ h-^ L^l^ G . (4.1) 

Such a field might act as a seed for galactic dynamo. 

If primordial vorticity is really not incompatible with standard cosmology, 
another interesting possibility is to generate primordial magnetic fields arises. 
It was noted by Vilenkin ||134|| that, as a consequence of parity violation in the 



Weinberg-Salam model of the electroweak interactions, macroscopic parity- 
violating currents may develop in a vortical thermal background. Vilenkin 
and Leavy ||135|| suggested that this currents may effectively give rise to 



strong magnetic fields. It was also recently noted by Brizard, Murayama and 
Wurtele ||136|| that in the presence of vorticity and a neutrino-antineutrino 



asymmetry, collective neutrino-plasma interactions may power astrophysical 
as well as cosmological magnetic fields. 

Clearly, in order to implement these scenarios, a suitable mechanism to 
produce the required amount of primordial vorticity has to be found. Among 
other exotic possibilities, generation of vorticity and magnetic fields by the 
anisotropic collapse of conventional matter into the potential well of pre- 
existing dark matter condensations in a texture-seeded scenario of structure 



formation ||132|| , and from rotating primordial black-holes ||135|| have been 
considered in the literature. 

Primordial phase transitions, may provide a more realistic source of vor- 
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ticity. The generation of primeval magnetic fields during some of these tran- 
sitions will be the subject of the next sections of this chapter. 

Another interesting possibility which is currently under study ||137 



IS 



that large scale vorticity and magnetic fields are generated at neutrino de- 
coupling in the presence of a large neutrino degeneracy. 

4.2 Magnetic fields from the quark- hadron 
phase transition 

It is a prediction of Quantum-Chromo-Dynamics (QCD) that at some very 
high temperature and/or density strongly interacting matter undergoes a de- 
confinement transition, where quark/gluon degrees of freedom are "melted- 
out" from hadrons. In the early Universe, due to the cosmological expansion, 
the process proceeds in the opposite direction starting from a "quark-gluon 
plasma" which at some critical temperature Tqqd condenses into colorless 
hadrons ||138|| . Lattice computations suggest that the QCD phase transition 



(QCDPT) is a first-order phase transition taking place at Tqcd ~ ISOMeV 
||139|| . Typically a first-order phase transition takes place by bubble nu- 
cleation. As the temperature supercools below Tqcd, sub-critical bubbles 
containing the hadronic phase grow as burning deflagration fronts releasing 
heat in the quark-gluon plasma in the form of supersonic shock fronts. When 
the shock fronts collide they reheat the plasma up to Tqcd stopping bubble 
grow. Clearly, up to this time, the transition is an out-of-equilibrium pro- 
cess. Later on, the growth of newly nucleated bubbles proceeds in thermal 
equilibrium giving rise to the, so-called, coexistence phase. The latent heat 
released by these bubbles compensate for the cooling due to the Universe ex- 
pansion keeping temperature at Tqcd- The transition ends when expansion 
wins over and the remaining quark-gluon plasma pockets are hadronized. 

The first step of the magnetogenesis scenario at the QCDPT proposed by 
Quashnock, Loeb and Spergel ||14(J| | consists in the formation of an electric 



field behind the shock fronts which precede the expanding bubbles. This 
is a consequence of the baryon asymmetry, which was presumably already 
present and which makes the baryonic components of the primordial plasma 
positively charged. At the same time, the leptonic component must be nega- 
tively charged to guarantee the Universe charge neutrality. The other crucial 
ingredient of the mechanism is the difference in the equation of state of the 
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baryonic and leptonic fluids. As a consequence, the strong pressure gradient 
produced by the passage of the shock wave gives rise to a radial electric field 
behind the shock front. Such a generation mechanism is usually known as a 
battery. Quashnock et al. gave the following estimate for the strength of the 
electric field 

^ - ^ \W/J \mo) 1^150 MevJ I £ J cm ' ^^'^^ 

where e represents the ratio of the energy densities of the two fluids, 5 = 
(iAp/p) is the pressure gradient and i is the average distance between nu- 
cleation sites. 

Small scale magnetic fields are generated by the electric currents induced 
by the electric fields. These fields, however, live on a very small scale (<C 
i )and presumably they are rapidly dissipated. Phenomenologically more 
interesting fields should be produced when the shock fronts collide giving 
rise to turbulence and vorticity on a scales of order i. Magnetic fields are 
produced on the same scale by the circulation of electric fields of magnitude 
given by Eq.( ^4.2| ). Then, using standard electrodynamics, Quashnock et 
al. found that the magnetic field produced on the scale i ~ 100 cm has a 
magnitude 

Bec^vE c^5G . (4.3) 



Following the approach which was first developed by Hogan ||105|| , the mag- 
netic field on scales L ^ i can be estimated by performing a proper volume 
average of the fields produced by a large number of magnetic dipoles of size 
i randomly oriented in space. Such an average gives 




Bl = BA-\ . (4.4) 



After the QCDPT the magnetic field evolves according to the frozen-in law 
( |1.24| ). It is straightforward to estimate the magnetic field strength at the 



recombination time on a given scale L. The smaller conceivable coherence 
length is given by the dissipation length at that time which, following the 



argument already described in Sec.pr4|, is found to be LdissiUec) — 5 x 10 
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cm, corresponding to 1 A.U. at present time. On this scale the magnetic field 
produced at the QCDPT is ~ 2 x 10~^^ G. This small strength is further 
dramatically suppressed if one considers scales of the order of the galactic 
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size ~ 10 kpc. Therefore, it looks quite unplausible that the magnetic fields 
generated by this mechanism could have any phenomenological relevance 
even if the galactic dynamo was effective. 

According to a more recent paper by Cheng and Olinto [|141|| stronger 
fields might be produced during the coexistence phase of the QCDPT. The 
new point raised by the authors of Ref . 1 141 1 is that even during such equi- 



librium phase a baryon excess build-up in front of the bubble wall, just as a 
consequence of the difference of the baryon masses in the quark and hadron 



phases. According to some numerical simulations ||142|| , this effect might en- 
hance the baryon density contrast by few orders of magnitude. Even more 
relevant is the thickness of the charged baryonic layer which, being controlled 
by baryon diffusion, is ~ 10^ fm rather than the microphysics QCD length 
scale ~ 1 fm. In this scenario magnetic fields are generated by the peculiar 
motion of the electric dipoles which arises from the convective transfer of the 
latent heat released by the expanding bubble walls. The field strength at the 
QCDPT time has been estimated by Cheng and Olinto to be as large as 

Bqcd ^ 10^ G (4.5) 

on a maximal coherence length /coh — Hqqjj. Once again, by assuming 
frozen-in evolution of the field, one can determine present time values: 

Bo ~ 10"^^ G /o - 1 PC . (4.6) 

Using Eq. ([4.4| ) Cheng and Olinto found that on the galactic length scale 
i?(kpc) ~ 10~^° G which may have same phenomenological relevance if the 
galactic dynamo is very effective. 

In a following work Sigl, Olinto and Jedamzik ||143|| investigated the pos- 



sible role that hydrodynamic instabilities produced by the expanding bubble 
walls may have in generating strong magnetic fields . Although is not clear 
whether these instabilities can really develop during the QCDPT, Sigl at al. 
claimed that this phenomenon is not implausible for a reasonable choice of 
the QCDPT parameters. By taking into account the damping due to the 



finite viscosity and heat conductivity of the plasma, the authors of Ref.[H3 
showed that the instability may grow non-linearly producing turbulence on 
a scale of the order of the bubble size at the percolation time. As a con- 
sequence, a MHD dynamo may operate to amplify seed magnetic fields and 
equipartition of the magnetic field energy with the kinetic energy may be 
reached. If this is the case, magnetic fields of the order of 10~^° G may be 
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obtained at the present time on a very large scale ~ 10 Mpc. Larger fields 
may be obtained by accounting for an inverse cascade (see Sec. |1.4|) . In the 
most optimistic situation in which the magnetic field was produced having 
maximal helicity at the QCDPT and equipartion between magnetic and ther- 
mal energy was realized at that time, Eqs.( p..34| , |1.35|) for the time evolution 
of the rms field strength and coherence length apply. By substituting in such 
equations the initial scale /(Tj) ~ ruiTQCD) ~ 30 km, and BrmrmsiTi) ~ 10^^ 
G, one finds the present time values 

S™s(To)~10-9 G Leoh~100 kpc. (4.7) 

Remarkably, we see that in this optimistic case no dynamo amplification is 
required to explain galactic, and probably also cluster, magnetic fields. 



4.3 Magnetic fields from the electroweak phase 
transition 

4.3.1 Magnetic fields from a turbulent charge flow 

Some of the ingredients which may give rise to magnetogenesis at the QGDPT 
may also be found at the electroweak phase transition (EWPT). As for the 
case of the QGDPT, magnetogenesis at the weak scale seems to require a 
first order transition. Although recent lattice computations performed in the 



framework of the standard electroweak theory ||144|| give a strong evidence 
against a first order transition, this remains a viable possibility if supersym- 
metric extension of the standard model are considered | p.45| ]. It is noticeable 
that a first order EWPT is also required for the successful realization of the 
electroweak baryogenesis scenario ||146|| . Indeed, as we shall see in the rest of 



this review, this is only one among several common aspects of baryogenesis 
and magnetogenesis. 

According to Baym, Bodeker and McLerran ||147|| strong magnetic fields 
can be generated by a first order EWPT Q via a dynamo mechanism. In 
this scenario seed fields are provided by random magnetic field fluctuations 

^At the time the paper by Baym et al. was written, a first order EWPT was thought 



to be compatible with the standard model. Therefore all computation in |147| were done 
in the framework of that model. 
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which are always present on a scale of the order of a thermal wavelength Q. 
The amplification of such seed fields proceeds as follows. When the Universe 
supercooled below the critical temperature (Tc ~ 100 GeV) the Higgs field 
locally tunneled from the unbroken SU(2) xU(l)y phase to the broken U{l)em 
phase. The tunneling gave rise to the formation of broken phase bubbles 
which then expanded by converting the false vacuum energy into kinetic 
energy. Although the bubble wall velocity is model dependent, one can find 
that for a wide range of the standard model parameters the expansion is 
subsonic (deflagration) which give rise to a supersonic shock wave ahead of 
the burning front. As the shock fronts collided turbulence should have formed 
in the cone associated with the bubble intersection. The Reynold number 
for the collision of two bubbles is 

n "Wfluid-Rbubble /. o\ 

Re ~ , (4.8) 

A 

where ffiuid ~ "Wwaii ~ 10~^ is the typical fluid velocity, -Rbubbic is the size of a 
bubble at the collision time and A is the scattering length of excitations in the 
electroweak plasma. The typical size of a bubble after the phase transition 
is completed is in the range 

-Rbubble ~ /b-f^ew (4-9) 

where 



_i mpi 



H-: ~ ^7^ ~ 10 cm (4.10) 



is the size of the event horizon at the electroweak scale, m-p\ is the Planck 
mass, g.f ~ 10^ is the number of massless degrees of freedom in the matter, 
and the fractional size /b is ~ 10~^ — 10~^. The typical scattering length A 
of excitations in the plasma is of order 

A ~ i- (4.11) 

where a„ is the fine structure constant at the electroweak scale, and (/ew ~ 9* 
is the number of degrees of freedom that scatter by electroweak processes. 



^It is worthwhile to observe here that thermal fluctuation in a dissipative plasma can 
actually produce stochastic magnetic fields on a scale larger than the thermal wavelength 



|4§ 
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By substituting these expressions Baym et al. found 

Re ~ IQ-^^a^l lna„| ~ 10^^ . (4.12) 

Such a huge Reynolds number means that turbulence fully develops at all 
scales smaller than -Rbubbie- Since conductivity is expected to be quite large 



at that time |^, magnetic fields followed the fluid motion so that a strong 
magnetic turbulence should also have been produced. In such a situation it 
is known that the kinetic energy of the turbulent flow is equipartitioned with 
the magnetic field energy. Therefore 

i?'(i?bubbie) ~ e(TcXd, (4.13) 

where e(Tc) ~ g^.T^ is the energy density of the electroweak plasma. 

In order to estimate the magnetic field strength on scale larger than 
-Rbubbie, Baym et al. treated the large scale field as a superposition of the 
field of dipoles with size -Rbubbie- This is similar to what done by other au- 
thors [|105| , |141|| (see what we wrote above for the QCDPT) apart for the fact 
that Baym et al. used a continuum approximation for the distribution of 
dipoles rather than to assume a random walk of the field lines. The density 
i/*(r) of dipoles pointing in the i-th direction was assumed to be Gaussianly 
distributed. This implies the following correlation functions for the density 
of dipoles 

V^(r)z/^(0)) = k5'^5^'^\y) (4.14) 



and for the magnetic field 

(B(r) • B(0)) ~ e^K f Sr^- ^^ • (4.15) 

J I r — Fd I I Td I 

The logarithmic divergence of the integral in these regions is cut off by the 
size of the typical dipole, f^H^, so that for r ^ f\^H~^, 

(B(r).B(0))~^ln('^V (4.16) 



By using this expression and the equipartition relation ( 4.13|) one finds that 
the strength of B^ measured by averaging on a size scale R is 
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This result can be better expressed in terms of the ratio r of (S^) to the 
energy p^ in photons which is a constant during Universe expansion in the 
absence of flux diffusion. From Eqs.( |4.17| )one gets 

(4.18) 

where Aew is the Hubble radius at the electroweak phase transition (~ 1 cm) 
times the scale factor, Tc/T^, where T^ is the photon temperature. 

From the previous results the authors of Ref. ||147|| estimated the average 




magnetic field strength at the present time. This is 

B{Um) ~ 10"^ - 10-^ G, (4.19) 

where /difr ~ 10 AU is present time diffusion length, and 

5(/gai) ~ 10-^' - 10-2° G, (4.20) 

on the galactic scale /gai ~ 10^ AU. 

4.3.2 Magnetic fields from Higgs field equilibration 

In the previous section we have seen that, concerning the generation of mag- 
netic fields, the QCDPT and the EWPT share several common aspects. How- 
ever, there is one important aspect which makes the EWPT much more inter- 
esting than the QCDPT. In fact, at the electroweak scale the electromagnetic 
field is directly influenced by the dynamics of the Higgs field which drives 
the EWPT. 

To start with we remind that, as a consequence of the Weinberg-Salam 
theory, before the EWPT is not even possible to define the electromagnetic 
field, and that this operation remains highly non-trivial until the transition is 
completed. In a sense, we can say that the electromagnetic field was "born" 
during the EWPT. The main problem in the definition of the electromagnetic 
field at the weak scale is the breaking of the translational invariance: the 
Higgs field module and its SU{2) and ?7y(l) phases take different values in 
different positions. This is either a consequence of the presence of thermal 
fluctuations, which close to T^ are locally able to break/restore the SU{2) x 
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f/y(l) symmetry or of the presence of large stable domains, or bubbles, where 
the broken symmetry has settled. 

The first generalized definition of the electromagnetic field in the presence 
of a non-trivial Higgs background was given by t'Hooft ||150|| in the seminal 
paper where he introduced magnetic monopoles in a SO (3) Georgi-Glashow 
model. t'Hooft definition is the following 

•^;:r = rG% + g-h'^'^r{D,^)\D,^r . (4.21) 

In the above G^^ = dW^ — dW^, where 

0« = ^JJl. (4.22) 

(r" are the Pauli matrices) is a unit isovector which defines the "direction" 
of the Higgs field in the S0{?>) isospace (which coincides with SU{2)) and 
{D^(f))"- = dfj,(f)°' + ge°''"^W'^(f)'^ , where Wj^ are the gauge fields components in the 
adjoint representation. The nice features of the definition (|4.21| ) are that it is 



gauge-invariant and it reduces to the standard definition of the electromag- 
netic field tensor if a gauge rotation can be performed so to have <^" = —6°'^ 
(unitary gauge). In some models, like that considered by t'Hooft, a topologi- 
cal obstruction may prevent this operation to be possible everywhere. In this 
case singular points (monopoles) or lines (strings) where 0° = appear which 
become the source of magnetic fields. t'Hooft result provides an existence 
proof of magnetic fields produced by non-trivial vacuum configurations. 

The Weinberg-Salam theory, which is based on the SU{2) x UyiX) group 
representation, does not predict topologically stable field configurations. We 
will see, however, that vacuum non-topological configurations possibly pro- 
duced during the EWPT can still be the source of magnetic fields. 

A possible generalization of the definition (|4.21|) for the Weinberg-Salam 



model was given by Vachaspati ||106|| . It is 

^""^ - -sin^vK^"^" ' "-'' ^^ 
. sin Ow 2 



F^em ^ -sinewrF^ + CO^ewFl, 



g $t$ 



where D^ = d^- i^r'^W^ - i^Y^. 



{D^<^yD,<!>-{D,<!>yD^<!> , (4.23) 



This expression was used by Vachaspati to argue that magnetic fields 
should have been produced during the EWPT. Synthetically, Vachaspati 
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argument is the following. It is known that well below the EWPT criti- 
cal temperature Tc the minimum energy state of the Universe corresponds 
to a spatially homogeneous vacuum in which $ is covariantly constant, i.e. 
D,y^ = D^cf)'^ = 0. However, during the EWPT, and immediately after 
it, thermal fluctuations give rise to a finite correlation length ^ ~ {eTc)~^. 
Therefore, there are spatial variations both in the Higgs field module |<l>| and 
in its SU{2) and f/(l)y phases which take random values in uncorrelated 



regions ^. It was noted by Davidson ||151|| that gradients in the radial part 



of the Higgs field cannot contribute to the production of magnetic fields as 
this component is electrically neutral. While this consideration is certainly 
correct, it does not imply the failure of Vachaspati argument. In fact, the 
role played by the spatial variations of the SU{2) and U{1)y "phases" of 
the the Higgs field cannot be disregarded. It is worthwhile to observe that 
gradients of these phases are not a mere gauge artifact as they correspond 
to a nonvanishing kinetic term in the Lagrangian. Of course one can always 
rotate Higgs fields phases into gauge boson degrees of freedom (see below) 
but this operation does not change Ff^ which is a gauge-invariant quantity. 

The contribution to the electromagnetic field produced by gradients of 0" 
can be readily determined by writing the Maxwell equations in the presence 
of an inhomogeneous Higgs background |152| 



ap pern 

^ fjLU 



+ 



sm 

i 



d" 



4 



flU 



9 



<|.t<|. 



{D^<^)^D,<^-D^<^{D,<^y 



(4.24) 



Even neglecting the second term on the righthand side of Eq.( |4.24j ), which 
depends on the definition of F^™ in a Higgs inhomogeneous background (see 
below), it is evident that a nonvanishing contribution to the electric 4-current 
arises from the covariant derivative of 0". The physical meaning of this 
contribution may look more clear to the reader if we write Eq.( 4.24] ) in the 
unitary gauge 






— le 



+ie \W^^ {D,W^) - W^ [D,W, 



a 



^t 



(4.25) 



^ Vachaspati |l06| did also consider Higgs field gradients produced by the presence of 
the cosmological horizon. However, since the Hubble radius at the EWPT is of the order 
of 1 cm whereas ^ ^ (eTc)^^ ~ 10^^^ cm, it is easy to realize that magnetic fields possibly 
produced by the presence of the cosmological horizon are phenomenologically irrelevant. 
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led^' (WIW, - W^Wl) 



Not surprisingly, we see that the electric currents produced by Higgs field 
equilibration after the EWPT are nothing but W boson currents. 

Since, on dimensional grounds, -Di/$ ~ v/S, where v is the Higgs field 
vacuum expectation value, Vachaspati concluded that magnetic fields (elec- 
tric fields were supposed to be screened by the plasma) should have been 
produced at the EWPT with strength 

B ~ sinOwgT'^ ^ 10^^ G . (4.26) 

Of course these fields live on a very small scale of the order of ^ and in 
order to determine fields on a larger scale Vachaspati claimed that a suitable 
average has to be performed (see return on this issue below in this section). 
Before discussing averages, however, let us try to understand better the 
nature of the magnetic fields which may have been produced by the Vachas- 
pati mechanism. We notice that Vachaspati's derivation does not seem to 
invoke any out-of-equilibrium process and indeed the reader may wonder 
what is the role played by the phase transition in the magnetogenesis. More- 
over, magnetic fields are produced anyway on a scale (eT)~^ by thermal 
fiuctuations of the gauge fields so that it is unclear what is the difference 
between magnetic fields produced by the Higgs fields equilibration and these 
more conventional fields. In our opinion, although Vachaspati's argument is 
basically correct its formulation was probably oversimplified. Indeed, several 
works showed that in order to reach a complete understanding of this phys- 
ical effect a more careful study of the dynamics of the phase transition is 
called for. We shall now review these works starting from the case of a first 
order phase transition. 

The case of a first order EWPT 

Before discussing the SU{2) x f/(l) case we cannot overlook some im- 
portant work which was previously done about phase equilibration during 
bubble collision in the framework of more simple models. In the context of 



a f/(l) Abelian gauge symmetry. Kibble and Vilenkin |p.53|| showed that the 
process of phase equilibration during bubble collisions give rise to relevant 
physical effects. The main tool developed by Kibble and Vilenkin to investi- 
gate this kind of processes is the, so-called, gauge-invariant phase difference 



defined by 

Ae= [ dx'' DkO , (4.27) 

where 6 is the U{1) Higgs field phase and D^O = d^6 + eA^ is the phase 
covariant derivative. A and B are points taken in the bubble interiors and 
/c = 1, 2, 3. A^ obeys the Klein-Gordon equation 

((9^ + m^)M = (4.28) 

where m = ev is the gauge boson mass. Kibble and Vilenkin assumed that 
during the collision the radial mode of the Higgs field is strongly damped so 
that it rapidly settles to its expectation value v everywhere. One can choose 
a frame of reference in which the bubbles are nucleated simultaneously with 
centers at {t,x,y,z) = (0, 0, 0, ±_Rc)- In this frame, the bubbles have equal 
initial radius Ri = Rq. Their first collision occurs at (tc, 0,0,0) when their 



radii are Re and t^ = JR^ — Rq- Given the symmetry of the problem about 
the axis joining the nucleation centers (2;-axis), the most natural gauge is the 
axial gauge. In this gauge 

e{x) = e{T, z), A"(x) = x'^air, z) , (4.29) 

where a = 0, 1, 2 and r^ = t^ — x^ — y'^. The condition Oair, 0) = fixes the 
gauge completely. At the point of first contact z = 0, r = tc the Higgs field 
phase was assumed to change from 9q to —9q going from a bubble into the 
other. This constitutes the initial condition of the problem. The following 
evolution of 9 is determined by the Maxwell equation 

d''F^,=j^ = -ev^ D^9 (4.30) 

and the Klein-Gordon equation which splits into 

dl9a + -0^9 - 819 + m^9 = , (4.31) 

r 

5^a + -dra - dla + m^a = . (4.32) 

r 



The solution of the linearized equations ( [4.31| ) and ( [4. 321 ) for r > tc then 
becomes 

n / X ^o^c f°° dk , / , ,1 , n\ / N 

9a{r, z) = / -—smkz [cosuj{t — tc) -\ smti;(r — tc) (,4-33) 

vrr J-oo k V ujtc / 
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OomHc r°° dk 

sin kz 



^ J-oo k 



+ ( \ T- ) smu;(r - t 

\uj uj'H 



T — tc , , 

— t--cosuj[t -tc) 

(4.34) 



where u;^ = /c^ + irP' . The gauge-invariant phase difference is deduced by the 
asymptotic behavior at z ^ ±00 

A^ = ^„(t, 0, 0, +00) -^„(t,0, 0,-00) (4.35) 



t 



cos m{t — tc) H sin m(t — tc) ) 



Thus, phase equihbration occurs with a time scale tc determined by the bub- 
ble size, with superimposed oscillations with frequency given by the gauge- 
field mass. As we see from Eq.( 4.34|) phase oscillations come together with 



oscillations of the gauge field. It follows from Eq. (|4.3CI|) that these oscillations 
give rise to an "electric" current. This current will source an "electromag- 
netic" field strength F^^, 0. Because of the symmetry of the problem the only 
nonvanishing component of F^i, is 

paz ^ x^'d.aiT, z) . (4.36) 

Therefore, we have an azimuthal magnetic field B''' = F^'' = pdza and a 
longitudinal electric field E^ = F°^ = —tdza = —{t/p)B'^{T,z), where we 
have used cylindrical coordinates (p, 0). We see that phase equilibration 
during bubble collision indeed produces some real physical effects. 

Kibble and Vilenkin did also consider the role of electric dissipation. They 
showed that a finite value of the electric conductivity a gives rise to a damping 
in the "electric" current which turns into a damping for the phase equilibra- 
tion. They found 

-at/2 ( ■ ^ 



Ae(t) = 2606-"*' ' cos mt + — sin mt (4.37) 

2m 



for small values of cr, and 



M{t) = 2^0 exp (-mH/a) (4.38) 



'*It is understood that since the toy model considered by Kibble e Vilenkin is not 
SU{2) X [/(l)y, Ffj^i, is not the physical electromagnetic field strength. 
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in the opposite case. The dissipation time scale is typically much smaller 
than the time which is required for two colliding bubble to merge completely. 
Therefore the gauge-invariant phase difference settles rapidly to zero in the 
overlapping region of the two bubbles and in its neighborhood. It is inter- 
esting to compute the line integral of DkO over the path ABCD represented 
in the Fig.[4.1|. From the previous considerations it follows that A6a_b = 0, 




Figure 4.1: Two colliding bubbles are depicted. The gauge invariant phase 



difference is computed along the path ABCD. (From Ref. ||153||). 



AG AD = AObc = and A9dc = 26*0. It is understood that in order for the 
integral to be meaningful, the vacuum expectation value of the Higgs field 
has to remain nonzero in the collision region and around it, so that the phase 
6 remains well defined and interpolates smoothly between its values inside 
the bubbles. Under these hypothesis we have 



D.edx' 



2ft 



ABCD 



• 



(4.39) 



The physical meaning of this quantity is recognizable at a glance in the 
unitary gauge, in which each A^ is given by a line integral of the vector 
potential A. We see that the gauge-invariant phase difference computed 
along the loop is nothing but the magnetic flux trough the loop itself 



0(5) = / AkOdx'' = - f DkOdx'' 

J ABCD e J ABCD 



290 
e 



(4.40) 
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In other words, phase equihbration give rise to a ring of magnetic flux near 
the circle on which bubble walls intersect. If the initial phase difference 
between the two bubbles is 2it, the total flux trapped in the ring is exactly 
one flux quantum, 27r/e. 

Kibble and Vilenkin did also consider the case in which three bubbles 
collide. They argued that in this case the formation of a string, in which 
interior symmetry is restored, is possible. Whether or not this happens is 
determined by the net phase variation along a closed path going through the 
three bubbles. The string forms if this quantity is larger than 27r. According 
to Kibble and Vilenkin strings cannot be produced by two bubble collisions 
because, for energetic reasons, the system will tend to choose the shorter of 
the two paths between the bubble phases so that a phase displacement > 27r 
can never be obtained. This argument, which was first used by Kibble [[154| 
for the study of defect formation, is often called the "geodesic rule" . 

The work of Kibble and Vilenkin was reconsidered by Copeland and Safiin 
155 1 and more recently by Copeland, Safiin and Tornkvist [|156|| who showed 



that during bubble collision the dynamics of the radial mode of the Higgs 
field cannot really be disregarded. In fact, violent fluctuations in the modulus 
of the Higgs field take place and cause symmetries to be restored locally, al- 
lowing the phase to "slip" by an integer multiple of 27r violating the geodesic 
rule. Therefore strings, which carry a magnetic flux, can be produced also 
by the collision of only two bubbles. Safiin and Copeland ||157|| went a step 
further by considering phase equilibration in the SU{2) x f/(l) case, namely 
the electroweak case. They showed that for some particular initial conditions 
the SU{2) X U{1) Lagrangian is equivalent to a U{1) Lagrangian so that part 



of Kibble and Vilenkin ||153|| considerations can be applied. The violation of 
the geodesic rule allows the formation of vortex configurations of the gauge 
fields. Safiin and Copeland argued that these configurations are related to 



the Nielsen-Olesen vortices [[158|| . Indeed, it is know that such a kind of non- 



perturbative solutions are allowed by the Weinberg-Salam model ||159|| (for a 



comprehensive review on electroweak strings see Ref. ||160|| ). Although elec 



troweak string are not topologically stable, numerical simulations performed 
in Ref. | |157|| show that in presence of small perturbations the vortices survives 



on times comparable to the time required for bubble to merge completely. 

The generation of magnetic fields in the SU{2) x f/(l)y case was not 
considered in the work by Safiin and Copeland. This issue was the subject 
of a following paper by Grasso and Riotto [ |152|| . The authors of Ref. [ [152| 
studied the dynamics of the gauge fields starting from the following initial 
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Higgs field configuration 

$i„(x) = i= ( .° ) + _ exp I -^^nV" 1 I , " , ,^„ 1 (4.41) 




which represents the superposition of the Higgs fields of two bubbles which 
are separated by a distance b. In the above n" is a unit vector in the SU{2) 
isospace and r" are the Pauli matrices. The phases and the orientation of 
the Higgs field were chosen to be uniform across any single bubble. It was 
assumed that Eq.( [4.41| ) holds until the two bubble collide (t = 0). Since n°'T°- 



is the only Lie-algebra direction which is involved before the collision, one 
can write the initial Higgs field configuration in the form | |157| ] 



x/2^^^ V ' ^ J V pi'^y 



$i.(x) = ^exp ( -^^n^r'^ ] ( ^,^,^",^(,) (4.42) 



In order to disentangle the peculiar role played by the Higgs field phases, the 
initial gauge fields W^^ and their derivatives were assumed to be zero at t = 0. 
This condition is of course gauge dependent and should be interpreted as a 
gauge choice. It is convenient to write the equation of motion for the gauge 
fields in the adjoint representation. For the SU{2) gauge fields we have 

D^F;^ = ^|p|2 e^^^D^ip^ip^ (4.43) 



where the isovector 0'^ has been defined in Eq.( [4.22D . Under the assumptions 



mentioned in the above, at t = 0, this equation reads 

S'^F;^ = -g\p\'^d^e{x) {rf - n^0"0^) . (4.44) 

In general, the unit isovector cf)"" can be decomposed into 

(f) = cos 6* <^o + sin 6* n X (^0 + 2 sin^ - m ■ 0o) n (4.45) 

where (J)q = —(0,0,1). It is straightforward to verify that in the unitary 
gauge, (f) reduces to 0o- The relevant point in Eq.( [4.42| ) is that the versor 
n, about which it is performed the SU{2) gauge rotation, does not depend 
on the space coordinates. Therefore, without loosing generality, we have the 
freedom to choose n to be everywhere perpendicular to 0o- In other words, 
(j) can be everywhere obtained by rotating (pQ by an angle 6 in the plane 
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identified by n and 0o- Formally, $ = cos^0o + sin^ n x (pQ, which clearly 
describes a simple U{1) transformation. In fact, since it is evident that the 



condition fi ± 0o also implies n ± 0, the equation of motion ( |4.44| ) becomes 



a'^F;^ = -g\p\^d^e{x)n'' . (4.46) 

As expected, we see that only the gauge field component along the direction 
n, namely A^ = n'^PF^, that has some initial dynamics which is created by 
a nonvanishing gradient of the phase between the two domains. When we 
generalize this result to the full SU{2) x f/(l)y gauge structure, an extra 
generator, namely the hypercharge, comes-in. Therefore in this case is not 
any more possible to choose an arbitrary direction for the unit vector n since 
different orientations of the unit vector n with respect to 0o correspond to 
different physical situations. We can still consider the case in which n is 
parallel to 0o but we should take in mind that this is not the only possibility. 
In this case we have 



5%. = '-p\x){dJ + dM, (4.47) 

d^Fl = -lp\x){d,e + d,v) (4.48) 

where g and g' are respectively the SU{2) and U{1)y gauge coupling con- 
stants. It is noticeable that in this case the charged gauge fields are not 
excited by the phase gradients at the time when bubble first collide. We 
can combine Eqs. ([4.47|) and (|4.48| ) to obtain the equation of motion for the 
Z-boson field 



d'F^. = '^ ^^ p\x) [dj + d.v) . (4.49) 

This equation tells us that a gradient in the phases of the Higgs field gives 
rise to a nontrivial dynamics of the Z-field with an effective gauge coupling 



constant Jg'^ + g"^. We see that the equilibration of the phase {6 + cp) can 
be now treated in analogy to the U{1) toy model studied by Kibble and 
Vilenkin [|153|| , the role of the U{1) "electromagnetic" field being now played 
by the Z-field. However, differently from Ref. ||153|| the authors of Ref. [ [152| 



left the Higgs field modulus free to change in space. Therefore, the equation 
of motion of p(x) has to be added to ( |4.49| ). Assuming the charged gauge 
field does not evolve significantly, the complete set of equations of motion 



that we can write at finite, though small, time after the bubbles first contact, 
is 



^^ 2cos9 



w 



p^{x) (du<^ + 



9 



2 cos 9w 



d^d^ (p( 



)e*t + 2A p\x) - -T]' p{x)e'^ = 



(4.50) 



where d,, 



df, + i 



9 7 

2 cos ewf^^ 



is the quartic coupling. Note that, in analogy with 153 



7] is the vacuum expectation value of $ and A 

a gauge invariant 



phase difference can be introduced by making use of the covariant derivative 
dfj_. Equations ( [4. 501 ) are the Nielsen-Olesen equations of motion | |158| |. Their 
solution describes a Z- vortex where p = at its core ||161|| . The geometry 
of the problem implies that the vortex is closed, forming a ring which axis 
coincide with the conjunction of bubble centers. This result provides further 
support to the possibility that electroweak strings are produced during the 
EWPT. 

In principle, in order to determine the magnetic field produced during the 
process that we illustrated in the above, we need a gauge-invariant definition 
of the electromagnetic field strength in the presence of the non-trivial Higgs 
background. We know however that such definition is not unique [p.62| ]. For 
example, the authors of Ref . [P-52|| used the definition given in Eq.( [4.23D to 
find that the electric current is 



a'^ FJ = 2 tan 9wd^ {Z^d, In p{x) - Z,d^ In p(x)) 



(4.51) 



whereas other authors | |163|| , using the definition 



J^r = - sin ^v^r F» + cos OwF^,, + 



flU 



fj.f 



flU 



'-^^^e'^'"^r{D,^T{Djy , (4.52) 
9 



found no electric current, hence no magnetic field, at all. We have to observe, 
however, that the choice between these, as others, gauge invariant definitions 
is more a matter of taste than physics. Different definitions just give the 
same name to different combinations of the gauge fields. The important 
requirement which acceptable definitions of the electromagnetic field have to 
fulfill is that they have to reproduce the standard definition in the broken 
phase with a uniform Higgs background. This requirement is fulfilled by 
both the definitions used in the Refs. ||152| and ||163|| . In our opinion, it is 
not really meaningful to ask what is the electromagnetic field inside, or very 



close to, the electroweak strings. The physically relevant question is what 
are the electromagnetic relics of the electroweak strings once the EWPT is 
concluded. 

One important point to take in mind is that electroweak strings are not 



topologically stable (see [p.6CI|] and references therein) and that, for the phys- 
ical value of the Weinberg angle, they rapidly decay after their formation. 
Depending on the nature of the decay process two scenarios are possible. 



According to Vachaspati ||164|| long strings should decay in short segments 



of length ~ m^. Since the Z-string carry a flux of Z-magnetic flux in its 

interior 

Att 47r 
(pz = — = — sm9w cos 9w ■ (4.53) 

a e 

and the Z gauge field is a linear superposition of the W^ and Y fields then, 
when the string terminates, the Y flux cannot terminate because it is a U{1) 
gauge field and the Y magnetic field is divergenceless. Therefore some field 
must continue even beyond the end of the string. This has to be the massless 
field of the theory, that is, the electromagnetic field. In some sense, a finite 



segment of Z-string terminates on magnetic monopoles ||159|| . The magnetic 
flux emanating from a monopole is: 

Att Att 

(pB = — tan 6w = — siii^ 9w ■ (4.54) 

a e 

This flux may remain frozen-in the surrounding plasma and become a seed 
for cosmological magnetic flelds. 

Another possibility is that Z-strings decay by the formation of a W- 



condensate in their cores. In fact, it was showed by Perkins ||165|| that while 
electroweak symmetry restoration in the core of the string reduces mw, the 
magnetic fleld via its coupling to the anomalous magnetic moment of the 
VT-fleld, causes, for eB > m^, the formation of a condensate of the W- 
flelds. Such a process is based on the Ambjorn-Olesen instability which will 
be discussed in some detail in the Chap.^ of this review. As noted in [p.52| 



the presence of an inhomogeneous VT-condensate produced by string decay 
gives rise to electric currents which may sustain magnetic flelds even after 
the Z string has disappeared. The formation of a VT-condensate by strong 
magnetic flelds at the EWPT time, was also considered by Olesen ||166| . 



We can now wonder what is the predicted strength of the magnetic flelds 
at the end of the EWPT. An attempt to answer to this question has been 
done by Ahonen and Enqvist [ 167| | (see also Ref.[ |168|| ) where the formation 
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of ring-like magnetic fields in collisions of bubbles of broken phase in an 
Abelian Higgs model were inspected. Under the assumption that magnetic 



fields are generated by a process that resembles the Kibble and Vilenkin [ p,53| 
mechanism, it was concluded that a magnetic field of the order of i? ~ 2x 10^° 
G with a coherence length of about 10^ GeV~^ may be originated. Assuming 



turbulent enhancement the authors of Ref . ||167|| of the field by inverse cascade 



| 5l| , a root-mean-square value of the magnetic field -Brms — 10 ^^ G on a 



comoving scale of 10 Mpc might be present today. Although our previous 
considerations give some partial support to the scenario advocated in [ p.67| 



we have to stress, however, that only in some restricted cases it is possible to 
reduce the dynamics of the system to the dynamics of a simple U{1) Abelian 
group. Furthermore, once Z-vortices are formed the non-Abelian nature of 
the electroweak theory shows due to the back-reaction of the magnetic field 
on the charged gauge bosons and it is not evident that the same numerical 



values obtained in ||167|| will be obtained in the case of the EWPT. 

However the most serious problem with the kind of scenario discussed in 
this section comes form the fact that, within the framework of the standard 
model, a first order EWPT seems to be incompatible with the Higgs mass 
experimental lower limit ||144|| . Although some parameter choice of the min- 



imal supersymmetric standard model (MSSM) may still allow a first order 
transition ||145||, which may give rise to magnetic fields in a way similar to 



that discussed in the above, we think it is worthwhile to keep an open mind 
and consider what may happen in the case of a second order transition or 
even in the case of a cross over. 

The case of a second order EWPT 

As we discussed in the first part of this section, magnetic fields generation 
by Higgs field equilibration share several common aspects with the formation 
of topological defects in the early Universe. This analogy holds, and it is even 
more evident, in the case of a second order transition. The theory of defect 
formation during a second order phase transition was developed in a seminal 
paper by Kibble ||154|] . We shortly review some relevant aspects of the Kibble 



mechanism. We start from the Universe being in the unbroken phase of a 
given symmetry group G. As the Universe cools and approach the critical 
temperature Tc protodomains are formed by thermal fluctuations where the 
vacuum is in one of the degenerate, classically equivalent, broken symmetry 
vacuum states. Let M be the manifold of the broken symmetry degenerate 
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vacua. The protodomains size is determined by the Higgs field correlation 
function. Protodomains become stable to thermal fluctuations when their 
free energy becomes larger than the temperature. The temperature at which 
this happens is usually named Ginsburg temperature To- Below To stable 
domains are formed which, in the case of a topologically nontrivial manifold 
M, give rise to defect production. Rather, if M is topologically trivial, 
phase equilibration will continue until the Higgs field is uniform everywhere. 
This is the case of the Weinberg-Salam model, as well as of its minimal 
supersymmetrical extension. 

Higgs phase equilibration, which occurs when stable domains merge, gives 



rise to magnetic fields in a way similar to that described by Vachaspati [ p.06| 
(see the beginning of this section). One should keep in mind, however, that 
as a matter of principle, the domain size, which determine the Higgs field 
gradient, is different from the correlation length at the critical temperature 



152|| . At the time when stable domains form, their size is given by the 
correlation length in the broken phase at the Ginsburg temperature. This 
temperature was computed, in the case of the EWPT, by the authors of 
Ref . |p.52|| by comparing the expansion rate of the Universe with the nucleation 
rate per unit volume of sub-critical bubbles of symmetric phase (with size 
equal to the correlation length in the broken phase) given by 

*-b 

where £b is the correlation length in the broken phase. S^^ is the high tem- 
perature limit of the Euclidean action (see e.g. Ref. ||169| |). It was shown that 
for the EWPT the Ginsburg temperature is very close to the critical tem- 
perature, Tg = Tc within a few percent. The corresponding size of a broken 
phase domain is determined by the correlation length in the broken phase at 
T = Tg 

j7^ = V"mTG)),Ta) (4.56) 

where V{(f), T) is the effective Higgs potential. (i{TG)'i is weakly dependent on 
Mh, ihiTc) ^ 11/Tg for Mh = 100 GeV and 4(^0) - 10/Tg for Mh = 200 
GeV. Using this result and Eq.( [4.23| ) the authors of Ref. [p.52|| estimated the 
magnetic field strength at the end of the EWPT to be of order of 

Bi ~ 4e-^ sin=^ OwiliTc) ~ 10^^ G , (4.57) 
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on a length scale ihiTc). 

Although it was showed by Martin and Davis ||170|| that magnetic fields 
produced on such a scale may be stable against thermal fiuctuations, it is 
clear that magnetic fields of phenomenological interest live on scales much 
larger than it,{TG)- Therefore, some kind of average is required. We are ready 
to return to the discussion of the Vachaspati mechanism for magnetic field 



generation ||106|| . Let us suppose we are interested in the magnetic field on a 



scale L = Ni. Vachaspati argued that, since the Higgs field is uncorrelated 
on scales larger than i, its gradient executes a random walk as we move along 
a line crossing A^ domains. Therefore, the average of the gradient -0^^$ over 
this path should scale as viV. Since the magnetic field is proportional to the 
product of two covariant derivatives, see Eg. ( [4. 231 ), Vachaspati concluded 



that it scales as 1/A^. This conclusion, however, overlooks the difference 
between (D^$''")(D^$) and (D^$^D^$). This point was noticed by Enqvist 
and Olesen ||107|| (see also Ref. [[L09|| ) who produced a different estimate for 
the average magnetic field , (-B)rms, l = B{L) ~ Bij \fN . Neglecting the 
possible role of the magnetic helicity (see the next section) and of possible 
related effects, e.g. inverse cascade, and using Eq. (|4.57|) , the hne-averaged 
field today on a scale L ~ 1 Mpc (A^ ~ 10^^) is found to be of the order 
Bo(lMpc) ~ 10-21 G. 

Another important point of this kind of scenario (for the reasons which 
will become clear in the next section) is that it naturally gives rise to a 
nonvanishing vorticity. This point can be understood by the analogy with 
the process which lead to the formation of superfiuid circulation in a Bose- 
Einstein fiuid which is rapidly taken below the critical point by a pressure 
quench [p.71|| . Consider a circular closed path through the superfiuid of length 



C = 2'kR. This path will cross A^ ~ C/£ domains, where £ is the characteris- 
tic size of a single domain. Assuming that the phase 9 of the condensate wave 
function is uncorrelated in each of the A^ domains (random-walk hypothesis) 
the typical mismatch of 9 is given by: 

59= f V9 ■ ds r^ y/N (4.58) 

where V9 is the phase gradient across two adjacent domains and ds is the line 
element along the circumference. It is well known (see e.g. ||172|| ) that from 



the Schrodinger equation it follows that the velocity of a superfiuid is given 
by the gradient of the phase trough the relation v^ = {h/m)'V9, therefore 
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( 4.58| ) implies 

Vs = i^M^4^ ■ (4-59) 



It was argued by Zurek ||171|| that this phenomenon can effectively simulate 



the formation of defects in the early Universe. As we discussed in the previ- 
ous section, although the standard model does not allows topological defects, 
embedded defects, namely electroweak strings, may be produced through a 



similar mechanism. Indeed a close analogy was showed to exist | 173 | between 
the EWPT and the ^He superfluid transition where formation of vortices is 
experimentally observed. This hypothesis received further support by some 
recent lattice simulations which showed evidence for the formation of a clus- 



ter of Z-strings just above the cross-over temperature ||174|| in the case of a 
3D SU{2) Higgs model. Electroweak strings should lead to the generation 
of magnetic fields in the same way we discussed in the case of a first order 
EWPT. Unfortunately, to estimate the strength of the magnetic field pro- 
duced by this mechanism requires the knowledge of the string density and 
net helicity which, so far, are rather unknown quantities. 



4.4 Magnetic helicity and electroweak baryo- 
genesis 

As we discussed in the introduction of our report, the cosmological magnetic 
flux is a nearly conserved quantity due to the high conductivity of the Uni- 
verse. In this section we will focus on another quantity which, for the same 
reason, is approximatively conserved during most of the Universe evolution. 
This is the so called magnetic helicity defined by 

n= f d^xA-B (4.60) 

where A is the electromagnetic vector potential and B is the magnetic field. 
In the presence of a small value of the electric conductivity a the time evo- 
lution of H is given by 

^ = -- fd^xB-WxB . (4.61) 

dt a J 

Besides for the fact to be a nearly conserved quantity, the magnetic he- 
licity is a very interesting quantity for a number of different reasons. The 
main among these reasons are: 
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• In a field theory language 7i coincides with the Chern-Simon number 
which is know to be related to the topological properties of the gauge 
fields. 

• Since 7i is a P (parity) and CP-odd function, the observation of a 
nonvanishing net value of this quantity would be a manifestation of a 
macroscopic violation of both these symmetries. 

• It is know from magneto-hydro-dynamics that the presence of magnetic 
helicity can lead to the amplification of magnetic fields and contribute 
to their self-organization into a large scale ordered configuration (see 
Sec. |1.4| ). The same phenomenon could take place at a cosmological 
level. 

In the last few years, several authors proposed mechanisms for the produc- 
tion of magnetic helicity in the early Universe starting from particle physics 
processes. Cornwall ||175|| suggested that magnetic helicity was initially stored 



in the Universe under the form of baryons (B) and leptons (L) numbers pos- 
sibly generated by some GUT scale baryogenesis mechanism. He assumed 
that an order one fraction of the total classically conserved B + L charge was 
dissipated by anomalous processes at the EW scale and showed that a small 
fraction of this dissipated charge, of the order of ub+lT'^, may have been 
converted into a magnetic helicity of the order of 

n ~ a-\NB + Nl) ~ lO^*^ erg cm . (4.62) 

Another possibility is that before symmetry breaking of a non-Abelian gauge 
symmetry vacuum configurations existed which carried nonvanishing winding 
number. It was shown by Jackiw and Pi ||176|| that after symmetry breaking. 



one direction in isospin space is identified with electromagnet ism, and the 
projection of the vacuum configuration becomes a magnetic field with non 
vanishing helicity. 

A different mechanism was proposed by Joyce and Shaposhnikov ||177 



In this case it was assumed that some excess of right-handed electrons over 
left-handed positrons was produced by some means {e.g. from some GUT 
scale leptogenesis) above a temperature Tr. At temperatures higher than Tr 
perturbative processes which changes electron chirality are out of thermal 
equilibrium (Tr ~ 3 TeV in the SM [|178|| ). Therefore, a chemical potential 
for right electrons fiR can be introduced above Tr. On the other hand, the 
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corresponding charge is not conserved because of the Abehan anomaly, which 
gives 

d,fR = -0 U'' ■ (4.63) 

In the above, f^^, and f^'^ are, respectively, the U{1)y hypercharge field 
strength and its dual, g' is the associated gauge coupling constant, and 
yji = —2 is the hypercharge of the right electron. As it is well known, 
Eg. ( 4.63 ) relates the variation in the number of the right handed electrons 



Nji to the variation of the topological properties (Chern-Simon number) of 
the hypercharge field configuration. By rewriting this expression in terms of 
the hypermagnetic By and of the hyperelectric Ey fields, 

d,fR = -^ By . Ey , (4.64) 

it is evident the relation of j^ with the hypermagnetic helicity. It is worth- 
while to observe that only the hypermagnetic helicity is coupled to the 
fermion number by the chiral anomaly whereas such a coupling is absent 
for the Maxwell magnetic helicity because of the vector-like coupling of the 
electromagnetic field to fermions. From Eq. (|4.63| ) it follows that the variation 
in Nji, is related to the variation in the Chern-Simon number, 

Ncs = -tI^ / rf'xe,,fc^-6fc, (4.65) 

by ANji = ^yj^ANcs- In the above bk represents the hypercharge field 
potential. The energy density sitting in right electrons is of order fiJiT"^ 
and their number density of order UrT'^. Such fermionic number can be 
reprocessed into hypermagnetic helicity of order g''^kb'^, with energy of order 
k'^b'^, where k is the momentum of the classical hypercharge field and b is 
its amplitude. Therefore, at 6 > T/g'"^ it is energetically convenient for the 
system to produce hypermagnetic helicity by "eating-up" fermions. It was 
showed by Joyce and Shaposhnikov | 177|| , and in more detail by Giovannini 



and Shaposhnikov ||179|| , that such a phenomenon corresponds to a magnetic 
dynamo instability. In fact, by adding the anomaly term to the Maxwell like 
equations for the hyperelectric and hypermagnetic fields these authors were 
able to write the anomalous magneto-hydrodynamical (AMHD) equations for 
the electroweak plasma in the expanding Universe, including the following 
generalized hypermagnetic diffusivity equation 

^r^ = V X (/iijBy) + V X (v X By) + -V^By . (4.66) 

OT 7ia a 
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In the above a is the Universe scale factor, a' = (y''^/47r, r = J a~^{t)dt is the 
conformal time, and a the electric conductivity of the electroweak plasma 
9[. By comparing Eq.( |4.66| ) with the usual magnetic diffusivity equation 
see e.g. Ref. |p.5||) one sees that the first term on the r.h.s. of Eq.( [4.66D 



corresponds to the so-called dynamo term which is know to be related to the 



vorticity in the plasma flux |TT|. We see that the fermion asymmetry, by 
providing a macroscopic parity violation, plays a similar role of that played 
by the fluid vorticity in the dynamo amplification of magnetic fields. In 
the scenario advocated in Ref. ||177| , |179|| , however, it is not clear what are 



the seed fields from which the dynamo amplification start from. Clearly, no 
amplification takes place if the initial value of the hypermagnetic field van- 
ishes. Perhaps, seeds field may have been provided by thermal fiuctuations 
or from a previous phase transition although this is a matter of speculation. 
According to Joyce and Shaposhnikov ||177|| , assuming that a large right elec- 



tron asymmetry fin/T ~ 10~^ was present when T = Tr, magnetic field of 
strength B ~ 10^^ G may have survived until the EWPT time with typical 
inhomogeneity scales ~ 10^/T. 



Another interesting point raised by Giovannini and Shaposhnikov [ |179| 



is that the Abelian anomaly may also process a preexisting hypermagnetic 
helicity into fermions. In this sense the presence of tangled magnetic fields 
in the early Universe may provide a new leptogenesis scenario. 

Indeed, assuming that a primordial hypermagnetic field By was present 
before the EWPT with some non-trivial topology (i.e. (By ■ V x By) 7^ 0) 
the kinetic equation of right electrons for T > Tc is 

^ (^) = - f ,— By • V X By - (F + F.,)^, (4.67) 

dt\T J A-K^aaT^ 88 r v ^ np; ^ , v ; 

where 

783 a'2 |By|2 

is the rate of the Nji nonconserving anomalous processes whereas F is the 
rate of the perturbative ones. In the case F„p > F, as a consequence of 



Eq.( [4.67| ), one finds 



vr^ /By -V X By 




"""-TSy |B,f ' l+Q|^l' («9) 
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Below the critical temperature the hypermagnetic fields are converted into 
ordinary Maxwell magnetic fields . Similarly to the usual EW baryogenesis 
scenario, the fermion number asymmetry produced by the Abelian anomaly 
may survive the sphaleron wash-out only if the EWPT is strongly first order, 
which we know to be incompatible with the standard model in the absence of 
primordial magnetic fields. However, Giovannini and Shaposhinikov argued 
that this argument might not apply in the presence of strong magnetic fields 
(we shall discuss this issue in Chap.^). If this is the case a baryon asymmetry 
compatible with the observations might have been generated at the EW 
scale. Another prediction of this scenario is the production of strong density 
fiuctuations at the BBN time which may affect the primordial synthesis of 



light elements ||180 



Primordial magnetic fields and the primordial magnetic helicity may also 
have been produced by the interaction of the hypercharge component of the 
electromagnetic field with a cosmic pseudoscalar field condensate which pro- 
vides the required macroscopic parity violation. This idea was first sketched 



by Turner and Widrow |^ in the framework of an infiationary model of the 
Universe which we shall discuss in more details in the next section. Turner 
and Widrow assumed the pseudoscalar field to be the axion, a particle which 
existence is invoked for the solution of the strong CP problem (for a review 
see Ref. [ |181| ). Although the axion is supposed to be electrically neutral 



it couples to electromagnetic field by means of the anomaly. Indeed, the 
effective Lagrangian for axion electrodynamics is 

£ = -U.ed^^e - \f^,F>^' + gaOF^^F'^'' , (4.70) 

where ga is a coupling constant of order a, the vacuum angle 9 = (pa/ fa, 
(pa is the axion field, fa is the Peccei-Quinn symmetry breaking scale (see 
Ref. ||181|| ), F^t, is the electromagnetic field strength and F^'^ is its dual. Since 
the axion field, as any other scalar field, is not conformally invariant (see the 
next section), it will be amplified during the infiationary expansion of the 
Universe starting from quantum fiuctuations, giving rise to {9^) ~ {Hq/ faY, 
which can act as a source term for the electromagnetic field |. 

^The careful reader may wonder what is the fate of axions in the presence of cosmic 
magnetic fields. Interestingly, it was showed by Ahonen, Enqvist and Raffelt |182| that 
although oscillating cosmic axions drives an oscillating electric field, the ensuing dissipation 
of axions is found to be inversely proportional to the plasma conductivity and is, therefore, 
negligible. 
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Carroll and Field [TEZ] reconsider in more details the idea of Turner and 



Widrow and found that the evolution of a Fourier mode of the magnetic field 
with wave number k is governed by the equation 

^+(F±,„.f.)f,^0 (4,71) 

where F± = a^{Bx ± iBy) are the Fourier modes corresponding to different 
circular polarizations, and r is the conformal time. One or both polariza- 
tion modes will be unstable for k < ga\d(f)a/dT\, whereas both polarization 
modes can becomes unstable to exponential growth if is oscillating. In this 
case it seems as if a quite strong magnetic field could be produced during 
inflation. However, such a conclusion was recently criticized by Giovannini 
184|| who noted that above the EWPT temperature QCD sphalerons [[185| 



are in thermal equilibrium which can effectively damp axion oscillations. In 
fact, because of the presence of QCD sphaleron the axion equation of motion 
becomes 

0, + (3/7 + 7)4 = (4.72) 



where fl85 



^ sphal _ "s-t /^ -^gN 

J a J a 

(where Og = g'^/An). Giovannini founds that sphaleron induced damp- 
ing dominates over damping produced by the expansion of the Universe if 
fa > 10^ GeV. Since astrophysical and cosmological bounds | ]125| | leave open 



a window 10^° GeV < fa < 10^^ GeV, it follows that no magnetic fields 
amplification was possible until QCD sphaleron went out of thermal equilib- 
rium. A very tiny magnetic helicity production from axion oscillations may 
occur at lower temperatures. In fact, Giovannini [|184|| estimated that in the 



temperature range 1 GeV > T > lOMeV a magnetic helicity of the order of 

10"=^^ (4.74) 



(B-VxB) _ ^., 



O-Pc 

may be generated, which is probably too small to have any phenomenological 
relevance. 

Generation of magnetic fields from coherent oscillations or roll-down of 
pseudoscalar fields different from the axion has been also considered in the 
literature. It is interesting that pseudoscalar fields with an axion-like anoma- 
lous coupling to the electromagnetic field appear in several possible extension 
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of the SM. Typically these fields have only perturbative derivative interac- 
tions and therefore vanishing potential at high temperatures, and acquire a 
potential at lower temperatures through non-perturbative interactions. The 
generic form of this potential is V{(j)) = Vq A{(f)/f), where A is a bounded 
function of the "Peccei-Quinn" scale / which in this case can be as large 
as the Planck scale. The pseudoscalar mass m ~ Vq/J may range from a 
few eV to 10^^ GeV. The amplification of magnetic fields proceeds in a way 
quite similar to that discussed for axions. The evolution equation of the elec- 
tromagnetic Fourier modes were derived by Brustein and Oaknin [ p.86|| who 
added to Eq. ( |4.71| ) the effect of a finite electric conductivity, finding 

If the pseudoscalar field is oscillating, the field velocity d(f)/dT changes sign 
periodically and both polarization modes are amplified, each during a dif- 
ferent semi-cycle. Each mode is amplified during one part of the cycle and 
damped during the other part of the cycle. Net amplification results when 
gkd(f)/dT > 0"^ and the total amplification is exponential in the number of cy- 
cles. Taking into account that a ~ lOT before the EWPT [4^, Brustein and 



Oaknin found that a huge amplification can be obtained from pseudoscalar 
field oscillations at T ~ 1 TeV, for a scalar mass m of few TeV and gf ~ 10. 
This is a particularly interesting range of parameters as in this case the pseu- 
doscalar field dynamics could be associated with breaking of supersymmetry 
at the TeV scale. Since only a limited range of Fourier modes are ampli- 
fied, with k not too different from T, and modes with k/T < tewO' ^ 1 
are rapidly dissipated, amplified fields may survive until the EWPT only if 
amplification occurred just before the transition. If this is a natural assump- 
tion may be a matter of discussion. However, it was pointed-out by Brustein 



and Oaknin | 186 | that if, depending on the form of the pseudoscalar field 
potential, this field rolls instead of oscillating, then the hypermagnetic fields 
would survive until EWPT even if the amplification takes place at higher 
temperatures before the transition. In this case there could be interesting 
consequences for the EW baryogenesis. This is so, because the pseudoscalar 
field may carry a considerable amount of helicity. According to what dis- 
cussed in the above, this number will be released in the form of fermions and 
baryons if the EWPT is strongly first order. Brustein and Oaknin argued 
that this mechanism could naturally generate the observed BAU. 
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Clearly, a more serious problem of this mechanism is the same of other 
EW baryogenesis scenarios, namely to have a strongly first order EWPT. An 
interesting possibility which was proposed by two different groups ||179| , p.87|| 
is that strong magnetic fields may enhance the strength of the EWPT (see 
Chap.|^). Unfortunately, detailed lattice computations [|188|| showed that this 
is not the case. Furthermore, in a recent work by Comelli et al. | 189 | it was 
shown that strong magnetic fields also increase the rate of EW sphalerons 
so that the preservation of the baryon asymmetry calls for a much stronger 
phase transition than required in the absence of a magnetic field. The au- 



thors of Ref . ||189|| showed that this effect overwhelms the gain in the phase 
transition strength (see Sec. 5.3.2]) . Therefore, the only way for the kind of 
EW baryogenesis mechanism discussed in the above to work is to invoke for 
extensions of the standard model which allow for a strong first order EWPT 

ra. 



Electroweak strings, which we introduced in the previous section, may 
also carry a net hypermagnetic helicity and act as a source of the observed 
BAU. One of the interesting features of this objects is that they should have 
been formed during the EWPT even if this transition is second order or just 



a cross over 



It is known that electroweak strings can carry magnetic 
helicity, hence a baryon number, which is related to the twisting and linking 
of the string gauge fields |p,6CI| , |T9( 



Several authors (see e.g. Ref.[ |191|| ) tried 
to construct a viable EW baryogenesis scenario based on these embedded 
defects. Many of these models, however, run into the same problems of more 
conventional EW baryogenesis scenarios. An interesting attempt was done 
by Barriola | |192| ] who invoked the presence of a primordial magnetic field at 
the EWPT time. Barriola observed that the production and the following 
decay of electroweak strings give rise, in the presence of an external magnetic 
field, to a variation in the baryon number AB in the time interval dt, given 
by 



AS 



327r2 



a 



dt / d-^x a cos 26wEz ■ Bz -\ sin'' 9 



w 



[e-b 



z + Ez ■ B^ 

(4.76) 

where Ez and Bz are the Z-electric and Z-magnetic fields of the strings 
whereas E and B are the corresponding standard Maxwell fields. The first 
term in the r.h.s. of ( |4.76| ) represents the change in the helicity of the string 
network, and the second and third terms come form the coupling of the string 
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fields with the external Maxwell fields. Whereas, the first term averages to 
zero over a large number of strings, the other terms may not. Clearly, some 
bias is required to select a direction in the baryon number violation. In 
the specific model considered by Barriola this is obtained by a CP violation 
coming from the extension of the Higgs sector of the SM, and from the 
dynamics of strings which are supposed to collapse along their axis. It was 
concluded by Barriola that such a mechanism could account for the BAU. 
Unfortunately, it is quite unclear if the invoked out-of-equilibrium mechanism 
based on the string collapse could indeed be effective enough to avoid the 
sphaleron wash-out. We think, however, that this possibility deserves further 
study. 

4.5 Magnetic fields from infiation 

As noted by Turner and Widrow ^^ inflation (for a comprehensive intro- 



duction to inflation see Ref . ||46|| ) provides four important ingredients for the 
production of primeval magnetic flelds . 

• Inflation naturally produces effects on very large scales, larger than the 
Hubble horizon, starting from microphysical processes operating on a 
causally connected volume. If electromagnetic quantum fluctuations 
are amplifled during inflation they could appear today as large-scale 
static magnetic flelds (electric fleld should be screened by the high 
conductivity plasma). 

• Inflation also provides the dynamical means to amplify these long- 
wavelength waves. If the conformal invariance of the electromagnetic 
fleld is broken in some way (see below) magnetic flelds could be excited 
during the de Sitter Universe expansion. This phenomenon is analo- 
gous to particle production occurring in a rapidly changing space-time 
metric. 

• During inflation (and perhaps during most of reheating) the Universe 
is not a good conductor so that magnetic flux is not conserved and the 
ratio r of the magnetic fleld with the radiation energy densities can 
increase. 

• Classical fluctuations with wavelength A < H^^ of massless, minimally 
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coupled fields can grow super-adiabatically, i.e. their energy density 
decrease only as ~ a~^ rather than a~^. 

The main obstacle on the way of this nice scenario is given by the fact that 
in a confornially fiat metric, like the Robertson- Walker usually considered, 
the background gravitational field does not produce particles if the underlying 
theory is confornially invariant [p.93|| . This is the case for photons since the 



classical electrodynamics is confornially invariant in the limit of vanishing 
fermion masses. Several ways out this obstacle have been proposed. Turner 



and Widrow ^^ considered three possibilities. The first is to break explicitly 
conformal invariance by introducing a gravitational coupling, like RA^A^ or 
Rfj^^A^^A^, where R is the curvature scalar, Rfj_u is the Ricci tensor, and A'^ 
is the electromagnetic field. These terms breaks gauge invariance and give 
the photons an effective, time-dependent mass. In fact, one of the most 
severe constraints to this scenario come from the experimental upper limit 
to the photon mass, which today is ?Ti^ < 2 x lO"'^^ eV [ p.94|| . Turner and 
Widrow showed that for some suitable (though theoretically unmotivated) 
choice of the parameters, such a mechanism may give rise to galactic magnetic 
fields even without invoking the galactic dynamo. We leave to the reader to 
judge if such a booty deserve the abandonment of the theoretical prejudice in 
favor of gauge invariance. A different model invoking a spontaneous breaking 
of gauge symmetry of electromagnetism, implying nonconservation of the 
electric charge, in the early stage of the evolution of the Universe has been 



proposed by Dolgov and Silk [p.95| 



The breaking of the conformal invariance may also be produced by terms 
of the form R^j^yXf^F^^ F''^'^ / m? or RF^^F^^, where m is some mass scale re- 
quired by dimensional considerations. Such terms arise due to one-loop vac- 
uum polarization effects in curved space-time, and they have the virtue of 
being gauge invariant. Unfortunately, Turner and Widrow showed that they 
may account only for a far too small contribution to primordial magnetic 
fields. The third way to break conformal invariance discussed by Turner and 
Widrow invoke a coupling of the photon to a charged field which is not confor- 
nially coupled or the anomalous coupling to a pseudoscalar. This mechanism 
was already illustrated in the previous section. 

The anomaly can give rise to breaking of the conformal invariance also 
in a different way. The kind of anomaly we are now discussing about is the 
conformal anomaly, which is related to the triangle diagram connecting two 
photons to a graviton. It is known (for a review see Ref. ||196|| ) that this kind 
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of diagrams breaks conformal invariance by producing a nonvanishing trace 
of the energy-momentum tensor 

TH = ^C^" "" ' (4-77) 

where a is the fine-structure constant of the theory based on the SU{N) 
gauge-symmetry with Nf fermion famihes, and 

(3=^N-'^Nf. (4.78) 



Dolgov [ |197|| pointed-out that such an effect may lead to strong electromag- 



netic fields amplification during inflation. In fact, Maxwell equations are 
modified by the anomaly in the following way 

d^F!: + (3^Fi: = (4.79) 

which, in the Fourier space, gives rise to the equation 

A" + j3-A' + A;M = , (4.80) 

a 

where A is the amplitude of the vector potential, and a prime stands for 
a derivation respect to the conformal time r. At the inflationary stage. 



when a' /a = — 1/r Dolgov found a solution of ( |4.80| ) growing like f -— I 



Since k~^ grows well above the Hubble radius during the de Sitter phase, 
a huge amplification can be obtained ii (3 > 0. Dolgov showed that for 
/3 ~ 1 the magnetic field generated during the inflationary stage can be large 
enough to give rise to the observed fields in galaxies even without a dynamo 
amplification. Unfortunately, such a large value of j3 seems to be unrealistic 
(/? ~ 0.06 for SU(5) with three charged fermions). The conclusion is that 
galactic magnetic fields might be produced by this mechanism only invoking 
a group larger than SU{5) with a large number of fermion families, and 
certainly no without the help of dynamo amplification. 

As we discussed in the above conformal invariance of the electromagnetic 
field is generally spoiled whenever the electromagnetic field is coupled to a 
scalar field. Ratra ||198|| suggested that a coupling of the form e'^'^F'^^F^j^, 



where k is a arbitrary parameter, may lead to a huge amplification of electro- 
magnetic quantum fluctuations into large scale magnetic fields during infla- 
tion. Such a coupling is produced in some peculiar models of inflation with 
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an exponential inflaton potential [ T^ ] . It should be noted by the reader that 
the scalar field coincide here with the inflaton fleld. According to Ratra, 
present time intergalactic magnetic flelds as large as 10~^ G may be pro- 
duced by this mechanism which would not require any dynamo ampliflcation 
to account for the observed galactic fields. Unfortunately, depending on the 
parameter of the underling model, the predicted field could also be as low as 
10-65 G ! 

A slightly more predictive, and perhaps theoretically better motivated, 
model has been proposed independently by Lemoine and Lemoine [pOO|| , and 



Gasperini, Giovannini and Veneziano ||201|| , which is based on superstring 
cosmology |P02 , 203 ], This model is based on the consideration that in string 
theory the electromagnetic fields is coupled not only to the metric (g^u), but 
also to the dilaton field 0. In the low energy limit of the theory, and after 
dimensional reduction from 10 to 4 space-time dimensions, such a coupling 
takes the form 

y^e-'^F^^'^F,, (4.81) 

which breaks conformal invariance of the electromagnetic field and coincides 
with the coupling considered by Ratra ||198|| if k = —1. Ratra, however, 
assumed infiation to be driven by the scalar field potential, which is not 
the case in string cosmology. In fact, typical dilaton potentials are much 
too steep to produce the required slow-roll of the infiaton (= dilaton) field. 
According to string cosmologists, this problem can be solved by assuming 
infiation to be driven by the kinetic part of the dilaton field, i.e. from 0' 
||202|| . In such a scenario the Universe evolves from a fiat, cold, and weakly 
coupled (0 = — oo) initial unstable vacuum state toward a curved, dilaton- 
driven, strong coupling regime. During this period, called pre-big-bang phase, 
the scale factor and of the dilaton evolve as 



a{T) ~ (— t)'^, 0(r) ~ Klna, r < 



-n 



(4.82) 



with j3 > 1 and k < 0. At r > —ti it begins the standard FRW phase with 
a radiation dominated Universe. In the presence of the non-trivial dilaton 
background the modified Maxwell equations takes the form 



V^F^" 



V„,0 F'^'^ = 



(4.83) 



Electromagnetic field amplification from quantum fiuctuations take place 
during the pre-big-bang phase when 0' = 6/t, where 6 = (3hi. By follow- 
ing the evolution of the electromagnetic field modes from t = — cxd to now. 
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Lemoine and Lemoine estimated that, in the most simple model of dilaton- 
driven inflation (with V{(j)) = p = p = 0) a. very tiny magnetic field is 
predicted today 



{B') 



2\l/2 



10 



-62 



llMpci 



-2.1 






PI 



0.07 



G 



(4.84) 



where Hi = 1/ti, which is far too small to account for galactic fields. 

Gasperini et al. [|201|| reached a different conclusion, claiming that mag- 
netic fields as large those required to explain galactic fields without dynamo 
amplification may be produced on the protogalactic scale. The reason for 
such a different result is that they assumed a new phase to exist between the 
dilaton-dominated phase and the FRW phase during which dilaton potential 
is nonvanishing. The new phase, called the string phase, should start when 
the string length scale Xs becomes comparable to the horizon size at the 
conformal time ts [p04|| . Unfortunately, the duration of such a phase is quite 
unknown, which makes the model not very predictive. 

Recently, several papers have been published (see e.g. Refs. ||205] , p06| . 



which proposed the generation of magnetic fields by fluctuations of 
scalar (or pseudo-scalar) fields which were amplified during, or at the end, of 
infiation. In some of those papers |p05| , |207|| the authors claim that magnetic 
fields as strong as those required to initiate a successful galactic dynamo 
may be produced. This conclusion, however, is based on incorrect assump- 
tions. According to Giovannini and Shaposhnikov 



], the main prob- 
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lem resides in the approximate treatment of dissipative effects adopted in 
Refs. ||205] , |207|] . Whereas in Ref. ||207|| dissipation of electric fields was totally 



neglected, in Ref. |p05|| an incorrect dependence of the electric conductivity 
on the temperature was used. Adopting the correct expression for the con- 
ductivity (see e.g. Refs.|^) the authours of Ref. |P06|| found that infiation 



produces fields which are too small to seed galactic dynamo. 

4.6 Magnetic fields from cosmic strings 

At the beginning of this chapter we shortly discussed the Harrison- Rees mech- 
anism for vortical production of magnetic fields in the primordial plasma. In 
this section we will shortly review as cosmic string may implement this mech- 
anism by providing a vorticity source on scales comparable to galactic sizes. 
Cosmic strings are one-dimensional topological defects which are supposed 
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to have been formed during some primordial phase transition through the 
Kibble mechanism ||154|| , which we already discussed in Sec.^]3| (for a re- 



view on cosmic strings see Ref. ||209|| ). The idea that cosmic string may 



produce plasma vorticity and magnetic fields was first proposed by Vachas- 



pati and Vilenkin [plCI|] . In this scenario vorticity is generated in the wakes 



of fast moving cosmic strings after structure formation begins. Differently 



from Harrison's ||131|| scenario (see Sec. |4.1| ), in this case vorticity does not 
decay with Universe expansion since the vortical eddies are gravitationally 
bounded to the string. Even if the mechanism we are considering is supposed 
to take place after recombination, a sufficient amount of ionization should 



be produced by the violent turbulent motion so that the Harrison-Rees |jT^ 
mechanism can still operate. The scale of coherence of the generated mag- 
netic fields is set by the scale of wiggles of the string and, for wakes created at 
recombination time it can be up to 100 kpc. The predicted field strength is of 
the order of ~ 10~^^ G, that could enough to seed the galactic dynamo. The 
main problem with this scenario is that is not clear whether stable vortical 
motion can be really generated by the chaotic motion of the string wiggles. 
An alternative mechanism has been proposed by Avelino and Shellard 



21 If . In their model, vorticity is generated not by the wiggles but by the 
strings themselves which, because of the finite dynamical friction, drag mat- 
ter behind them inducing circular motions over inter-string scales. Unfor- 
tunately, the magnetic field strength predicted by this model at the present 
time is very weak ~ 10~^^ G, which can only marginally seed the galactic 
dynamo. 

Larger fields may be produced if cosmic strings are superconducting. 



String superconductivity was first conceived by Witten ||212|| . The charge 



carriers on these strings can be either fermions or bosons, and the critical 
currents can be as large as 10^° Ampere. If primordial magnetic fields pre- 
exist, or formed together with, cosmic string they may play a role charging 
up superconducting string loops and delaying their collapse [ pi3|| . Other- 



wise, superconducting cosmic strings can themselves give rise to magnetic 
fields in a way similar to that proposed by Avelino and Shellard. An impor- 
tant difference, however, arises with respect to the non-superconducting case 
discussed in Ref.[ ^ll|] . As shown by Dimopoulos and Davis | |214| |, supercon- 



ducting strings networks may be more tangled and slower than conventional 
cosmic strings because of the strong current which increases dynamical fric- 
tion. It was shown by Dimopoulos ||215|| that if the string velocity is small 
enough the gravitational influence of the string on the surrounding plasma 
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becomes a relevant effect. As a consequence, plasma is dragged by the string 
acquiring substantial momentum. Such a momentum may induce turbulence 
which could generate magnetic fields on scale of the order of the inter-string 
distance. Such a distance is smaller than the conventional cosmic string dis- 
tance. Quite strong magnetic fields may be produced by this mechanism. 
The only known constraint come from the requirement that the string net- 
work does not produce too large temperature anisotropies in the CMBR. By 
imposing this constraint, Dimopoulos estimated that present time magnetic 
fields as large as 10~^^ G with a coherence scale of ~ 1 Mpc may be produced. 
In a recent paper by Voloshin it was claimed that the generation of large 
scale magnetic fields by domain walls is not possible |pi(: . 
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Chapter 5 

Particles and their couplings in 
the presence of strong magnetic 
fields 



In the previous chapter we have seen that, very strong magnetic fields could 
have been produced in the early Universe. Here we investigate the effects 
of such strong fields on bound states of quarks and on condensates created 
by spontaneous symmetry breaking. For example, we have already seen 
in Chap.|] that strong magnetic fields can affect masses and decay rates of 
charged particles and modify the rate of weak processes. As we already 



discussed in Sec. |1.4| another crucial issue concerns the stability of strong 
magnetic fields. QED allows the existence of arbitrary large magnetic fields 
provided matter constituents have spin less than ^ ||38[| . 

This is so because the Lorentz force cannot perform any work on charged 
particles so that real particle-antiparticles free pairs cannot be produced. 



In Sec. |3.1| we have already seen that quantum corrections do not spoil this 
classical argument. We have also seen that although pair production can 
be catalyzed by strong magnetic fields at finite temperature and density, 
in this case it is the heat-bath that pays for the energy cost of the effect. 
This situation changes when, in the presence of very strong fields, QCD 
and electroweak corrections cannot be disregarded. We shall see that the 
QCD and electroweak field allow for the formation of condensates of charged 
particles with no energy cost. This may lead to screening of magnetic, or 
hypermagnetic, fields resembling the Meissner effect in superconductors. 



108 



5.1 Low Lying States for Particles in Uniform 
Magnetic Fields 

Following Refs.piT], pi8|| in this section we will consider questions related 



to the mass shifts and decay of bound states of quarks. Mass shifts occur 
both due to the effect that magnetic fields have on the strong binding forces, 
and due to the direct interactions of charged spinning particles with external 
fields. The modifications of the strong forces are such as to close the gap 
between the proton and neutron masses and ultimately make the proton 
heavier. A delicate interplay between the anomalous magnetic moments of 
the proton and neutron drives the mass shifts due to the direct interactions 
in the same direction. For B > 1.5 x 10^^ G the neutron becomes stable and 
as the field is increased past 2.7 x 10^^ G the proton becomes unstable to 
decay into a neutron, positron and neutrino. 

The quantum mechanics of a Dirac particle with no anomalous magnetic 
moment in a uniform external magnetic field is straightforward. We shall 
present the results for the case where particles do have such anomalous mo- 
ments. In reality, in fields so strong that the mass shifts induced by such 
fields are of the order of the mass itself one cannot define a magnetic mo- 



ment as the energies are no longer linear in the external field. Schwinger ^7\ 
calculated the self energy of an electron in an external field and we shall use 
here his results. We cannot follow this procedure for the proton or neutron 
as we do not have a good field theory calculation of the magnetic moments of 
these particles, even for small magnetic fields; all we have at hand is a phe- 
nomeno logical anomalous magnetic moment. However, for fields that change 
the energies of these particles by only a few percent, we will consider these 
as point particles with the given anomalous moments. In the Sec. |5.1.7] we 
will discuss possible limitations of this approach. 

5.1.1 Protons in an External Field 

The Dirac Hamiltonian for a proton with a uniform external magnetic field 
B is 

iJ = a • (p - eA(r)) + /5M, - ^ (| - l) /55] ■ B . (5.1) 
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The vector potential A(r) is related to the magnetic field by A(r) = |r x B 
and gp = 5.58 is the proton's Lande g factor. We first solve this equation for 
the case where the momentum along the magnetic field direction is zero and 
then boost along that direction till we obtain the desired momentum. For B 



along the z direction and Pz = the energy levels are ||3§ 



Er, 



2eB(n + 



eBs + M„ 



2Mr. 



9p 



l]Bs 



(5.2) 



In the above, n denotes the Landau level, m the orbital angular momentum 
about the magnetic field direction and s = ±1 indicates whether the spin is 
along or opposed to that direction; the levels are degenerate in m. n = and 
s = +1 yield the lowest energy 



E = M„ 



M„ 

^ 2M„ 



9p 
2 



DB. 



(5.3) 



As we shall be interested in these states only we will drop the n and s 
quantum numbers. The Dirac wave function for this state is 



^: 



m,pz 



=o(rj 






\o J 



(pm{x,y) 



(5.4) 



's are the standard wave functions of the lowest Landau level; 



(pm{x,y) 



||ei?| 



m + l 



vvrm! 



[x + iy]"^ exp 



4 



\eB\{x'^ + y^) 



Boosting to a finite value of pz is straightforward; we obtain 



(5.5) 



EmiPz) = \JPz^ + M2 



(5.6) 



with a wave function 
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(5.7) 
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where 26, the rapidity, is obtained from tanh26' = Pzl Emipz)- 
In the non-relativistic hmit the energy becomes 

2 



Em{Pz 

and the wave function reduces to 
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2M 
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(5.8) 



(5.9) 



5.1.2 Neutrons in an External Field 

For a neutron the Dirac Hamiltonian is somewhat simpler 

/f = «.p + ,aM„-^(|);3E.B. 

with Qn = —3.82. Again for Pz = the states of lowest energy, the ones we 
shall be interested in, have energies 



(5.10) 



E{p^,Pz = 0) 
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Boosting to a finite Pz "we obtain 

E{v) = \/^(p±,P. = 0)'+p.2, 
The wave functions corresponding to this energy are 



V^p(r) 



(2vr)^ 



■u[p,s 



(5.11) 



(5.12) 



(5.13) 



where u(p, s = —1) is the standard spinor for a particle with momentum p. 



energy -i/p^ + M„^ (not the energy of Eq. (|5.12| )) and spin down. 
In the non-relativistic limit 



E(p) = M„ + 



2M„. 



and the wave functions are 



V^p(r) 
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(27r 



(5.15) 



111 



5.1.3 Electrons in an External Field 



We might be tempted to use, for the electron, the formahsm used for the 
proton with the Lande factor replaced hj Qe = 2 + a/ir. However as we shall 
see for magnetic fields sufficiently strong as to make the proton heavier than 
the neutron the change in energy of the electron would appear to be larger 
than the mass of the electron itself. The point particle formalism breaks down 
and we have solve QED, to one loop, in a strong magnetic field; fortunately 
this problem was treated by Schwinger[P^]. The energy of an electron with 
Pz = 0, spin up and in the lowest Landau level is 



E,^ 



m,pz=0 



M, 



■ a , f2eB\ 

IH In — ^ 

27r \M,^J_ 



(5.16) 



For field strengths of subsequent interest this correction is negligible; the 
energy of an electron in the lowest Landau level, with spin down and a 
momentum of Pz is 

' ' (5.17) 
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and with wave function similar to those of the proton 
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(5.18) 



where the boost rapidity, 26, is defined below Eq. ( |5.7] ) while the Landau 
level wave function is defined in Eq. ( ^.51) . The reason the complex conjugate 
wave function appears is that the electron charge is opposite to that of the 
proton. 



5.1.4 Effects of Magnetic Fields on Strong Forces 

We must be sure that shifts due to changes of color strong forces will not 
shift states in the opposite direction. The best method to study masses of 



QCD bound states is the use of sum rules ||219|| . This method uses the SVZ 
|220|| generalized short distance expansion that includes not only perturba- 
tive pieces, but also higher dimensional operators like the chiral and gluon 
condensates refiecting the non Abelian nature of the vacuum. Fortunately 
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the proton has a simple structure ||219|| which reflects the fact that if chiral 
symmetry is restored the proton and neutron masses vanish. 

1 
Mp^n = 3a {qq)^ + small corrections (5.19) 

where a is a constant. Meson mass terms are more involved; for example the 
p mass is 

Mp = b (perturbative terms) + c {G^^G^") + d {qq) . (5.20) 

b, c and d are constants of comparable magnitude ||219|| . As we shall show 
it is only the change of (qq) due to external magnetic fields that may be 
obtained in a reliable manner. 

In the presence of external fields we expect the chiral condensates for 
quarks of different charges to vary and Eq. ( |5.19| ) becomes 

M^ = a (2 (mm) + (rfJ)) , 

M^ = a (2{dd) + {uu)) . (5.21) 

To first order in condensate changes we find 



•-P 



9 V {uu) (dd) J 



Combining 



A simple method for studying the behavior of the chiral condensates in the 
presence of external constant fields is through the use of the Nambu-Jona- 



Lasinio model ||221|| . This has been done by Klevansky and Lemmer [p22| 
and a fit to their results is 



{qq)iB) = {qq)iO) 



e,B^^- 



(5.24) 



A2 
with A = 270 MeV and e^ the charge on the quark. To lowest order we find 

^<'«> - ' ^''^^' . (5,25) 
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and 

5M, - 5M„ = ^ (^^ j ; (5.26) 

As in the previous section, these corrections are such as drive the proton 
energy up faster than that of the neutron. One can understand the sign of 
this effect; the radius of a quark-antiquark pair will decrease with increasing 
magnetic field thus making the condensate larger. As the u quark has twice 
the charge of the d quark, its condensate will grow faster and as there are 
more u quarks in the proton than in the neutron its mass will increase faster. 
The fact that our estimate of the sign of the neutron-proton mass differ- 
ence is the same as that due to electromagnetic effects is crucial. QCD sum 
rules and our method of evaluating the chiral condensates are both crude 
and the magnitude of the mass difference is uncertain. Had the sign of the 
hadronic correction been opposite, cancellations could have occurred and the 
argument for a narrowing of the mass difference and ultimate reversal could 
not have been made. We are quite sure of the sign of the chiral change. The 
magnetic field acts in the naive way in the spin of the scalar bound state in 



the condensate so it is sure that goes as calculated [223 



5.1.5 Proton Life Time 

First we study the decay kinematics. 

Combining Eq. ( p.3|) , Eq. ( |5.14|) and Eq. ( |5.26| ) we find the proton-neutron 



energy difference as a function of the applied magnetic field, 

A(5) = (-1.3 + 0.385i4 + O.llBfA MeV . (5.27) 



i?i4 is the strength of the magnetic field in units of 10^^ T (1 T = 10^ G). 
The neutron becomes stable for B > 1.5 x 10^^ T and the proton becomes 
unstable to /3 decay for B > 2.7 x 10^^ T. We shall now turn to a calculation 
of the life time of the proton for fields satisfying the last inequality. 

With the wave functions of the various particles in the magnetic fields we 
may define field operators for these particles. For the proton and electron we 
shall restrict the summation over states to the lowest Landau levels with spin 
up, down respectively; for magnetic fields of interest the other states will not 
contribute to the calculation of decay properties. For the same reason, the 
neutron field will be restricted to spin down only. The proton and neutron 
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kinematics will be taken as non-relativistic. 



'^pij)=T.j dPz 
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(5.28) 



with (j)m{x,y) defined in Eq. { f).5\) and the energy, E^ipz) in Eq- ( ^-81) . 
O'miPz) is the annihilation operator for a proton with momentum p^z and 
angular momentum m; bm{Pz) is the same for the negative energy states. For 
the neutron the field is 



^Jr) 



d^p 



a[p) 



1 


voy 



ipr 



(27r) 



I + b\p) 






viy 



g-ipr 



(2^: 



(5.29) 



with an obvious definition of the annihilation operators. For the electron we 
use fully relativistic kinematics and the field is 



*eW-E/*.Jf 



am{p 



+bLiPz 
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V sinh^ y 



3«PzZ 



/27r 
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(5.30) 



5.1.6 Decay Rates and Spectrum 

The part of the weak Hamiltonian responsible for the decay p 

is 

G 



n + e^ + Ue 



H = ^J d^'x^nl.il + 75)^p^.7'^(l + 75)^e . (5.31) 
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For non-relativistic heavy particles the matrix element of this Hamiltonian 
between a proton with quantum numbers p^ = , m = rrii, a neutron 
with momentum p„, a neutrino with momentum p^^ and an electron in state 
m = ruf and with pz,e is 



1 



2Gp E,+pz^e Y . .. ,.. M 



/ 



dx dy (f)*^^{x,y)(j)mi{x,y) exp[-i{p±^r^ + P±,..) ■ r- J ; (5.32) 



6i, is the azimuthal angle of the neutrino. The integral in the above expression 
can be evaluated in a multipole expansion. Note that the natural extent of the 
integral in the transverse direction is 1/y/eB whereas the neutron momenta 



are, from Eq. (|5.27|) , of the order of JeB{0.12 + 0.04i?i4); thus setting the 



exponential term in this integral equal to one will yield a good estimate for 
the rate and spectrum of this decay. The positron spectrum is given by 

dr _4G^M,E.+p.,^^_^^^3. (5.33) 



dp,,e 3 (27r)6 Ee 

where A is defined in Eq. ( p.27| ). For A ^ Mg the total rate is easily obtained 



_ 2 G|M 4 



For 5 = 5 X 10^^ T, the lifetime is r = 6 s. 



5.1.7 Caveats and Limitations 



For general magnetic fields we expect the masses of particles to be nonlinear 
functions of these fields. Such an expression has been obtained, to order 



a, for the electron[0. For small fields this reduces to a power series, up 
to logarithmic terms, in B/Bc, where Be is some scale. For the electron 
Be = m1/e. For the hadronic case the value of Bq is uncertain. B^ = 
Mp/e = 1.7 X 10^^ T is probably too large and Mp should be replaced by 
a quark constituent mass and e be e^; in that case i?c = (2 — 4) x 10^^ T, 
depending on the quark type. This is also the range of values of A^/cg in 
Eq. ( p.24| ). The effects we have studied need fields around a few xlO^^ T or 



an order of magnitude smaller than the lowest candidate for B^.. Eq. ( 5.27 ) 
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bay be viewed as a power series expansion up to terms of order (B/Bc)"^; as 
the coefficient of the quadratic term was obtained from a fit to a numerical 



solution, logarithms of B/ Be may be hidden in the coefficient. As in Ref. p7| , 
even powers will be spin independent and the odd ones will be linear in 
the spin direction and may be viewed as field dependent corrections to the 
magnetic moment. We cannot prove, but only hope, that the coefficient of the 
(B/Bc)^ term, the ffist correction to the magnetic moment is not unusually 
large; should it turn out to be big and of opposite sign to the linear and 
quadratic terms, the conclusions of this analysis would be invalidated. These 
arguments, probably, apply best to the field dependence of the magnetic 
moments of the quarks rather than the total moment of the baryons. We 
may ask what is the effect on these magnetic moments due to changes in 
the "orbital" part of the quark wave functions. To ffist order we expect no 
effect as all the quarks are in S states and there is no orbital contribution to 
the total moment. The next order perturbation correction will be down by 
(rfy/rc)^ compared to the leading effect; Vh is a hadronic radius and r^ is the 
quarks cyclotron radius. This again contributes to the {B/ B^^ term in the 
expansion for the energy of a baryon. 



Another limitation is due to the results of Ref. ||218|| where it is shown 



that fields of the order of a few x 10^^ T are screened by a changes in chiral 
condensates. In fact, as the chiral condensate will, in large fields, point in 
the charged vr direction, the baryonic states will not have a definite charge. 
Whether the proton-neutron reversal takes place for fields below those that 
are screened by chiral condensates or vice versa is a subtle question; the 
approximations used in this paper and in Ref. [|218|| are not reliable to give 



an unambiguous answer. The treatment of the effects of magnetic fields on 
the strong force contributions to the baryon masses relies on the Nambu- 
Jona-Lasinio model and in Ref. [ 218 1 the variation of f^^ with magnetic field 



was not taken into account. It is clear from this discussion and the one from 
the previous paragraph that we cannot push the results of this calculation 
past few xlO^^ T. 

Experimental Consequences 

The mass evolution of protons, neutrons and electrons in magnetic fields and 
due to electromagnetism alone, will force a proton to decay in a very intense 
field. Including the effects of chiral condensates diminishes the field even 
further. Qualitatively, it is clear that the effect enhances the electromagnetic 
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contribution but its exact value depends on the model. This points to a novel 
astrophysical mechanism for creation of extra galactic positrons. 

5.1.8 Proton neutron mass difference by a lattice cal- 
culation 

One can take another approach, in principle exact, to calculate the mass 
difference using lattice gauge theory. We introduce the magnetic field in the 
lattice by multiplying a link U^{x) with a phase U^{x). Fixing the field, for 
convenience, in the direction z we then set: 

U^{x) = exp{—ieBa?yLx) for x = L^, — 1; = otherwise (5.35) 

Uy{x) = exp{ieBa'^{x — xq) (5.36) 

where xq is an offset for the magnetic field. 

Consequently the plaquette in the x-y plane is 

P^{x) = exp{ieBa^{-L^Ly + 1)) for x = L^ - l,y = Ly - 1 (5.37) 

P^{x) = exp{ieBa^) otherwise (5.38) 

The magnetic field is homogeneous only if B is is quantized as o^eB = 
2'Kn{LxLy)~^ This is a troublesome condition since the field is very large for 
reasonable lattice size. In the simulations we ignored this condition inducing 
some inhomogeneity. The results are very preliminary: at this stage it is 
difficult to find these effects with present lattice technology, but in principle 
it is possible. For details of the calculations see [^25||. This method is still 



not efficient given the present state of the art in simulations. 

5.2 Screening of Very Intense Magnetic Fields 
by chiral symmetry breaking 

Now we discuss another interesting phenomenon if very strong fields could 
be created. 

In very intense magnetic fields, B > 1.5 x 10^^ G, the breaking of the 
strong interaction SU{2) x SU{2) symmetry arranges itself so that instead 
of the neutral a field acquiring a vacuum expectation value it is the charged 
vr field that does and the magnetic field is screened. 
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In the previous section we discussed that fields with comphcated inter- 
actions of non-electromagnetic origin can induce various instabilities in the 
presence of very intense magnetic fields. By very intense we mean 10^^ G 



to 10 G. Fields with anomalous magnetic moments |224| or fields coupled 



by transition moments [|217|| may induce vacuum instabilities. The usual 
breaking of the strong interactions, chiral symmetry [xSM) , is incompatible 
with very intense magnetic fields. Using the standard SU{2) x SU{2) chiral a 
model we show that magnetic fields B > Be with Be = \/2m.„fT^ are screened; 
/^ = 132 MeV is the pion decay constant and m^ is the mass of the charged 
pions. This result is opposite to what occurs in a superconductor; in that 
case it is weak fields that are screened and large ones penetrate and destroy 
the superconducting state. 

As the magnetic fields are going to be screened we must be very careful 
in how we specify an external field. One way would be to give /j^ a spatial 
dependence and take it to vanish outside some large region of space. In the 
region that /^r vanishes we could specify the external field and see how it 
behaves in that part of space where chiral symmetry is broken. This is the 
procedure used in studying the behavior of fields inside superconductors. In 
the present situation we find this division artificial and, instead of specifying 
the magnetic fields, we shall specify the external currents. Specifically we 
will look, at first, at the electromagnetic field coupled to the charged part of 
the a model and to the current J in a long straight wire. From this result it 
will be easy to deduce the behavior in other current configurations. We will 
discuss a solenoidal current configuration towards the end of this work. 

The Hamiltonian density for this problem is 

H = -Vo- • Vcr + -Vtto • Vtto + (V + eA)7r'^ • (V - eA)n 

+ ^(^2 + ^ . ^ _ j2)2 ^ ^2 ^^^ _ or) + 1(V X A)2 - j . A ;(5.39) 

j is the external current. We have used cylindrical coordinates with p the 
two dimensional vector normal to the z direction. We will study this problem 
in the limit of very large g, where the radial degree of freedom of the chiral 
field is frozen out and we may write 



a 


= IttCOSX, 


TTO 


= U sin X COS 6^ , 


TTx 


= /vr sin X sin 6 cos 
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% = fn sin X sin 6^ sin . (5.40) 

In terms of these variables the Hamiltonian density becomes 

H = f ( Vx)' + ^ sin^ xi^O)' + Y sin^ x sin^ e{V<P - eAf 
+ m2/2(l-cosx) + ^(Vx A)2-j-A. (5.41) 

The angular field can be eliminated by a gauge transformation. For a 
current along a long wire we have 

j = /5(p)z; (5.42) 

The vector potential will point along the z direction, A = Az and the fields 
will depend on the radial coordinate only. The equations of motion become 

— V^x + sin X cos xC^^Y + c^ sin x cos x sin^ OA^ + m^ sin x = , 

V(sin^xV6') + e^sin^Xsin6'cos6'A^ = 0, 
-V^A + e^fl sin^ x sin^ 9 A - 15{p) = (16.43) 

In the absence of the chiral field the last of the Eqs. (|5.43| ) gives the classical 
vector potential due to a long wire 

A = — In^, (5.44) 

with a an ultraviolet cutoff. The energy per unit length in the z direction 
associated with this configuration is 

7-2 p 

E = —\n-, (5.45) 

Air a 

where R is the transverse extent of space (an infrared cutoff). 

Before discussing the solutions of ([5. 431) it is instructive to look at the case 



where there is no explicit chiral symmetry breaking, m^r = 0. The solution 
that eliminates the infrared divergence in Eq. (|5.45|) is x = ^ = ^/2 and A 
satisfying 

- VM + e^f^A - U{p) = . (5.46) 

For any current the field A is damped for distances p > l/ef-,^ and there is no 
infrared divergence in the energy. (Aside from the fact that chiral symmetry 
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is broken explicitly, the reason the above discussion is only of pedagogical 
value is that the coupling of the pions to the quantized electromagnetic field 
does break the SU{2) x SU{2) symmetry into SU{2) x U{1) and the charged 
pions get a light mass, m,r ~ 35 MeV ||194|| , even in the otherwise chiral 
symmetry limit.) 



The term in Eq. ( 5.41 ) responsible for the pion mass prevents us from 
setting X = 7'"/2 everywhere; the energy density would behave as vr/^m^i?^, 
an infrared divergence worse than that due to the wire with no chiral field 
present. We expect that x will vary from 7r/2 to as p increases and that 
asymptotically we will recover classical electrodynamics. Although we cannot 
obtain a closed solution to Eqs. ( ^.43 ), if the transition between x = ^/2 
and X = occurs at large p, we can find an approximate solution. The 
approximation consists of neglecting the (Vx)^ term in Eq. ( |5.41|) ; we shall 
return to this shortly. The solution of these approximate equations of motion 
is 



X 

e 

A 



2 





for p < po 
for p> po 



-2^ [Ko{ef^p) - IoieUp)Ko{ef^po)/IoieUpo)] for p < 
77- In -^ for p > 



S.47) 



Po is a parameter to be determined by minimizing the energy density of Eq. 



(|5.41| ). Note that for p > po the vector potential as well as the field return to 
values these would have in the absence of any chiral fields and that for p < po 
the magnetic field decreases exponentially as -B ~ exp(— c/ttp). The physical 
picture is that, as in a superconductor, near p = a cylindrical current sheet 
is set up that opposes the current in the wire and there is a return current 
near p = po; Ampere's law insures that the field at large distances is as 
discussed above. The energy density for the above configuration, neglecting 
the spatial variation of x? is 



H 



4:11 



KojeUpo) 
loieUpo) 



+ ln(e/^po) 



+ nmif^p' + 



(5.48) 



where the dots represent infrared and ultraviolet regulated terms which are, 
however, independent of po. For po > l/e/,r the term involving the Bessel 
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functions may be neglected and minimizing the rest with respect to po yields 

This is the main result of this work. 

We still have to discuss the validity of the two approximations we have 
made. The neglect of the Bessel functions in Eq. ( ^.48|) is valid for ef-^pQ > 1 



which in turn provides a condition on the current /, el/rriT, > 2\/2ti or more 
generally 

I/m^ » 1 . (5.50) 

The same condition permits us to neglect the spatial variation of x around 
p = Po- Let X vary from 7r/2 to in the region p — d/2 to p + d/2, with 1/d 
of the order of /^ or ttt,^ . The contribution of the variation of x to the energy 
density is AH = n^f^pod. Eq. ( ^.50| ) insures that AH is smaller than the 
other terms in Eq. ( |5.48| ). 

Eq.( |5.49|) has a very straightforward explanation. It results from a com- 
petition of the magnetic energy density ^B"^ and the energy density of the 
pion mass term m^/^(l — cosx)- The magnetic field due to the current / is 
B = I/2np and the transition occurs ai B = B^ with B^. = \/2m^/^. The 
reader may worry that the magnetic fields very close to such thin wires are so 
large as to invalidate completely the use of the chiral model as a low energy 
effective QCD theory. In order to avoid this problem we may consider the 
field due to a solenoid of radius R. The field is zero outside the solenoid, 
B = B{p)z inside with B{R) = Bq. At no point does the field become un- 
boundedly large. Using the same approximations as previously we obtain the 
following solutions of the equations of motion (for Bq > B^) 

( ioi p> R 
X = < f for i? > p > Po 
[ for p < Po , 

e =-, 

2 ' 

( ior p> R 

B = < aiKoief^rp) + a2/o(e/,rP) for i? > p > po (5.51) 

[ Be for p < Po ; 

Continuity of the vector potential determines the coefficients ai and 02, 

ai = [Boh{eUpo)-BMeUR)]/D{R,po), 

as = [BoKi{eUR)-BeKi{eUpo)]/D{R,po); (5.52) 
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D{R,po) = Ki{ef^R)Ii{eUpo) - i^i(e/^po)A(e/^i?) and po is determined, 
once more, by minimizing the energy. For {R, po) > ^/^fn 

Po = R-^\n (Bo/B,) . (5.53) 

For Bq < Be, X = everywhere and B{p) = Bq in the interior of the solenoid. 
Thus, for any current configuration, the chiral fields will adjust themselves 
to screen out fields larger than B^. Topological excitations may occur in the 
form of magnetic vortices; the angular field of Eq. (|5.4CI|) will wind around 
a quantized flux tube of radius l/ef^ [|172 |. 



5.3 The effect of strong magnetic fields on 
the elect roweak vacuum 



It was pointed out by Ambjorn and Olesen ||224|| (see also Ref. ||226|| ) that 



the Weinberg-Salam model of electroweak interactions shows an instability 
at i? ~ 10^*^ Gauss. The nature of such instability can be understood by 
looking at the expression of the energy of a particle with electric charge e, 
and spin s, moving in homogeneous magnetic field B directed along the z- 



axis. As we already discussed in Sec.|3.1|, above a critical field B^ = m /e 



particle energy is discretized into Landau levels 

El = kl + {2n + l)e|B| - 2eB • s + m^ . (5.54) 

We observe that energy of scalar [s = 0) and spinor {sz = ±1/2) is always 
positive, and indeed no instability arise in QED (it is possible to verify that 
quantum one-loop corrections do not spoil this conclusion). In the case of 
vector particles (s^ = 0, ±1), however, the lowest energy level (n = 0, k^ = 
0, Sz = +1) becomes imaginary for B > Be, which could be the signal 
of vacuum instability. The persistence of imaginary values of the one loop 
corrected lowest level energy [ p24|| seems to confirm the physical reality of 
the instability. 

As it is well known the Weinberg-Salam model contains some charged 
vector fields, namely the W^ gauge bosons. The coupling of the W^ field to 
an external electromagnetic field A^f* is given by 

An* = -l\F'^:f ^iD.W.-D^W.l'-ml^WlW^-zeF'^^'W^W (5.55) 
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with 

D^ = d,-zeAf . (5.56) 

The important term in the previous expression is the "anomalous" magnetic 
moment term ieFf^^^W^W^ , which arises because the non-Abehan nature of 
the SU{2) component of the Weinberg-Salam model gauge group structure. 
Due to this term the mass eigenvalues of the W Lagrangian becomes 

w? = m^ ± eB . (5.57) 

As expected from the considerations in the above, a tachyonic mode ap- 
pears for B > Be- The corresponding eigenvector for zero kinetic energy is 
determined by solving the equation of motions 

DiW, - D,Wi = z, J = x, y , (5.58) 

where Wi^2 = Wx^iWy. Ambjorn and Olesen argued that a suitable solution 
of this equation is 

|iy(x,y)| = e-3'"^(^''+?'') , (5.59) 

corresponding to a vortex configuration where VT-fields wind around the z- 
axis. This configuration corresponds to the Nielsen-Olesen vortex solution 
||158|| . A similar phenomenon should also take place for Z bosons. Given the 
linearity of the equations of motion it is natural to assume that a superpo- 
sitions of vortices is formed above the critical field. This effect resemble the 
behaviour of a type-II superconductor in the presence of a critical field mag- 
netic field. In that case U{1) symmetry is locally broken by the formations of 
a lattice of Abrikosov vortices in the Cooper-pairs condensate through which 
the magnetic field can flow. In the electroweak case this situation is reversed, 
with the formation oi a W condensate along the vortices. Concerning the 
back-reaction of the W condensate on the magnetic field, an interesting effect 
arises. By writing the electric current induced by the W fields 

j^{W) = 2ie {W^D^W - WD^W^) , (5.60) 

Ambjorn and Olesen noticed that its sign is opposite to the current induced 
by the Cooper pairs in a type-II superconductor, which is responsible for the 
Meissner magnetic field screening effect. Therefore, they concluded that the 
PF-condensate induce anti- screening of the external magnetic field. 



124 



Although the Higgs field $ does not couple directly to the electromagnetic 
field (this is difi^erent from the case of a superconductor where the Cooper- 
pairs condensate couples directly to A^^*) , it does through the action of the 
W condensate. This can be seen by considering the Higgs, W potential in 
the presence of the magnetic fields: 

V{(f), W) = 2 (eB - mlr) |Vr|^ + /^Wr - 2A0^0^ + 2g'^\W\^ + X (^t^o) ■ 

(5.61) 
In the above 0o and 0+ are respectively the Higgs field vev and charged 
component, g is the SU(2) coupling constant, and A is the Higgs the self- 
interaction coupling constant. We see that the Vl^-condensate influences the 
the Higgs field at classical level due to the 0^|Vrp term. It is straightforward 
to verify that if eB < m^ = \g'^4>'i the minimum of 1^(0, W) sits in the stan- 
dard field value = 0o with no W condensate. Otherwise a W condensate 
is energetically favoured with the minimum of the potential sitting in 

'Pmin = 00 —2 Zir (5-62) 



my -m^y 



where 



rriH = 4A0O, m^, = -g'^cpl . (5.63) 

We see that the Higgs expectation value will vanish as the average magnetic 
field strength approaches zero, provided the Higgs mass is larger than the W 
mass. This seems to suggest that a PF-condensate should exist for 

m^ < eB < m]j , (5.64) 

and that the SU{2) x t/y(l) symmetry is restored above H'^^ = m'jj/e. 
Thus, anti-screening should produce restoration of the electroweak symme- 
try in the core of W vortices. If rriH < friw the electroweak vacuum is 
expected to behave like a type I superconductor with the formation of ho- 
mogeneous VT-condensate above the critical magnetic field. The previous 
qualitative conclusion have been confirmed by analytical and numerical com- 
putations performed for mfj = rriw in Ref . [p24| | , and for arbitrary Higgs mass 
inRefs.p^Pl. 

A different scenario seems, however, to arise if thermal corrections are 
taken into account. Indeed, recent finite temperature lattice computations 



||188|| showed no evidence of the Ambjorn and Olesen phase. According to 
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Skalozub and Demchik [2^] such a behaviour may be explained by properly 



accounting the contribution of Higgs and gauge bosons daisy diagrams to the 
effective finite temperature potential. 

In conclusion, it is quite uncertain if the Ambjorn and Olesen phenomenon 
was really possible in the early Universe. 

5.3.1 The electroweak phase transition in a magnetic 
field 

We shall now consider the possible effects of strong magnetic fields on the 
electroweak phase transition (EWPT). As it is well known, the properties 
of the EWPT are determined by the Higgs field effective potential. In the 
framework of the minimal standard model (MSM), taking into account ra- 
diative corrections from all the known particles and for finite temperature 
effects, one obtains that 

Kff(0,T) ~ -^iii^ - aT^)<P^ - T5cf + ^(A - 5At)0' . (5.65) 

where is the radial component of the Higgs field and T is the temperature 



(for the definitions of the coefficients see e.g. Ref. [[L87]] ). 

A strong hypermagnetic field can produce corrections to the effective 
potential as it affects the charge particles propagators (see below). There 
is, however, a more direct and simpler effect of magnetic and hypermag- 
netic fields on the EWPT which was recently pointed-out by Giovannini and 



Shaposhnikov | |179|| and by Elmfors, Enqvist and Kainulainen ||187|| . The au 



thors of Refs. [|179| , |187|| noticed that hypermagnetic fields affect the Gibbs 



free energy (in practice the pressure) difference between the broken and the 
unbroken phase, hence the strength of the transition. The effect can be un- 
derstood by the analogy with the Meissner effect, i.e. the expulsion of the 
magnetic field from superconductors as consequence of photon getting an ef- 
fective mass inside the specimen. In our case, it is the Z-component of the 
hypercharge f/(l)y magnetic field which is expelled from the broken phase. 
This is just because Z-bosons are massive in that phase. Such a process has 
a cost in terms of free energy. Since in the broken phase the hypercharge 
field decomposes into 

A^ = cos6'^v4^ - sin6'^Z^ , (5.66) 
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we see that the Gibbs free energy in the broken and unbroken phases are 



Gu 



V{cl>)--cos^e.^{B\. 



1 

2 



V{id)-UB^^'f 



(5.67) 
(5.68) 



where i?f?'* is the external hypermagnetic field. In other words, compared 

to the case in which no magnetic field is present, the energy barrier between 

unbroken and broken phase, hence the strength of the transition, is enhanced 

X 
by the quantity - sin^ OwiBy 
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According to the authors of Refs. ||179 , 

this effect can have important consequence for baryogenesis. 

In any scenario of baryogenesis it is crucial to know at which epoch do 
the sphaleronic transitions, which violate the sum {B + L) of the baryon and 
lepton numbers, fall out of thermal equilibrium. Generally this happens at 
temperatures below T such that 



E{f) 
T 



> A 



(5.69) 



where E{T) is the sphaleron energy at the temperature T and A c^ 35 — 45, 
depending on the poorly known prefactor of the sphaleron rate. In the case 
of baryogenesis at the electroweak scale one requires the sphalerons to drop 
out of thermal equilibrium soon after the electroweak phase transition. It 
follows that the requirement T = T^, where T^ is the critical temperature, 
turns eq. ( p.69| ) into a lower bound on the Higgs vacuum expectation value 
(VEV), 

^>1. (5.70) 



As we already discussed, it is by now agreed ||144|| that the standard model 
(SM) does not have a phase transition strong enough as to fulfill Eq.( ^.70|) , 
whereas there is still some room left in the parameter space of the minimal 
super symmetric standard model (MSSM) ||145 |. 

The interesting observation made in Refs. ||179| , |187|| is that a magnetic field 
for the hypercharge U{1)y present for T > T^. may help to fulfill Eq.( p.70 ). 
In fact, it follows from the Eqs.( |5.67D , that in presence of the magnetic field 
the critical temperature is defined by the expression 



\/(O,T,)-\/(0,T, 



-sin^e^{Br{Tc)f 



(5.71) 
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This expression implies a smaller value of T^ than that it would take in the 
absence of the magnetic field, hence a larger value of the ratio ( |5.70| ). 



Two major problems, however, bar the way of this intriguing scenario. 
The first problem is that by affecting fermion, Higgs and gauge field prop- 
agators, the hypermagnetic field changes the electroweak effective potential 
in a nontrivial way. Two different approaches have been used to estimate 



the relevance of this kind of effects based either on lattice simulations [144 



or analytical computations [ p29|| . Both approaches agreed in the conclusion 
that for a Higgs field mass compatible with the experimental constraints 
{niH > 75 GeV), and for field strengths B, By < 10^^ G, the standard 
model EWPT is second order or a cross-over. Although this negative result 
could, perhaps, be overcome by adopting a supersymmetrical extension of 
the standard model (see e.g. Ref. ||145| |), a second, and more serious problem 



arises by considering the effect of the magnetic field on the anomalous pro- 
cesses (sphalerons) which are responsible for lepton and baryon violation at 
the weak scale. This effect will be the subject of the next section. 

5.3.2 Sphalerons in strong magnetic fields 

The sphaleron, is a static and unstable solution of the field equations of the 
electroweak model, corresponding to the top of the energy barrier between 
two topologically distinct vacua ||231|| . In the limit of vanishing Weinberg an- 



gle, 9^u — > 0, the sphaleron is a spherically symmetric, hedgehog-like configu- 
ration of SU{2) gauge and Higgs fields. No direct coupling of the sphaleron 
to a magnetic field is present in this case. As 9^ is turned on, the f/y(l) field 
is excited and the spherical symmetry is reduced to an axial symmetry. A 
very good approximation to the exact solution is obtained using the Ansatz 



by Klinkhamer and Laterveer ||232|| , which requires four scalar functions of r 
only, 

g'a,dx'={l~fo{0)F,, 

gWta'' dx' = (1 - fimFi^' + ^^2^') + (1 - f3(0)F3<j' , 

^^U m ) ' ^'-''^ 

where g and g' are the SU{2)l and U{1)y gauge couplings, v is the Higgs 
VEV such that Mw = gv/2, Mu = V2Xv, ^ = gvr, a" (a = 1, 2, 3) are the 
Pauli matrices, and the -F^'s are 1-forms defined in Ref. ||231|| . The boundary 
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conditions for the four scalar functions are 



/(O , fsiO , HO - /o(0 - 
/(O , fsiO , HO , /o(0 



for ,^ 
for ^ 





oo - 



(5.73) 



It is known [p31| , p32|] that for 6*^^, 7^ the sphaleron has some interesting 
electromagnetic properties. In fact, differently from the pure SU{2) case, in 
the physical case a nonvanishing hypercharge current Jj comes-in. At the 
first order in 9^, Ji takes the form 



J, 



(1) 



j.h\o[i-fm 



-q V 



^3ij-^^ 



3 ' 



(5.74) 



where h and / are the solutions in the ft 

moment 

,(1) _ ^_^ 



limit, giving for the dipole 



/i^ 



3 g^v 
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/ d^enmi-fio]- 

Jo 



(5.75) 



The reader should note that the dipole moment is a true electromagnetic 
one because in the broken phase only the electromagnetic component of the 
hypercharge field survives at long distances. 



Comelli et al. ||189|| considered what happens to the sphaleron when an 
external hypercharge magnetic field, i?f?*, is turned on. They found that the 
energy functional is modified as 



with 



and 



^n 



d^x 



E — Eq — Edip 



ii^^.i^^. + ^^^- + (A^)nA$) 



Edip = J d^xJiAj = -j d^xfijfi 



y($) 



V 



(5.76) 



(5.77) 



(5.78) 



ext 



with fij = diAj — djAf . A constant external hypermagnetic field By 
directed along the X3 axis was assumed. In the 9^ ^ limit the sphaleron 
has no hypercharge contribution and then E^^^ = 0. At 0(9^), using (|5.74|) 
and ( |5.75| ) the authors of Ref. |p.89|| got a simple magnetic dipole interaction 
energy 



771(1) _ ..Wnext 



(5.79) 
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Figure 5.1: The VEV at the critical temperature, v{Tc), and the sphaleron 
energy vs. the external magnetic field for M^ = M-^. We see that even if 
v{Tc)/Tc > 1 the washout condition E/Tc > 35 is far from being fulfilled. 
FromRef.p9| 



In order to assess the range of validity of the approximation ( |5.79|) one needs 
to go beyond the leading order in 6^ and look for a nonlinear i?y^*-dependence 
of E. This requires to solve the full set of equations of motion for the gauge 
fields and the Higgs in the presence of the external magnetic field. Fortu- 
nately, a uniform B^^ does not spoil the axial symmetry of the problem. 
Furthermore, the equation of motion are left unchanged (i9j/?^* = 0) with 
respect to the free field case. The only modification induced by By^ resides 
in the boundary conditions since - as i^ ^ oo - we now have 



/(0,M0 



1 



/3(0,/o(0 



1 



5f^*sin2^^ 



e 



(5.80) 
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whereas the boundary condition for ^ ^ are left unchanged. 

The solution of the sphaleron equation of motions with the boundary con 
ditions in the above were determined numerically by the authors of Ref. ||189 



They showed that in the considered i?f?'*-range the corrections to the linear 
approximation 

AE-Zi^cos^^Efi^* 

are less than 5%. For larger values of -By^* non-linear effects increase sharply. 
However, as we discussed in the previous sections, for such large magnetic 
fields the broken phase of the SM is believed to become unstable to the 
formation either of ly-condensates |p24|| or of a mixed phase [p.88|| . In such 



situations the sphaleron solution does not exist any more. Therefore, it is 
safe to limit the previous analysis to values By^ < 0.4 T^. 

The reduction of sphaleron energy due to the interaction with the field 
By has relevant consequences on the sphaleronic transition rate which is 
increased with respect to free field case. As a consequence, in an external 
magnetic field the relation between the Higgs VEV and the sphaleron energy 
is altered and Eq. (|5.70|) does not imply (|5.69|) any more. We can understand 



it by considering the linear approximation to E, 

E ^ E{B^y = 0) -/i«i?r COS0. ^^(eo- '^^^^^m'^A (5.81) 

^ V 9 v^ J 

where m^^^ is the 0{9w) dipole moment expressed in units of e/awMw{T). 
From the Fig. |5.1| we see that even if v{Tc)/Tc > 1 the washout condition 
E/Tc > 35 is far from being fulfilled. 

It follows form the previous considerations, that even if a strong magnetic 
fields might increase the strength of the EWPT, such an effect would not help 
baryogenesis. 
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Conclusions 



In this review we have analyzed a large variety of aspects of magnetic fields 
in the early Universe. Our exposition followed an inverse-chronological order. 
In the first part of Chap .|T] we discussed what observations tell us about recent 
time fields and their evolution in galaxies and clusters of galaxies. As we have 
seen, a final answer about the origin of these fields is not yet available. Several 
arguments, however, suggest that galactic and cluster fields were preexisting, 
or at least contemporary to, their hosts. The main reasons in favor of this 
thesis are: the ubiquity of the fields and the uniformity of their strength; the 
theoretical problems with the MHD amplification mechanisms, especially to 
explain the origin of cluster fields; the observation of fiG magnetic fields in 
high-redshift galaxies. It is reassuring that new ideas continuously appear to 
determine these fields at all times. For example, a very recent one by Loeb 
and Waxman proposes to look for fluctuations in the radio background from 
intergalactic synchrotron emission of relativistic electrons interacting with 
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A consistent and economical mechanism which may naturally explain the 
early origin of both galactic and cluster magnetic fields, is the adiabatic 
compression of a primeval field with strength in the range Bq ~ 10^^ — 10^^° 
Gauss (-Bo is the intensity that the primordial field would have today under 
the assumption of adiabatic decay of the field due to the Hubble expansion). 
If this was the case, two other interesting effects may arise: a) magnetic 
fields may have affected structure formation perhaps helping to solve some 
of the problems of the CDM scenario; b) magnetic fields can have produced 
observable imprints in the CMBR. Given the current theoretical uncertainties 
about the MHD of galaxies and clusters, and the preliminary status of N-body 
simulations in the presence of magnetic fields, the most promising possibility 
to test the primordial origin hypothesis of cosmic magnetic fields comes from 
the forthcoming observations of the CMBR anisotropies. 
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In Chap.0 we reviewed several possible effects of magnetic fields on the 
CMBR. In first place, we showed that magnetic fields may affect the isotropy 
of the CMBR. A magnetic field which is homogeneous through the en- 
tire Hubble volume, would spoil Universe isotropy giving rise to a dipole 
anisotropy in the CMBR. On the basis of this argument it was shown that 
COBE measurements provide an upper limit on the present time equivalent 
strength of a homogeneous cosmic magnetic field which is roughly 3 x 10~^ G. 
More plausibly, magnetic fields are tangled on scales much smaller than the 
Hubble radius. In this case the effect on the Universe geometry is negligible 
and much more interesting effects may be produced on small angular scales. 
Some of these effects arise as a consequence of MHD modes appearing in 
the magnetized photon-baryon plasma in place of the usual acoustic modes. 
The amplitude and velocity of the MHD modes depend on the magnetic 
field intensity and spatial direction. Some of these modes are quite different 
from standard scalar and tensor modes which are usually considered in the 
theoretical analysis of the CMBR distortions. For example, Alfven waves 
have the peculiar property of not beeing depleted by the Universe expan- 
sion inspite of their vectorial nature. These modes are well suited to probe 
perturbations as those generated by cosmic defects and primordial phase 
transitions. Tangled magnetic fields, whose production is predicted by sev- 
eral models and, which are observed in the intercluster medium, are also 
expected to produce Alfven waves. Another interesting aspect of this kind 
of isocurvature perturbations, is that they are not affected by Silk damping. 
The polarization power spectrum of CMBR can also be affected by primor- 
dial magnetic fields. This a consequence of the Faraday rotation produced 
by the field on the CMB photons on their way through the last scattering 
surface. Magnetic fields with strength Bq > 10^^ G may have produced a 
detectable level of depolarization. Furthermore, it was shown that because of 
the polarization dependence of the Compton scattering, the depolarization 
can feed-back into a temperature anisotropy. It was concluded that the best 
strategy to identify the imprint of primordial magnetic fields on the CMBR 
is probably to look for their signature in the temperature and polarization 
anisotropics cross-correlation. This method may probably reach a sensitivity 
of ~ 10~^° G for the present time equivalent magnetic field strength when 
the forthcoming balloon and satellite missions data is analyzed. Some results 



from the Boomerang and Maxima experiments [jT^ are already out with sur 



prising results. It is too early to decide the reasons why, if experimentally 



confirmed, the second peak is low. We have verified, together with Edsjo |65| 
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that magnetic fields can only decrease the ratio of amplitude of the second 
peak with respect to the first. Therefore the effect cannot be explained in 
terms of primordial magnetic fields. Interesting constraints on the strength 
of these fields will be available only when the amplitude of several peaks and 
the polarization of CMBR will be measured. 

Another period of the Universe history when primordial magnetic fields 
may have produced observable consequences is the big-bang nucleosynthesis. 
This subject was treated in Chap.|^. Three main effects have been discussed: 
the effect of the magnetic field energy density on the Universe expansion; the 
modification produced by a strong magnetic field of the electron-positron 
gas thermodynamics; the modification of the weak processes keeping neu- 
trons and protons in chemical equilibrium. All these effects produces a vari- 
ation in the final neutron-to-proton ratio, hence in the relative abundances 
of light relic elements. The effect of the field on the Universe expansion rate 
was showed to be globally dominant, though the others cannot be neglected. 
Furthermore, the non-gravitational effects of the magnetic field can exceed 
that on the expansion rate in delimited regions where the magnetic field 
intensity may be larger than the Universe mean value. In this case these 
effects could have produced fluctuations in baryon to photon ratio and in the 
relic neutrino temperature. Apparently, the BBN upper bound on primordial 
magnetic flelds, which is -Bq ^ 7 x 10^^ G, looks less stringent than other 
limits which come from the Faraday rotation measurements (RMs) of dis- 
tant quasars, or from the global isotropy of the CMBR. However, we showed 



in Sec. |3.4| that this conclusion is not correct if magnetic flelds are tangled. 
The reason is that BBN probes length scales which are of the order of the 
Hubble horizon size at BBN time (which today corresponds approximately 
to 100 pc) whereas CMBR and the RMs probe much larger scales. Further- 
more, constraints derived from the analysis of effects taking place at different 
times may not be directly comparable if the magnetic fleld evolution is not 
adiabatic. 

In Chap.^ we reviewed some of the models which predict the genera- 
tion of magnetic flelds in the early Universe. We flrst discussed those models 
which invoke a flrst order phase transition. This kind of transitions naturally 
provide some conditions, as charge separation, out-of-equilibrium condition, 
and high level of turbulence, which are known to be important ingredients 
of magnetogenesis. We discussed the cases of the QCD and the electroweak 
phase transitions (EWPT). In the case of the EWPT some extension of the 
particles physics standard model has to be invoked for the transition to be 
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first order. Magnetic fields may be generated during the EWPT from the 
non-trivial dynamics of the gauge fields produced by the equilibration of the 
electrically charged components of the Higgs field. This effect resembles the 
Kibble mechanism for the formation of topological defects. It is interesting 
that such a mechanism gives rise to magnetic fields, though only on very 
small scales, even if the phase transition is second order. In general, since 
the production of magnetic fields during a phase transition is a causal phe- 
nomenon, the coherence length scale of these fields at the generation time 
cannot exceed the horizon radius at that time. Typically, once this length 
is adiabatically re-scaled to present time, one gets coherence cell sizes which 
are much smaller than those observed today in galaxies and the inter-cluster 
medium. This problem may be eased by the effect of the magnetic helicity 
which is expected to be produced during primordial phase transitions. Helicity 
may help the formation of large magnetic structures starting from small ones 
(inverse cascade). Indeed, this is a quite common phenomenon in MHD. Some 
estimates of the quantitative relevance of this effect have been given, for ex- 
ample, in Sec. |1.4| . We have seen that the QCD phase transition might indeed 
give rise to phenomenological interesting values of the present time magnetic 
field strength and coherence size but only assuming quite optimistic condi- 
tions. Magnetic fields produced at the EWPT might have played a role in 
the generation of galactic magnetic fields only if they were amplified by a 
galactic dynamo. The problem with the small coherence scale of magnetic 
fields produced in the early Universe may be circumvented if the produc- 
tion mechanism was not-causal. This may be possible if magnetic fields were 
produced during infiation by the superadiabatic amplification of preexisting 
quantum fiuctuations of the gauge fields. This phenomenon, however, can 
take place only if the conformal invariance of the electromagnetic field is 
broken. In Sec.E^ we have discussed several interesting mechanisms which 



have been proposed in the literature to avoid this obstacle. Unfortunately, 
although some results are encouraging, at the present status of art, none of 
these model seems to offer any firm prediction. Further work on the subject 
is, therefore, necessary. 

Even if magnetic fields, produced during the electroweak phase transitions 
or before, are nor the progenitor of galactic magnetic fields, they may still 
have had other interesting cosmological consequences. Perhaps the most 
intriguing possibility is that magnetic fields played a role in the generation 
of the baryon asymmetry of the Universe (BAU). The magnetic fields may 
influence electroweak baryogenesis at two levels. At a first level, magnetic 

135 



fields can play an indirect role on electroweak baryogenesis by modifying 
the free energy difference between the symmetric and broken phases, the 
Higgs effective potential, and the rate of sphaleron, baryon number violating, 
transitions. In Chap.|^ we showed, however, at this level, no significative 
modifications to arise with respect to the standard scenario. Magnetic fields, 
or better their hypermagnetic progenitors, may have played a much more 
direct role in the generation of the BAU if they possessed a net helicity. 
Indeed, it is well known from field theory that the hypermagnetic helicity 
coincides with the Chern-Simon number which can be converted into baryons 
and leptons numbers by the Abelian anomaly. The origin of the primordial 
magnetic helicity is still matter of speculation. Among other possibilities 
which we have review in Chap.^ one of the more discussed in the literature 
is that a net hypermagnetic helicity of the Universe arise by an anomalous 
coupling of the gauge fields to an oscillating pseudoscalar field. The existence 
of pseudoscalar fields of this kind is required by several extensions of the 
particle physics standard model. However, it must be admitted that the 
mechanisms for generation of fields that are large and extended at the same 
time are far from being fully understood. 

Large magnetic fields would also have a profound effect on chirality but 
it is also quite far from observability. As a rule of the thumb, particle physics 
effects will appear at the earliest at B = m"^, the tt being the lightest hadron. 
This makes these effects difficult to test. Large fields are required. The only 
hope is the existence of superconductive strings. 

Are the observed fields, so widespread, of early origin or some seeds fields 
were rapidly enhanced by a dynamo mechanism? This question remains 
unanswered but high precision CMBR acoustic peak measurements may very 
well provide a breakthrough. 
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